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Les méthodes qui se rattachent aux transformations intégrales 
trouvent un vaste champ d'applications en analyse mathématique 
ces derniers temps. Elles ont été utilisées avec succès dans la résolu- 
tion d'équations différentielles et intégrales, dans l’étude des fonc- 
tions spéciales et le calcul d’intégrales. L'avantage de ces méthodes 
est dans le fait qu’elles permettent de construire les tables des 
transformations directes et inverses des diverses fonctions que l’on 
rencontre dans les applications. 

Dans cet ouvrage on considère les transformations intégrales 
les plus répandues. Une première partie de cinq chapitres est con- 
sacrée aux principes de la théorie. Le premier chapitre expose des 
éléments de théorie des transformations de Fourier et certaines de 
leurs applications. Le second chapitre qui est le plus important et 
le plus long traite de la transformation de Laplace, ainsi que de celle 
de Mellin. L'objet du troisième chapitre est la transformation inté- 
grale de Bessel, ou plus exactement les transformations intégrales 
ayant une fonction de Bessel pour noyau: les transformations de 
Hankel, Meijer, Kontorovitch-Lébédev. Dans le chapitre quatre les 
auteurs passent brièvement en revue d’autres transformations inté- 
grales. Le cinquième est réservé aux principes de la théorie du calcul 
opérationnel. 

Le calcul symbolique ou opérationnel, on le sait, fut systémati- 
quement élaboré dans la moitié du siècle passé. A la fin du XIX!, 
Heaviside l’applique avec bonheur à la résolution de certains pro- 
blèmes se rattachant à la théorie des oscillations électromagnétiques. 
Le large développement du calcul symbolique est à l’origine de 
nombreux travaux dont l’objet était précisément de le justifier. 
Le point de vue opérationnel de Heaviside fut évincé par les travaux 
de Carson, Doetch, Van der Pol, etc., qui se servirent de la transfor- 
mation de Laplace et de l’intégrale de Mellin. 

Cependant une telle situation ne pouvait pas durer trop long- 
temps, car l’essor de l’analyse fonctionnelle et, notamment, la 
théorie des opérateurs linéaires contribuèrent au développement des 


\ 
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méthodes opérationnelles en analyse mathématique. Dans les tra- 
vaux [44], [195], on trouvera un exposé opérationnel du calcul sym- 
bolique basé sur la transformation de Laplace. 

Le retour au point de vue opérationnel initial est l’œuvre de 
Mikusinski [165]. Ce dernier donne en effet une justification opéra- 
tionnelle rigoureuse du calcul symbolique sans faire usage de la 
transformation de Laplace. A ces fins il introduit diverses notations 
pour la fonction et sa valeur en un point. Il désigne la fonction par 
{f (t)} et sa valeur en un point t par f (t). Exemple: 2 est un nombre, 
{2} est une fonction prenant la valeur constante 2. 

Au chapitre V, contrairement à Mikusinski, on définit la convolu- 
tion de sorte à ne pas faire de distinction entre une constante et la 
fonction d’une constante. Comme dans nombre de cas l’usage de 
l'intégrale de Laplace simplifie singulièrement les divers calculs 
et transformations relatifs à l'établissement des formules opération- 
nelles, ici on dégage le lien existant entre le calcul opérationnel 
construit et la transformation de Laplace. Dans ce chapitre, on 
étudie également la transformation généralisée de Laplace et on 
énonce ses principales propriétés. En fin de chapitre on examine 
succinctement le calcul symbolique de l'opérateur de Bessel et l’on 
établit son lien avec la transformation de Meijer. 

La deuxième partie du livre est occupée par les tables des formu- 
les des transformations intégrales. En résolvant des problèmes con- 
crets, on établit des formules des transformations intégrales qui 
peuvent être généralisées à d’autres problèmes. Aussi les tables des 
formules des transformations intégrales possèdent-elles un vaste 
terrain d'applications qui englobe les domaines les plus variés 
des connaissances: mathématiques, physique, mécanique, électro- 
technique, etc. Les tables des formules sont précédées des notations 
des fonctions spéciales et de certaines constantes (chapitre VI). Dans 
les autres chapitres on étudie les transformations de Fourier, Laplace- 
Carson, Mellin, Hankel, Meijer, Kontorovitch-Lébédev, Mehler- 
Fock, Hilbert, etc. Pour dresser ces tables on s’est servi essentielle- 
ment de travaux portant sur le même sujet, parmi lesquels il importe 
de faire une mention spéciale à: Erdelyi A., Magnus W., 
Oberhettinger F., Tricomi F.G., Tables of integral 
Transforms, 1954; Oberhettinger F., Tabellen zur Fourier- 
Transformation, 1957. 

11 n’est pas exclu que des erreurs se soient glissées au cours de la 
manipulation d’une telle quantité de formules, aussi les auteurs 
seraient-ils reconnaissants aux lecteurs de bien vouloir leur commu- 
uiquer les incorrections relevées. 

V. Ditkine, A. Proudnikov 


PRINCIPES DE LA THÉORIE 


CHAPITRE PREMIER 


TRANSFORMATION DE FOURIER 


$ 1. Prologue à la théorie des séries 
de Fourier 


Toute fonction f ({) peut être représentée sous des hypothèses 
assez générales par une série infinie de la forme 


fO=R+ D. (an cos nt +basinnt), (1.1) 
n=1 
où 
an = + | f(t)cosntdt (n—0,1,2,...), (1.2) 
bn = — JfGsinntdt (n—1,9,...). (1.3) 


Cette série s’appelle série trigonométrique de Fourier, les nombres a, 
et b, coefficients de Fourier de la fonction f (t). Les membres de la 
série (1.1) étant tous 2x-périodiques, on peut limiter l'étude de cette 
série à tout intervalle de longueur 21. Dans le cas d’un intervalle 
de longueur arbitraire 2/, la série de Fourier 


[= + > (an cos TE +6, sin), (1.4) 


n=1 


où 


l 
an = + | F0 cos TE dt (n=0,1,2,...) (1.5) 
CL | 
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et 


L 
bn =+ | sin + dt (n=1,2,...), (1.6) 
1 


est une fonction 2/-périodique. 
Compte tenu de l'identité d'Euler 


eiF—cosp+isino, (1.7) 
on peut mettre la série (1.1) sous la forme complexe 


f (t) = Ÿ Cneint, (1.8) 


n= — 00 


x 
Cn = | f@e-tt dt (n=0, +1, +92, ...). 1(1.9) 
—7 
Parfois, au lieu de l'intervalle J—x, xl ou J0, 2xl il est plus commo- 


de de considérer un intervalle de longueur 1, par exemple, JO, 1i. 
Les coefficients de Fourier s’écrivent alors 


1 
Cn = | f(the-2nint dt (n=0, +1, +2, ...). (1.10) 
0 


Si f (ft) est une fonction paire i.e. f (—t) = f (£) et si elle est 
intégrable sur l’intervalle ]—Z, I, “ii 


oa=2 | f(#) dt. 

|!" 
De façon analogue, si f (ft) est une fonction impaire, i.e. j (—t) — 
— —f(t), alors 


l 
[f@a=0. 
l 


Si f (£) est une fonction paire sur l'intervalle ]—x, ni, la fonction 
f (£) cos nt est paire et la fonction f (t) sin nt impaire quel que soit n. 
Si l’on définit les coefficients de la série de Fourier d’une fonction 
paire f (t) d’après les formules (1.2) et (1.3), on obtient} 


I 


= | Î (£) cos nt dt (1.11) 
‘ 0 


et b, — 0. La série de Fourier de laÿ fonction j (f) ne contient que 
des cosinus ; les coefficients de cette série sont donnés par la formu- 
le (1.11). 
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Si f(t) est une fonction impaire, les fonctions f (t) cos nt et 
f (€) sin nt seront respectivement impaire et paire. La série de Fourier 
d’une fonction impaire f (f) ne contient que des sinus et ses coeffi- 
cients sont définis par 
T 
2 ; 
ta = + | (sin ni dt. (4.12) 


0 


Donc, toute fonction intégrable entre O0 et x peut formellement 
être représentée sur cet intervalle par une série de Fourier de sinus 
ou de cosinus, sans aucune hypothèse sur sa parité, son imparité, 
sa périodicité et voire même sur sa définition en dehors de cet inter- 
valle. 

Voici un théorème fondamental de la théorie des séries de Fourier. 

Théorème 1 (Riemann-Lebesgue). Si une fonc- 
tion f t) est intégrable sur un intervalle Ja, bl, alors 


b b 
j f (#) cos M dt + 0, | f(t)sin At dt+0 


lorsque À —+> co. 

Voir démonstration dans {106], [270]. Le théorème de Riemann- 
Lebessue possèdent des corollaires importants. 

Corollaire 1. Les coefficients de Fourier de toute fonction 
intégrable tendent vers zéro. 

Corollaire 2. Le comportement d'une série de Fourier en un 
point t dépend uniquement du comportement de la fonction en un voisi- 
nage immédiat de ce point (principe de‘localisation). 

Formulons le critère de convergence le plus usité d’une série de 
Fourier. Si f est une fonction à variation bornée, sa série de Fourier 


converge en tout point { vers la valeur + € ( + 0) — f (£ — 0)}. 


Si f est de plus continue sur un intervalle Ja, b{, la série de Fourier 
est uniformément convergente sur cet intervalle (cf. [270]). 


$ 2. Formule intégrale de Fourier 


Supposons qu’une fonction f (f) 2l-périodique est représentée par 
la série (1.4). En.portant dans (1.4) les expressions (1.5) de a, et (1.6) 
de b,, on obtient 


l co l 
fO=+ f(x) + À | f(r)cos Æ(r—# dr! (113) 


Si l’on pose: = = À, = — AA et que l” on passe formellement à la 


limite lorsque Z —+ co, la somme se transforme en une intégrale et 
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nous obtenons la formule intégrale de Fourier 
n OO © 
fo=+| dà | f (x) cos À (tr — +) dx. (1.14) 
Û —00 
Cette formule définit une fonction sur l’intervalle ]J—o, +ol 


de la même façon que la série de Fourier définit une fonction de 
période finie. La formule (1.14) se ramène à la forme suivante: 


1 = lim 3e L Çeun ( e-iatf (x) dx dh (145) 
-{ 00 
ou 
f (= lim + Î EUCD }; (x) dr. (1.16) 


L'intégrale du second membre de la formule (1.14) s'appelle intégrale 
double de Fourier, la formule (1.15) forme complexe de l'intégrale de 
Fourier et la formule (1.16) représentation de la fonction f ({) par 
l'intégrale simple de Fourier. Les conditions classiques de validité 
des formules précédentes s’établissent à l’aide du théorème suivante. 

Théorème 2 [32]. Supposons qu'une fonction f (t) est inté- 
grable-Lebesgue sur l'intervalle ]—o, + oo *) et qu'elle est à variation 
bornée sur tout intervalle fini. Les formules (1.14), (1.15), (1.16) sont 
valables si l’on remplace leurs premiers membres par 5 {f (t + 0) + 


+ f(t — 0)} aux points de discontinuité de f (t). Les formules (1.14), 
(1.15), (1.16) peuvent également s'écrire sous la forme suivante: 


f() = | La (u)costu+-b(u) sin tu] du, (1.17) 
0 
où 
au)=+ | f{)cosutdt, b(u)=+ | f (#) sin ut dt. 


Si f (t) est une fonction paire, la formule (1.17) s'écrit alors 


= {+ | cos tu du | j (x) cos ut dr. (1.18) 
0 0 


*) Ce qu'on note brièvement: f (ft) € L] —0o0, + ol. 
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Le développement (1.18) s'appelle développement de Fourier de f (t) 
en série de cosinus. 


De façon analogue, si f ({) est une fonction impaire, on obtient 
le développement de Fourier de f (t) en série de sinus 


[ae 


f=<+ lu sin {u du | f(x) sin ut dt. (1.19) 
0 0 


S 3. Propriétés fondamentales 
des transformées de Fourier 


1. Les formules considérées ici conduisent à des relations récipro- 
ques ou duales entre couples de fonctions. 
Si l’on pose 


En L C —iut 9 
F (u) = ETS É f{t)e dt, (1.20) 
la formule (1.15) donne 


f(t) = 


DE f F (u) eit" du, (1.21) 


où l'intégrale du second membre est prise dans sa valeur principale, 
i.e. comme la limite 


[ 
lim | F(u)eitu du. 
l— 


La fonction F (u) porte le nom de fransformée de Fourier de f (t). 

Si la fonction f ({) est intégrable sur l'intervalle ]J—oco, +oo, 
la fonction F (u) existe pour tous les t. Si les fonctions F (u) et f (f) 
sont les transformées de Fourier l’une de l’autre, on les appelle 
couple de transformées de Fourier. En posant 


F.(&)= y + Î (6) cos ut dt, (4.22) 
on déduit de (1.18) 
{(t) = y Z F4 (u) cos tu du. (1.23) 


Les fonctions ainsi liées entre elles s’appellent couple de cosinus- 
transformées de Fourier. De façon analogue, de la formule (1.19) 
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on obtient le couple de sinus-transformées de Fourier : 


F(u)= 1)" 2 f f(t)sinut dt, (1.24) 


f (#) — = | F, (u) sin tu du. (1.25) 
0 


Si f (t) est paire, 

F (u) = Fe (u); 
si f (£) est impaire, 

F (u) = iF, (u). 

2. Supposons que les fonctions F (u) et G (u) sont les transformées 
de Fourier respectivement des fonctions f (f) et g (t), définies par les 
formules (1.20), (1.21). 

On a formellement 


© 


= | F(u)G(u)e-it" du — : 


| 4 24 


1F (u)e-itu du | g (x) eivt dr — 


—œ 


= | gtodr | F(ue-int-0 qu = 


—©œ 


Li 
PE 


Lee] 


= | emit-T&, (1.26) 


—œo 


i.e. les fonctions 

F(u)G(u) et hk(t)=—— = | |  EOTE—A dx 
forment un couple de transformées de Ro La fonction h (f) est 
appelée produit de convolution (ou de composition) des fonctions f (t) 
et g (£). 

Théorème 3. Supposons que f(t) est la transformée de 
Fourier de la fonction F (t) € L ]—o, +oof et que la fonction g(t) € 
€ LI—oo, + oo (de sorte que sa transformée de Fourier G (t) est une fonc- 
tion bornée). Ceci étant le produit V 9x F (t) G (t) appartient à 
L\—oco, +oof et sa transformée de Fourier est 


h(= À sfr. 
La validité de ce théorème est entraînée par celle de l’inversion 
de l’ordre d'intégration dans (1.26), inversion qui découle de la 
convergence absolue. 


$ 3 PROPRIÊTES FONDAMENTALES DES TRANSFORMÉES DE FOURIER 17 


Théorème 4. Supposons que f (t) et g(t)E L ]—oco, +oof. 
Alors h (t) E L ]—co, +oof ef sa transformée de Fourier est La fonction 
V2n F(t) G (t). 

Pour les propriétés sus-indiquées du produit de la convolution 
moyennant d’autres hypothèses, voir [238], [21]. Par analogie avec 
(1.26), dans le cas des cosinus-transformées, on aura 


| F.(u) G. (u) cos tu du = 
0 


C0 


y ? ( F.(u) cos tu du Ter) COS Tu dT — 
0 


0 


— ( g(T) dx | Fu) [cos|{—Ttlu+ cos (f+ T)u]du — 
0 0 


=7 Setif(t-rl)+ft+od (4127 
0 


et dans le cas des sinus-transformées 
| F, (u) G, (u) sin tu du — 
0 


n'4 Î F,(u) sin tu du ( g (x) sin ut dt — 
0 0 


oo © 
= g(r) dr | F,(u)[cos|{—Tt|u—cos (f + Tt) u] du — 
7e 8 LAS | 


27 
0 


=+ | e@t(lt-sl)—f(t+vidr. (1.28) 
0 


On peut effectuer r fois le produit de convolution. On obtiendra 
alors 


7 | FOR... Fear 
| [ În (Tn) dtn | fn-1 (Tn-1) dtn1 - - 


en? < Le 


| | ht)f(—T—...—tv)dm (1.29) 


2—0394 
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3. Formules de Parseval. Soient f (t) € L ]—oo, +oof, g (t) une 
fonction intégrable sur tout intervalle fini et 
l 


LE tt 
G (t) ei lim jee dt 
pour tous les £{ et, de plus, G (t) est partout finie et appartient 
à L ]—oco, +oco[. On a alors 

_ it dr — 
(] F(t)G(t) dt = = ÿ G(t) dt J f (x) eïtr dx 


O0 


= Lo dx | G (t) eitr di — | f()e(—Tar. (1.30) 


Si en particulier f = g, on obtient 


| IF Pdt= | |f(@) Pat. (1.31) 

Les formules (1.30) et (1.31) contiennent de toute évidence les 

formules des cosinus et sinus-transformées de Fourier. En effet, 
si les fonctions sont paires, on a 


[ Fe (6) G.(t)dt = (4) g(t) dt (1:32) 
à 
et 
[ 1r. (OP dt = { v (OP dt ; (1.33) 
si elles sont a | 
[r. (#) G, (8) dt = ] f () g (6) dt (1.34) 
et 
2 (6) P dt = ( v (dt. (1.35) 
Ces formules sont ar à la de 
Liv (HP dt = + 0 +> (ai + 6?) (1.36} 


de la théorie des séries de Fourier et s che formules de Parseval. 
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&. Transformées de Fourier de fonctions analytiques. Soit f (2) 
une fonction analytique, régulière dans la bande (—a << y << b) 
où a >0, b>=>0. Supposons que dans chaque bande intérieure 
à (—a << y< b) 

Ofet(h-e)x] (x —> — oo), 


f (2) -{ Ofe-à-erx] (x > + ©), 


où & >> 0; À et u étant des nombres positifs arbitraires. La fonction 


(1.37) 


F(u)= ÿ (D) eft dt (1.38) 


vérifie les mêmes conditions par changement de a, b, À, uen à, u, 
a, b respectivement, et 


fais Î F (w) ei" du (1.39) 


Es 
pour tous les z de la bande (—a << y << b). 


$ 4. Transformées multiples de Fourier 


On a par définition 


F (vw, 1)=# {f(x y)] = 271 


OO © 


[ [ eitus+anf (x, y)dz dy. (1.40) 


00 — 00 


La fonction F (w, À) porte le nom de transformée de Fourier de la 
fonction de deux variables f (x, y). Si les fonctions f (x, y) et F (w, À) 
appartiennent à L, on a la formule d’inversion suivante 

f(x v)=F'IF (@, = ( Îe-iurtF (w, À) du da. (1.41) 


—-œo —00 


Si f(x, y) et g (x, y) appartiennent à ZL, alors existe l’intégrale 


oO © 


f (z, ES x ee p=—— | [ f(z, y)8(x—£8, y—n) di dn, 
7 (1.42) 
et, de plus, 
F (f(x, y) X * g(z, y)]= F (w, à)G (o, À). (1.43) 


La généralisation à un plus grand nombre de variables est immédiate. 
On trouvera un exposé complet de la théorie des transformées muiti- 
ples de Fourier dans l’ouvrage de Bochner and Chandrasekaran [21]. 
Voir par ailleurs [421]. 


2% 
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$ 5. Quelques applications de la transformation 
de Fourier 


La transformation de Fourier joue un rôle très important dans 
la résolution d’un grand nombre de problèmes de physique mathé- 
matique, tels les problèmes aux limites posés pour l'équation de 
Laplace, de Helmholtz et de Fourier dans un domaine en forme d’une 
bande infinie, de demi-bande, de cylindre infinie, de demi-cylindre, 
etc. ! 

En particulier, la transformation de Fourier est intéressante dans 
les problèmes qui se ramènent à une intégration d'équations de la 
forme 


+ L(u)=f(& y), 


où Z (u) est un opérateur différentiel linéaire ne contenant pas la 
variable x et f(x, y) une fonction donnée. Examinons quelques 
problèmes de physique mathématique solubles par la transformation 
de Fourier. 

1. Considérons un problème d’hydrodynamique [35], [132] se rame- 
nant à la résolution de l'équation de Laplace 


0°u ou 
Ze tas 0 (y<D0) (1.44) 
sous les conditions initiales et aux limites suivantes: 
= + pour y=0, (1.45) 
u= (x), = 0 pour y=0 et t—0. (1.46) 


Soit 


Ul(w, y, t)=#{u(z, y, t)] = u(z, y, t)e-iux dr. 


ES) 


Siu—+0et _ —+ 0 lorsque [x] — oc, on obtient 


F | = ]= — EU. 


CE 
Et l'équation (1.44) est remplacée par 
ŒU o 


Sa solution qui tend vers 0 lorsque y —+ —o est de la forme 
U =c(w, t)elwiv, 
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où c(w, t) = U|,-,. Compte tenu de cette dernière relation, en 
appliquant la transformation de Fourier à l'équation (1.45), on 
obtient 


c(w, t) = A(w)eiVéloit+ B (w)e-iVelul, 
De (1.46) il vient 
A() = B(o) =7® (o). 


où ® (w) est la transformée de Fourier de la fonction (x). Donc, 
U—=O{(w)cos (V glo]t) ele 
et 


u(z, y, t)— = L D (w) COS Ve | o | t-eloly+iex do. 
2. Soit à déterminer la solution u (x, t) de l'équation de la chaleur 
ôu o?u 


"dt — or? (1-47) 
qui vérifie la condition initiale 
u (x, 0) = f(x) (—o <zxz< +0). (1.48) 


Soit 


U(w,t)=Flu(xz, De Î u(x, theiwdz. (1.49) 


— 00 


Moyennant les mêmes conditions que dans 1, il vient 


F [  |= — QU (w, t). 


"Êx? 
Donc (1.47) et (1.48) se ramènent à l'équation 


OÙ (@, t) 2 
DE — © U (w, t), 


d'où 
U (o, t)= A(w)e-%*t. 
En faisant t—0, on obtient 


OO 


77 } 


À (©) — f(z)ev dr F (w). 


Donc, 
U (w,t)=F(w)e-®* 


te 
tw 
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et 
u (x, D=—— [ F (w) e-%*t-ixe Go. 
Par ailleurs _ 
u(z, t) = —— ï et 0 do ( Î (È) eivé dE — 


nd tél 


— 00 00 
(x—E6)3 


4. 
TE j f@e" # &. 
3. Soit à déterminer une solution u (x, t) de l'équation (1.47), 
telle que 
u(z, 0) =0 (20), u(0, 2 =f() (> 0). 
Posons 


U, (w, t) — A u (x, t)sin wzx dx. 
0 
Alors (1.47) se ramène à la forme 


700 — ne L wf(t)—@U, (w, t). 
D'où 


ue t 
U,(&, = A(u)e-c®t J/ Zoe-ett [e-usf (x) dr. 
0 
Comme U,(w, t) = 0 Éd t = 0, il vient À (wo) = 0. Donc, 


u(xz, t) = 3: Ee-£&*t sin Ex dE | et"tf (+) dr = 
0 


=+ {io | Eett(r-9 sin Er dE = 

‘ | | { 
= f(r) G—+) 2e AU-T) dr — [ V(z,t—7T) f(x) dr, 
2 ; ] 
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« 


ou : 
A TT de 
+? 
(x, t)= Ce, 
X(Z, 1) = —— 
7. 


Pour justifier les solutions formelles précédentes, il suffit de suppo- 
ser, par exemple, que toutes les fonctions envisagées appartiennent 
à L ]—oo, +ol. 

&. Dans certains cas la transformation de Fourier permet de 
résoudre des équations intégrales de la forme 


pa f(0)+ | k(z— y) pu) dy, (1:50) 


— œ© 


où f (x) et k (x) sont des fonctions données, q (x) la fonction cherchée. 
L'application de la transformation de Fourier donne 


Du) = ÿ {f(e)+ jee — y) @ (y) dy } ei** dr = 
=FU+- J sta | ken dr 
=F(U)+—= e) p (y) dy J #6 eiu+nu an — 

| = F(u)+V2x O(u) K (u). 

D'où | 

… F (u 
DE TETE TA 
et 


F(u) 
p(G)— FE) 1—V 21 K (u) 


D'un autre côté 


e”'*4 qu. 


p(z)— f(x) = —— = j (ere FU} et du= 


_ K(u) ir 
j F0) are" du. 
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En introduisant la notation 
_ K (u) 
Fe (4) 1—V 2x K (u) ? 
il vient 


pa)=f(a)+ | r(2—E) j(E) dE. 

Si dans (1.50) les fonctions œ (x), f (x) et k (x) sont nulles pour 
les valeurs négatives de l'argument, on arrive alors à l’équation 
pa)=f(+Tk(G-HDpEOSE (>0) 

0 
Dans ce cas la solution est de la forme 
p(z)=f(2)+ | Or (5 d8, 
0 
où r (x) = 0 lorsque z << 0 est une fonction dont la transformée 
de Fourier est 
K (u) 

1—V2xKk(u) 
On remarquera que ceci entraîne 


R(u) = 


r (2) =k(2)+ | #O r (xD d8. 
0 


Une autre équation intégrale soluble par la même méthode est 


© 


fa= À k(z—y) p (dy. (1.51) 
Par analogie à ce qui précède, on aura formellement 
Fe | au [ke-De@&- 
== | pd | k(—5) eau = 
n Le] Lee] 


= | POSE | Ame dn= V2r O(u)K (u), 


-0 —œo 
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d’où 


_ 1 C F(u) -ixu 
p(x)= [ Re" du. 

On remarquera que pour que cette méthode soit rigoureuse, il 
faut que toutes les fonctions figurant dans l’équation vérifient des 
conditions spéciales. On a, par exemple, le théorème : 

Supposons que f(x) E L? ]—o, +of, k(x) E L ]—o, +ol. 
Une condition nécessaire et suffisante pour que l'équation (1.51) possède 
une solution œ (x) appartenant à L°? ]—o, +ool esf que re € 
EL? 1—o, +ol. 

Voyons encore une équation intégrale de la forme 


pa=f(a)+ | k(+ 0) page 


Comme précédemment on a 


DU= FES | eds | E(e+u) pu) dv = 


= F(u)+ 7 puy) dy | k(z+y)et ar = 


=FU)+ 


] 
j @ (y) dy j k (n) eitn-uu an — 


FEV O(—u) K (u). 
En inversant le signe de u, on obtient 
O(—u)=F(—u)+V2rO@(u) K(—u). 
Les dernières égalités entraînent 


__F(u)+V2xF(-—u) K(u) 
DU=——- KE kCS 
Donc, 


_ Q FG)+HV2RF(—0) K(U) ice 
Poe | or du. 


Dans nombre de problèmes de physique mathématique on a à con- 
sidérer des équations intégrales de la forme 


PU =f(z)+ | k(z— y) 6 dv 
0 
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à noyau symétrique, dépendant du module de Ia différence des deux 
arguments. L'étude de cette équation se fait également à l’aide des 
intégrales de Fourier (cf. [74]). 

5. Calcul des intégrales. Pour calculer des intégrales on se sert 
des formules de Parseval (1.30) à (1.35). Pour le calcul de certaines 
intégrales contenant des fonctions trigonométriques et exponentielles, 
on peut se servir des formules d’inversion correspondantes. Suppo- 
sons que 


S = | e "2 J; (x) 45. 
0 
En vertu de (1.24), on a 
-ax 2 
f@)=e, F(Uu)- y À 7 PT d 
0 pour 0<u<i, 
1 
g(z)=2J, (x), G,(u) — 2 + 
1 
r(z—) 
(-h<res<i). 
En appliquant la formule de Parseval (1.34), on obtient 


(u2—1) "72 pour u>1 


VE 
fe ea T, (x) dr = Vif CT NP REENE"E Cu 
d r (5) (x—1) ” 
: PCs. Pi at 1 | 
arf) nina Vars 


Le proongement étant analytique, ce résultat est valable pour 
RevV>—+ L (2381. 
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CHAPITRE II 


TRANSFORMATION DE LAPLACE 


$ 1. L'intégrale de Laplace 
et ses propriétés fondamentales 


Soit f (£) une fonction de variable réelle # € [0, col, intégrable- 
Lebesgue sur tout intervalle J0, AÎ. Soit un nombre complexe p — 
= 6 + it. L'expression 

fe (p)= À e-Ptf (0) dt = &1f (E) (2.1) 
0 
s’appelle l’intégrale de Laplace, et la fonction f* (p) la transformée 
de Laplace de la fonction f (t). Voici les propriétés fondamentales 
de l'intégrale de Laplace. 

1°. Si l’intégrale (2.1) est convergente en un point p,, elle le 
sera en tous les points p tels que Re (p — p,) > 0. 

Trois cas sont à distinguer : 

1) L'intégrale est partout divergente. 

2) L'intégrale est partout convergente. 

3) I] existe un nombre o. tel que l’intégrale est convergente pour 
Re p >> ©. et divergente pour Re p << ©.. 

Dans le plan complexe, la droite Re p = o. est appelée droite de 
convergence, et le nombre ©. l’abscisse de convergence de l'intégrale (2.1). 

20. Si l'intégrale (2.1) est absolument convergente en un point 
Po = Oo + ito, elle le sera absolument et uniformément dans le 
demi-plan Re p > ©. 

On définirait de façon analogue la droite de convergence absolue 
Re p = 6, et l’abscisse de convergence absolue o,. De toute éviden- 
ce, 6, > cet il est aisé d’exhiber des exemples où 6, > ©.. 

3°. Si l'intégrale (2.1) est convergente en un point pos = 60 + 
+ ito, elle le sera également dans un domaine A défini par l’inégalité 


[P— pol Sk(o—09) eX6-00, 5200 
où Q > 0 et k > 1 sont des constantes. 


49. Si ©. << ©, l’intégrale (2.1) est une fonction analytique de 
la variable p en tous les points du demi-plan Re p > 0. et 


dM*(p) y pen et 
_ “1 1" f(t)ePt dt. 
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5°. Soient f; (p) et f: (p) les transformées de Laplace des fonc 
tions f. (£) et f, (t). Si les deux intégrales de Laplace sont convergentes 
en un point Po et 


fi @o + nl) = f2 (Po + nb), 


où l'est une constante >>0 et n — 0, 1, 2, . . ., alors on aura partout 
fa (€) = fe (*). 

Cette propriété entraîne que la transformée f* (p) de Laplace 
définit la fonction f ({) de façon unique, à un ensemble de mesure 
nulle près. 

6°. Si l'intégrale (2.1) converge en un point po = Oo + it 
O9 > 0, alors 


t 
lim e-cot | f (u) du =0, 
{—00 0 


{ 
1.e. \ Î (u) du —o (est) lorsque t —+ oo. 


0 
7°. Si: a) f (t) est bornée inférieurement, i.e. il existe un nombre 
positif C tel que f (t) > —C pour tous les t > 0, b) existe l’une des 
limites 
’ e co 
lim — | f(t) dt ou limo | f(t)e-st dt = lim of” (o), 
e0 € 0 C—o 0 O—+>00 


l’autre existera également et 
’ e 
lim Fu | f(t) dt= limof” (o). 
0 go 


8°. Si: a) f (f) est bornée inférieurement et b) existe l’une des 
limites 
._ AS C L 
lim + | f(t)dt ou lim o | f(t)e-°t di, 
î U Fe Ç 


£— 00 


l’autre existera également et 


€ 
lim + | f (t) dt = lim of* (o). 
e-o € À o—0 

Les deux dernières propriétés de l'intégrale de Laplace découlent 
de la théorie générale des théorèmes taubériens [256], [257], (2581. 
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Une condition nécessaire et suffisante de convergence de l’inté- 
grale (2.1) est que pour un 0, > 0 et £ —> co 


fi (#) = | 7 (u) du =o (est), 
(ul 
i.e. 


Le 
lim e-5ot Î f (u) du — 0. 
{00 


0 


On a déjà dit que la transformée de Laplace définissait f (f) de façon 
unique (à un ensemble de mesure nulle près). Voyons maintenant 
comment déterminer f (t) si l’on connaît f* (p). 

Théorème 1 (théorème d’inversion). Si l'in 
tégrale (2.1) admet ©, << © pour abscisse de convergence, alors il existe 


 . O0 pour t<<0, 
mi (y 
sn 2ri jo (p) + gr [ f (u) du pour t>0, (2-2) 
0 
où Y>O,, Y> (0. 
Donc, pour presque tous les t 
a ; Y+io0 
[= [ f* (p (2.3) 
Y—ico 


où l'intégrale est prise dans sa valeur principale. 
Remarque. La propriété 6° entraîne 
| Aer, 


4 
où f1 (Ë) = ff), o>o, os >0etp=0+ix. Il existe une 


0 
constante Q telle que |f, (f) | << Qest (o > 0.) pour tous les tf. 
Par conséquent, 


[ER f° (p) PRIT _ (2.4) 


Donc, si 


t 


= [per ar, > et fi(#)= | f(u) du, 
0 
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la transformée de Laplace de la fonction j, (t) sera —— P @ 


et, de plus, 
l'intégrale de Laplace convergera absolument pour © > o.. Donc, si 
{ En 


fn (9) = | dt | Ba (Gi = 
0 0 


0 


t 
[OU &, 
0 


alors 
fs) _( 7 
DA = [re dt 

et 

Y+io0 

1 + 
f (t) = _ Ori Î Î Y>O. (2.5) 
Y— io 


L’inégalité (2.4) entraîne que l’intégrale de (2.5) convergera 
absolument et uniformément pour rz — 3 sur tout intervalle [a, bl]. 
Il est évident que cette intégrale convergera d’autant plus rapide- 
ment que le nombre n sera grand. 

En général, on calcule les intégrales de (2.3) et (2.5) en déformant 
de façon convenable le chemin d'intégration; souvent il est possible 
de se servir des lemmes suivants. 

Lemme 1 (Jordan). Soit C, l'arc de cercle 


3: ; 
|21=Rns <-<arg 2; lim R, = co. 
Eu n—00 


Si sur Ch, une fonction © (z) d’une variable complexe z tend uniformé- 
ment en arg z vers O lorsque n —> ©, alors 


lim [ O(z)edz=0 pour t>0. 


Ca 


Lemme 2. Supposons que lim R, = œ et que Cr et Cr 


ñn+00 


sont respectivement les arcs de cercle 
|zl= R,, 0O<KRez<Yy, Imz>0, 
|zl= Rr, 0O<Rez<Y Imz<0. 


Si sur Les arcs Ci, Cn (n = 1, 2, 3, ...) une fonction © (z) est 
uniformément bornée et tend vers O lorsque n —> co, alors 


lim | @(z)et 40 et lim UE 
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Lemme 3. Soit ® (z) une fonction analytique régulière dans 
le demi-plan Re z > y. Si sur les arcs 


[zl= Rs, — + —<arg:<—-(limR, — co), 


la fonction ® (z) tend vers O lorsque n —- oo uniformément en arg z, 
alors pour t << 0 


Y+iw Y+ioc 
lim ( O(z)e dz— | ® (z) ee‘! dz—0. 
nue Yv-io y ic 


Voyons comment on calcule une intégrale du type (2.3). Soit 


£1f(t, ]= = e-"Vr*), 
1.0. 
a+ico; 
ft = ( in 7 (2.6) 


En vertu du théorème de Cauchy, l’intégration le long de la 
droite L(a—ib, a + ib) équivaut 
à une intégration le long du con- 
tour composé des arcs Cret Cr du 
cercle | p| = R, des bords B, et B, 
de la coupure et de l’arc de cercle c,: 
[Jpl=r(—x < argp<n) (l'inté- 
gration a lieu dans le sens indiqué 
par les flèches sur la figure 1). Comme 


. et VP 
tT>0, la fonction ————— —+ 0 lors- 
V’P 
que À —+ oo sur les arcs Cr et CR. 
Donc, en vertu du lemme de ‘Jordan 


et Vr+rt Fig. 1 


l'intégrale de l’expression 


P 
prise le long de CR et Cr tend vers 0 lorsque t > 0 et R —+ co. 
Par conséquent, 


f(#, D= lim — | +(+ [} be-tVrtpt 7 
B2 V 


B1 oc, 


*) S'agissant des fonctions multivoques 2% — ealnz, În z, arctgz, etc., 
sauf mention du contraire, on considère toujours leur détermination principale, 


i.e. In 1 = O0 (arg 1 = 0), arctg 1 — . etc. 
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Le long du bord B,, on a p =#zre-it, Wp = —i Vz; le long du 
bord B,, on a p = re, Vp = i V x. Donc, 
r Fe [R 
— itVx-xt = = — it Vz-xt 
Î [ e =" | [ e 
B:1 R 2 ? 


Sur l'arc c,, lorsque r —+ 0, on a 


Cr 


dz 
iVz° 


2sr 
< M 7 


Donc pour t>0 


f(£, t) =+ | e* cost V/ x TE = + [ e-tu® cos tu du = 
0 0 


(2.7) 


1 , 
où z=u. Pour t<<0, les intégrales de -7=€ *VP+Pt sur CR 
et Cr (R — co) tendent vers 0, i.e. f (&, t) = 0 pour £ < 0. Donc 


A ne ete Vo 
AE = | (2.8) 

Théorème 2. Si l'intégrale (2.1) est absolument convergente, 
alors lim f* (o + it) = O0 et la convergence est uniforme pour 


T— +00 
tous les o (0 > O1 >> Go). 

Théorème 3. Si l'intégrale (2.1) est absolument convergente, 
H (z) une fonction analytique au voisinage de chaque point z = f* (p) 
et H (0) = 0, alors la fonction O (p) = H [f* (p)] est représentable 
par une intégrale de Laplace absolument convergente sur le demi- 
plan Re p > o4. 

Il existe d’importants critères qui nous permettent d'établir si 
une fonction donnée (analytique dans le demi-plan Re p > y) est 
une transformée de Laplace ou non. Dans de nombreux cas, le théorè- 
me 3 nous fournit une réponse à cette question. Par exemple, l’inté- 
grale | ePi dt — - est absolument convergente lorsque Re p > (0. 

Ô 


D’après le théorème 3 il vient que | 


Vr+1 
ment convergente de Laplace. Ceci entraîne la représentabilité de la 


— 1 est intégrale absolu- 
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i 
fonction — 1) eVr+i et ainsi de suite. On a en parti- 
Vr+i 

culier le 

Théorème 4. Une fonction analytique régulière au voisinage 
d’un point situé à l'infini et nulle en ce point est représentable par une 
intégrale de Laplace absolument convergente. 

Enonçons d’autres théorèmes de même nature. 

Théorème 5. Supposons que dans le demi-plan Re p > y 
une fonction analytique f* (p) remplisse les conditions: 

4°. lim F(oHit) _ 0, 


Te +00 O—+ IT 


Co >; 


la convergence étant uniforme dans le demi-plan. 
2°. Pour tous les tE]—oo,—+oo [, il existe 
; o+1iw pe 
lim —— ( LL pt q —©O (t). 
sé ne : i 


3°. La fonction ® (f) est absolument continue et existe l'intégrale 


F (p)= | D" (t)e-Pt dt. 
0 

Alors f* (p) = F (p) et, par conséquent, f* (p) est une transformée 
de Laplace. 

Théorème 6. Si une fonction f* (p) est analytique dans le 
demi-plan Re p >> y et bornée dans chaque demi-plan Re p > © > 7 
et si pour 6 >> y existe l'intégrale 

| [ff (o+it)l dt<oo, 1<r<2, 
alors f* (p) est représentable par l'intégrale de Laplace dans le demi- 
plan Rep > Y. 

Théorème 7. Siune fonction f* (p) analytique dans le demi- 

plan Re p > 7 vérifie la condition 


sup [ [f" (o+ it) |” dr < 00, 
O>Y su 
où 1Lr << 2, alors f* (p) est représentable par l'intégrale de Laplace 
dans ce demi-plan. 
Théorème 8. La condition 
sup Î | f* (o + it) [? dt << 00 
O>Y Se 


3—0394 
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est nécessaire et suffisante pour que la fonction f* (p), analytique dans 
le demi-plan Re p >> y, soit la transformée de Laplace d'une fonction 
f (t) telle que 

| [f(t)[2e-21t dt < 00. 

0 


Théorème 9. Supposons que 
1°. f* (p) est une fonction régulière dans toute région finie bis plan 


de la variable complexe p, Que aux pôles P1, Pes P 
(IP1<|Pel <Psl<... LIml<: y de la fonction Ÿ* (p). 
où Re p, < 6. pour tous Les : R 
20. Existe la 
1 ep) 1 VTT pe (p) 
Dee t 
im | Ebetap= | LBetap, 
Y—-1v0 Y—ico 
V>Ga Y>0. 


30. Existe une suite de contours simples C, d'extrémités les points 
y + iB, et y — iB, de la droite Re p = y. (Ces contours sont contenus 


Fig. 2 


dans le plan Re p < y et ne passent pas par les pôles p,.) Chaque 
contour C, renferme l'origine des coordonnées et les n premiers pôles 


P1 De: Ps . 3, Pn (fig. 2). 
4°. Pour tous lest >0 


Lu ar | Î" (P) = dp=0. 
Alors l'intégrale est égale à la somme d'une série convergente, soit 


Y+io0 
7m | ÉPetan= Z rt» 


Y—ic 
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où r, (t) est le résidu de la fonction r CP ) ePt qu point p = Pr (n — 
— 4, 2, ...)et r,(t) le résidu au er 0. 


Remarque. Si la fonction Li (p) vérifie les conditions 


des lemmes 1 et 2, il semble naturel de choisir pour C, des arcs de 
cercles centrés à l’origine des coordonnées. 


Fig. 3 


Si existent un nombre @ >> 0 et une suite de nombres positifs f, 
et Ôh > 0 tels que 


1) 
limf,—=, limô,—=0. 
2) 
f* (0 + iPn) f°(—Bn+tr) 
DT CR < Ôn; BF, Fix <Q, 


—P,So<y et [T|<h:, 


alors pour contours C, on peut prendre le contour en forme de cro- 
chets représenté sur la figure 3. 


S 2. Théorèmes de convolution 


On appelle produit de convolution ou de composition des fonctions 
a (t) et b (t) de la variable réelle £ la fonction c (t) définie par 


t 
c(t) = | a (t—7+) b (x) dr. 
0 


Ce produit est noté symboliquement 
c(t) =a(t) x b (+). 
3® 
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Voici quelques propriétés du produit de convolution. 
a) Commutativité 

a(t) X b(t) = b (t) x a (1). 
b) Associativité 

(a kb)kc—=axk(bxc). 


*« 


c) Distributivité par rapport à la somme 
la) + b(G)xkc (Er) = a (8) xc (E) + b (6) xc (+). 


d) Théorème 10 (Titchmarsh). Si le produit de 
convolution de fonctions a (t) et b (t) continues sur l'intervalle [0, œl 
est identiquement nul, l'une au moins de ces fonctions l’est également. 

Ce théorème a été démontré par Titchmarsh en 1924 [238], [236]. 
PA démonstrations ont été proposées par la suite [61], [204], 

Théorème 11 *) (théorème de convolution). Si les intégrales 


ft) | heat et fi(p)= | fi(b eat 
| 


0 


sont absolument convergentes pour Rep>o,, alors f*(p) = 
— fi (p) f5 (p) est la transformée de Laplace de la fonction 
t 
1@)= À AGE) fe (0 dr 
0 
et l'intégrale 


f'(p)= | FE) ert ai 
0 


est absolument convergente pour Re p > ©. 
Enonçons ce théorème sous la forme suivante. 


Théorème 11’. Si I) ED y LME Sont les transfor- 


P 
mées de Laplace respectivement des fonctions f (t), g (t) et h (t), alors 
presque partout 


t 


d 
k (t) =+ | f(© z(t—7) dr. 
0 
Faisons une remarque relativement à ce théorème. Supposons 


ue 
. Bi 


fe fer#h@d, > (2.9 
œi 
*) Ce théorème est parfois appelé théorème de multiplication ou de Borel. 
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et 
B2 
fi(p)—=\ e'fo(t)dt, Br >. (2.10) 
Il vient alors ” 
Bi+B2 
FO) (p)= | ef (0) de, (2.11) 
ai+a2 


min(B1: t-@z2) 
1O= | hoktG-Ddr. (2.12) 

max(œ1; {—f2) 
La démonstration de (2.11) et (2.12) utilise l’intégrabilité absolue 

de (2.9) et (2.10). 

Théorème 12. Soient données deux fonctions f (t) et g (t) 

d'indice de croissance s, et s:, i.e. 

[F()I<AMent, |g(t)|< Met. 


La transformée de Laplace du produit f (t) g (t) est la fonction 


a+ico 
mr | le (p—2) ds, 


a—ioco 


où 


a>s, et Rep>s+a, 


O0 


jp e"tf(dt, g*(p)= [ e-Pt g(t)dt. 
0 0 


A. Efros a démontré en 1935 un théorème de multiplication 
généralisé qui est d’une grande importance. Enonçons-le. 
Théorème 13*). Soient f* (p) — | e-vtf (t) dt et g* (p) et 
0 


q (p) des fonctions analytiques telles que 


g"(pe-tin= À e-Plg(s, © dt. 
0 


Alors 


Pate" (= [era | et, Dr. 
0 0 


*) Les hypothèses exactes de ce théorème sont données au chapitre V, $ 12. 
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Si en particulier on pose g(p)= p, | e-Pltg(t, t)dt—e-Pg* (p), 


: ) 
1.0. g(t, t)=g(t—7+), on obtient (pour t>t, g(t—Tt) =0) 


f" (p)g* (p) — | e-Pt dt | f(x) g(t—T) dr — 
Ô 0 


= [e-n a f(r)g(t— 7) dr. 
0 0 


$ 3. Quelques propriétés de la transformée de Laplace 


Voici quelques propositions simples constituant l'appareil de la 
méthode opérationnelle. Dans toute la suite on adoptera les notations 


fe (p)= | e-ntf(t) dt= & 1j (+1, (2.13) 


0 
F(p)=p | e-ntf() dt = C1f(b). (2.14) 
0 


19. Propriétés de linéarité. Soit 


f(t)= 2 cufn (é), 


k=1 


où c, sont des constantes (complexes) arbitraires. Alors 


LIf()1= £ [2 Chfn (t)] 2 Ch£ [fx (£)] 2 Chfk (P) = f° (b). 


(2.15) 
D'après (2.15) on a formellement 
d .ÀA+d})—ft{(t, À 
£ [+ 1 €, 1) |=£[- EEE | 
f” (p; À + d\)—#f* (p. }.) — d 
= ——%@—— — <= f (P; À), (2.16) 


A2 


f(£, à) dà | = | LI (, À] d\= 


Ai 


f°(P, A)dA. (2.17) 


z| 


Des propriétés analogues ont lieu pour la transformée de Laplace- 
Carson (2.14). 


Le > 
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2. Propriété de similitude. Pour toute constante «, on a 


s[(2)]= fr (4)e-manal menace ar (ap, (248 
0 0 

C Ë (£)}=(: (=) e-Ptdt=— ap À st e-er de = 7 (ap). (2.19) 
0 0 


3°. Transformée de Laplace de dérivées. Une intégration par par- 
ties nous donne aussitôt 


Z (9 ()1= p'f° (p)— pr tf (0) — p-?f" (0) — p-*f" (0) — 


2. — pfin-2 (0)— ft (0), (2.20) 
Cf (1 = pr (p)— p'f(0)— p°-1f (0)— pr-2f" (0) —… 
2. — pfn-2 (0) — pjt-1 (0), (2.21) 


où x est un entier positif. 

La propriété duale de 3° est 

40. Différentiation de la transformée de Laplace. Pour n entier 
positif, on a 


OT) | mf(the-rtdt=(—1} L{f(#)], (2.22) 
0 
- 
Et —(—1ÿC Er ()—n em-1f (t) dt | | (2.23) 


5°. Transformée de Laplace d'intégrales. Pour n entier positif on a 


£ É dr dt + f (tn) dtn- | = re | (2.24) 
0 0 0 


6°. Intégration de la transformée de Laplace. Si l'intégrale 


| f* (g) da est convergente, elle est la transformée de Laplace de la 
d 
fonction 


10 » 1.0. 
(r (4) d=£ [2]. (2.25) 
P 


Pour tout nr entier positif, on a de toute évidence 


(da | an … À fe (qu) dan = £ [ L2 x (2.26) 
P q an 
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Citons encore quelques formules de la même nature. Une inté- 
gration de l’expression 


f" (&p) = TE e-Pt dt 
0 


par rapport à & entre O0 et 1 donne 
1 00 1 1 
Jrenanenals(1)Es(|(2)E] 


En posant ap=q, t—= ax, il vient 


£[ (ie dr | =+ { f* (g) da. (2.27) 
De façon analogue, | | 
2 La |= D] f* (g) da. (2.28) 
D'où 
£[ Ï = [r (a) dg, (2.29) 
IQ H= f* (a) da. (2.30) 


7°. Pour tout + positif et puisque f (£ — rt) = 0 lorsque t << 7, 
il est aisé d’obtenir 


Lij(—n1= | f(t—rTe-rtdt= 


A2 9 


— | f (u) e-Ptu+0 du = e-r | f(u) e-rv du, 
0 0 


i. e. 
LIf(t—7T)]=e-rÿ* (p). (2.31) 
8°. Pour tout q complexe, on a 


f(p—g)= | f@e-v-atat- | 1f()esle-rtdt= £ [f()e], (2.82) 
0 0 
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f(p—g)=(p—0) | f(t)e-@-di= p | [f (t)ett]e-Pt dt — 
0 0 


t 
—<p fifbetie-nta=c{[iet- | med]. (233 
0 0 


Formulons encore deux théorèmes importants qui nous seront 
très utiles pour résoudre un grand nombre de problèmes pratiques. 

Théorème 14 (premier théorème de déve- 
loppement). Si une fonction f* (p) est régulière au point infini 
[132] et admet en son voisinage le développement de Laurent 


f" (p)= À) > + 
k=1 ph” 
alors 
Ch Ch ,k- 
2=-2(>Dr Le (2.34) 
h—1 k—1 
Ceci étant 
f(t) — De T— ant 


est une fonction entière. 

Théorème 15 (deuxième théorème de dé- 
veloppement). Soit une fonction f* (p) vérifiant les conditions 
suivantes : 

40. f* (p) est méromorphe et régulière dans un demi-plan Re p > so. 

20. IT existe un système de cercles emboîïtés 


Ch: Ipl= Re R<Re<... Ri +00, 


sur lequel f* (p) tend uniformément vers O par rapport à arg p. 
30. Quel que soit a >> 50, est absolument convergente l'intégrale 


a+ioco 


Î f'(p)dp, 


et 


f*(p}=Æ£(Dres,, f° (p)e°!), (2.35) 
Ph 


où la somme des résidus est calculée en tous les points singuliers p4 
de la fonction f* (p) dans l'ordre de non-décroissance de leurs modules. 
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Corollaire. La fonction homographique 


1, MP) _amPM+am3p" + ...+ap+ao 
f(p)= NP) bnPn—bn-ypt +... bip+ bo ” Se (2:08) 


est la transformée de Laplace de la fonction 
- 1 "HR S 
= >» —— lin ———— ({f — kePiy, (2.37 
( ) 2 (rx —1)! Dh db R=i {f (P) (P Pr) } ( ) 


où pa sont les pôles de f* (p), rx leurs multiplicités, la somme des 
résidus est calculée en tous les pôles. 

En particulier, si tous les pôles de f* (p) sont simples, en calcu- 
lant les résidus en ces pôles, on obtient 


M 


Si les polynômes M(p) et N oi possèdent des coefficients réels, 

alors 
nn M{(?) _ M (px)_ pat M (px) op, 

rh=pé=£(>eer+2RD Er en), (2.38) 
où la première somme s’étend à toutes les racines réelles de N (p), 
et la seconde à toutes les racines complexes à partie imaginaire 
positive. 

On remarquera que dans la formule (2.37) le terme qui corres- 
pond à la racine complexe p;, = 6, + it, a pour expression 


M (Ph) e°ht 
N° (Pr) 


(cos T4t + i sin tAt). 
$ 4. Transformées de Laplace de quelques 
fonctions élémentaires 
Par définition de la fonction eulérienne L' on a 
T (k+ 1) = | et dt, Rek>—1. 
0 
En posant t — pt, on obtient 


UE) | e—Ptr* dr. 
0 


ph#i 
Donc 


L () = mL Rek>—1. (2.39) 
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En particulier, pour les 4 entiers non négatifs 
£ (t)= A (n= 0, 1, 2, ...). (2.40) 
Une dérivation par rapport au paramètre k de (2.39) donne 
£ (* In = Lb(&+1)—Inp], Rek>—1, 


où Ÿ est la dérivée logarithmique de la fonction T. En faisant 4 = 0, 
on obtient 


£ (nt) = —+(C+np), 


où 
C=—ÿ(1)= lim (1++ FR .++—inn)= 0,577215665 


est la constante d'Euler. 
En se servant des développements 


a n #2? . æ tin+i 
cos { — pi (— 1) Gr: sint= >, (—1) Cri? 
n=0 T==( 
et compte tenu de (2.40), on obtiendrait 
L (cos t) — FE s “L (sin = Tr. (2.41) 


h 
Cherchons l’image de la fonction t © J, (2 Vt), où 4 est un nombre 
complexe, Re 4 > —1, et ee (t) la fonction de Bessel d'ordre k 


HnE=(z) à rer). 


En posant t—2V/z,on a 


es] 


— {y en 
n(2V2)=:7 >, 1 


nl  T(n+k+t): 


n=0 
—iyn entk 
£ (7 (2VT) = 2(3 CT TRIER 
1 
=; CN De 7, Rek>—1. (2.42) 


prth+i phti 


Introduisons les ne. de Laguerre 


Lntt)= (pet) (n—=0,1,2, ...). 


n! 
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Compte tenu de (2.33) et (2.40), il vient 
- nt — n 78 
L {g (t)}= £ {{ € } — (p+ 1)7+1 (P). 
Comme g (0) = g’ (0) = g"(0) = ...—gt"-12(0) —0, moyennant (2.20) 


on trouve 


an ,n -t … n\pn 
£ { dtn (4 € )} = (p+ 1)7+1 x 
Et enfin, en se servant encore de (2.33) 


L {Ln @}=< (1-2) (n=0, 1, 2,...). (2.43) 


$ 9. Calcul d’intégrales 


Citons quelques exemples concrets. Soit 


f(t)= 1. usine, 10. (2.44) 
0 
En appliquant (2.41) et (2.33), il vient £{usintu}= Er, 
u° du A _ 1 A -1 
2 {@)= arte rite}. 
Donc, f (t)=< et. Considérons l’intégrale 
I Ï RO, kt1>z 0. (2.45) 


En posant u—2Vt, on trouve 


RTE 


CEREVD 
— JR 2R= RENE dt. 
0 


Introduisons l'intégrale Z (À) — Lui nan a 


0 t 2 


En vertu de 


(2.42) et (2.18), il vient 
1 


k … tb \RHi 
LEP RCVRy=(G) er. 
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Donc 
Lu pass, r() ,Jr(G) 1 
LUON= x |! e a LE ; 
d’où 


nn + "(9 


PR PR HT 


Calculons enfin l'intégrale 


C Jo (t) — cost 
Ÿ 200 à 


où J, est la fonction de Bessel d'ordre zéro. On a 
L {Jo (t) — cos t} — Von Pi 
Donc, en vertu de (2.30) 


jeta | (Er) dr 
— [in er 2 [=in 2. 


$ 6. Application de la transformation 
de Laplace à la résolution des équations différentielles 
et intégrales 


1. Soit donnée l’équation différentielle 
anut9 (t)+a,nmum-9(t)+ ... +<au’(t)+aou(t)=f(t), (2.46) 
où u (t) est la fonction inconnue de la variable indépendante f, f (t) 
une fonction « perturbatrice » donnée et a; (i = 0, 1, 2, ..., n) 
des coefficients constants [171]. En multipliant les deux membres 
de (2.46) par e-?* et en intégrant par rapport à t de zéro à l'infini, 
il vient : 
a* (p)u* (p)—b"° (p)=f" (p); (2.47) 
où 
a* (P)=anP" +an1P + ... + @P+0@0, (2.48) 
b*(Pp)=bnap + ... + b1P + bo, (2.49) 
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et 
On = nu (0) , 
On-2 = nu” (0) + an-au (0), 
bn-3 = nu” (0) + a;-1u" (0) + an-eu (0), 


b, = anutr-2) (0) + an utr-3 (0) + ... + au (0), 
bo=anum- 1) (0) +anum-2 (0) + ... 
. + au" (0) + au (0). 


En résolvant l'équation (2.47) par rapport à u* (p), on obtient 
la formule suivante : 


1 1° ()--b* (p) 
PERS (2.51) 
Posons 
(D) = 0) 50. (2.52) 
Donc 
u°® (p) =f" (p)r° (p)+s* (p). (2.53) 


Les quantités r* (p) et s* (p) sont des fractions rationnelles qui, 
par des méthodes standard, peuvent être décomposées en fractions 
élémentaires. Le théorème de convolution donne 


u (t) = | f(n)r(t— +) di+s(t) (2.54) 


0 


Nous avons obtenu la solution générale de l’équation (2.46). Cette 
solution contient x constantes arbitraires qui sont les valeurs ini- 
tiales de la fonction inconnue v (t) et de ses nr — 1 dérivées. L’allure 
de cette solution dépendra de la nature des racines caractéristiques 
de l’équation 

a* (p) = 0. (2.55) 


1°. Toutes les racines de l’équation (2.55) sont réelles et dis- 
tinctes : 
a* (p) = An (P — P1) (P — P2) - .. (p — Ph). (2.56) 


On a alors 


s = A 
rer 
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où les coefficients constants r, et s, se déterminent par les formules 


1 ge —-0? (Pn) 
a" (pr) À a’ (px)? 


TR — 


b* (px) =u(0) 2 api-t+u (0) À api? + — 
..+u-2 (0) > api-"+i + um-t (0) a. 
I=n-1 
Donc, 


r (t) => rePn!, s(t) > srePh'. (2.57) 


En portant (2.57) dans (2.54), il vient 


t 


u(t=S on || fe a de+b" (p)]. (2.58) 


20. L’équation (2.55) possède des racines nulles : 


a* (p) = anp" (2.59) 
et 
: 1 
r er 
bn 4 , bn 1 
Fe Da Don rie _ 
Alors 
ini 
RUE TE: 


in—2 pn-1 


D +. +2 + GT 


Dans ce cas PR (2.54) s'écrit 


s(t) = Pn= 


t 
1 — T)n- 1 bn-1- k 
= E [ re Eat S de (2.60) 


30. Les racines de (2.55) sont réelles et confondues, i.e. 


a* (p) = an (P — Pi)”; (2.61) 
donc 
: _ 1 
CO) = pr 


°(p) = 2° @) En + ni _4 
| = PT (= = pe pet + Ho À 
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où c4 sont des fonctions homogènes linéaires des données initiales 
que l’on définit par des méthodes standard en décomposant des 
fractions rationnelles en fractions simples. On trouve 


pe t t ù Ch . 
=" à s(t)= er1 2-T" 1 
La formule (2.54) se détaille comme suit 
t 
t— 7)" 1, ,P4( 
u (= | f(x) CT + 2 aër mr. (62 
0 


Dans le cas général, le polynôme a*(p) se présente sous la forme 


a* (p) = an (p — pr) - + (p — pa) (P — Pan)" P° + bip + ci) -.. 


. (° + 0x p + cn) (P° + usa P + Che)" +. . (2.63) 
Donc 
$ ! k un 
. VEEY 
de = à P—Pa Pa rer (P— + y PT bvyP+ Cy + 


CsPp-- Ds 
+ nn FE NE 
> (P° + bns1P + Ch41) LÉ 


LA 
PDT rer 
+ P°+b;,p+cy 


re Ha 


(P— Pi4 NU 


C!p+ D: 
+ > + , 
(P? + bn41p + Ch+1) 
où les coefficients A, B8, Cy: Co; D); y A2: B8, C, Ca, D,, Ds 
sont des constantes, et en outre À, Bh, C:, Cs, D, C5 sont des 
fonctions homogènes linéaires des données initiales. Les fonctions 
r(t) et A Mur comme suit: 


r(t)— S Age! + 3 HW a Dm im 


a=1i 


s 
_- —t 2D,—C;b, 
+ D € [C cos Pt + — 5 — sin Prt | + 


To 
 Ÿ 08 gra x 
0 
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X SID Pa+r (Te — T1) - -. SIN a+ ( — Tor) dT . . . Toi + 


DD de t To Ta 
É —— [ | EE Î Si Px+1T1 - SN Ph+s (To — TU) . 
Ph+1 0 0 0 


. Sin Px+1 (£ — Ts) dr; é a dt | rss 
l B; 
Ha AePa' +2 (B— NT 010 ePisi! + 


ai 
k es 
- > ; : — C! B | 
+2e ? [ Cy cos Prt+ — PT sin Pyt |+ 
_ be s ARR Vo-1 Te 
+e à: EEE me | cos qu X 
Pt Ù 0 . 
X SiN Pa+1 sd ... SÈN Qué (É — Tor) dE . .. dt + 
t To-1 Ta 
Coby+ \ ; ; 
20° MERE | -.. | SIN Pa+1Ta -SÈN Pr+r (Ta — T1) . . - 
R+1 0 9 ü 


. sin Pau (— To) dt dieu | Fo 


où 


R+1 


Pr+1 = Ch+1 — 


En portant les expressions de r (t) et s (f) dans la formule (2.54) 
on obtient la solution générale de l’équation différentielle étudiée. 

2. Soit un système d'équations différentielles linéaires à coef- 
ficients constants a;,, dont les termes f, (t) sont des fonctions don- 
nées du temps 


dz 
= ui + GeTae + - .. + dinTn + fa (): 


+ = Ao]Ty + Goolo + .. + ZonTn + Î2 (£), 
(2.64) 


_ = Any21 + donTe +... + AnnTn + fn (t). 


h —0394 
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En multipliant chaque équation du système par e-?{ et en intégrant 
par rapport à t de O0 à co, on obtient 
(au — P) Ti (P) + Gie72 (D) +... + GinTn (P) = 
= — [fi (p) + 21(0)], 
GT 1 (P) + (Ass — P) T5 (P) + . .. + Gent (P) = 


Gn1Ti (D) + AneT° (P) + . . . + (Ann — P) Zn (P) = 


= — [fr (P) + zn (0)]. 

La résolution de ce système donne 
A 
TR = ND (= 1; 2:50); (2.66) 
où 
(&11 — P) do in 
ben ee 
Œn1 An2 ... (aan — p) 


est le déterminant principal du système d'équations (2.65), 


A: = 2 (—1)*"f;Aun (p) 2 (— 1)" 2; (0) An (p), 


et 
(a — D) di12 ce.  ,h-1  i,h# Gin 
da (ass — D) Ao,h-1 se, k+1 on 
Ah (p) — di-1,1 dix, ° Œj-1,h-1 Ci, k+1 --. Œinion 


i+1,2 --. its, kh1 ita, ki ce. Cipan 


Œn1 ne ... n, k=—1 ah, R+1 :--: (aan — D) 


est le mineur obtenu par suppression de la i-ième ligne et la k-ième 
colonne. Donc, la formule (2.66) peut encore s’écrire 


TR (P) = à fi Cp) Di (p) + à zi (0) Di (p) (2.67) 


* — (—1\ith# Air (P) 
Da (p)=(— 154 ARC 
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sont des fractions rationnelles de p dont le degré du numérateur 
Ar (p) est inférieur d’au moins une unité à celui, égal à n, du déno- 
minateur À (p). Pour décomposer D, (p) en fractions simples, il 
faut connaître les racines de l'équation À (p) = 0. Une fois qu’on 
a déterminé la fonction x; (p) à partir de (2.67) et ensuite x, (t),on a 

n tt n 

(= | fi (x) Din (é—) dr + S 23 (0) Dan (+). 

i=1 0 i=1 
Cette méthode peut être appliquée à des systèmes d’équations liné- 
aires d'ordre supérieur 


di d 
D. {on TE + un DE + cyan } —= J\t (V=T, 2, :.:,n), 
=1 


Zh(0)— an. Th (0)=Bx (k—=1, 2,...,n). 


Appliquons la transformation de Laplace aux deux membres de cette 
équation. On obtiendra l'équation suivante: 


à (avr D? + byrp + Cyr) Ti (P) = f$ (P) + à [(avrP + dur) Gr + GvrBr]. 


D'où l’on déduit z* (p) par les méthodes de l’algèbre linéaire. Sans 
exhiber de formules générales, considérons à titre d'exemple le 
système de deux a __ du second ordre 


d°xz { 


de et An r Aire > = 1121 + AuoTe + fi (4), 
Ars , (2 68) 
TR + An — _ 1. = A91T1 + GooTo + fo (t), 


où 7, et x, sont les de: inconnues de t, f1 (t) et f, (t) des fonc- 
tions « perturbatrices » du temps données, A1, A9, Aots Â9es Gris 
Gje, Go) ee des coefficients constants. 
En multipliant chaque équation du système (2.68) par e-7' et 
en intégrant par rapport à { de O0 à oc, on obtient 
p°2i — pti (0) — x; (0) + PAuTi — 41121 (0) + A12 PTS — 
— Ajo%e (0) = auf + 1272 + fi, 
P?z5 — PTe (0) — x, (0) + pAoyTT — Ao1T1 (0) + A22 PTS — 
— A30To (0) — GT + GooT5 + f5. 


Après réduction des termes semblables, il vient 
(P°+ PA — @is) 2f + (A12P — Que) 22 = 
= fi + (P + Aux) T1 (0) + A12%2 (0) + x; (0), 
(A21P — es) 2 + (P?+ PÂ2e — G2e) 22 = 
= f5 + Ao1T4 (0) + (p + Az) 2 (0) + x: (0). 


(2.69) 


4® 


52 TRANSFORMATION DE LAPLACE [CE I 


Donc, 
Zi (k—1, 092), (2.70) 


où À est le déterminant principal du système d’équations (2.69) 
(P°+ A11P — Gu) (A19D — Gi) 
(431P — @21) (P?+ A22P — G92) |” 
Les quantités A, et À, se calculent à l’aide des formules 
Ai fivii—fovia+ ÔT, A = f2Vèe — fi + 05, 


Vi1 = D? + Ao2P — Goo,  V12 = AyeP — Go, 
V1 = AnP—Gns V2 = P°+ AuP—ü@u; 
6 = 21 (0) [(P + An) Yt1 — AnVi2] + 
+ Lo (0) LA1oYT1 — (P + A9) Vi2] + Zi (0) VT1 — Z2 (0) VÉ2, 
05 = Zi (0) [A21V22 — (P + Au) Y51] + 
+ Za (0) [(P + Ào2) V$2 — A12V21] + Zi (0) V51 + z2 (0) V5. 


Introduisons les notations: 


* R+n Vén k ô$ æ 
Pin =(—1) TA Di = (k, n = 1, 2). (2.71) 


Les formules (2.70) s’écrivent alors 
TR = fil + lie + Di (k=1, 2). (2.72) 


Les quantités l'é,, D? (k, n = 1, 2) sont des fractions rationnelles 
du paramètre p dont le degré du numérateur est inférieur à celui du 
dénominateur, qui est égal à quatre. En décomposant ces fractions 
en fractions simples, on déduit aussitôt 


t 
Th &= | LA (T) Lu (CT) + far) Tuo(t—t)ldt+D:(t) (2.73) 
0 


(k=—1, 2). 


3. Il existe des classes d'équations différentielles dont on peut 
représenter les solutions par des intégrales de Laplace, où la variable 
indépendante figure sous le signe somme comme paramètre. Consi- 
dérons l’équation { 


(an + bat) 2 (t) + (ans + Dnat) 2 (E) +... + (Go + bot) z (t) = 0. 
5 # (2.74) 
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Soit 
z(t)— | ePty (p) dp, 
l'intervalle d’intégration n'est assujetti à aucune condition pour 
l'instant. Ceci étant 
ab (t)= À ertpho(p) ap, te%(t) = | teP'phv (p) dp= 
Û d 
= [eP' pv (p)] — | eP° D [pv (p)] dp. 


En portant ces expressions dans (2.74), on obtient 


| ep 16 expo (P)— À ba LP*v (p) | dp+ 


k=0 


+ ZX dx le” pho (p)1 = 0. (2.75) 


k=0 


Cette équation est vérifiée si l’accolade de (2.75) s’annule, et on 
obtient alors une équation différentielle du premier ordre pour la 

détermination de la fonction v (p). Le second terme doit être égalé 

| zéro moyennant un choix convenable de l'intervalle d’intégration. 
oit 


"(+ (a+ b+t)z (i) + ax (t) = 0. 
La transformation de Laplace, appliquée à la fonction zx (t) = 
D | eP‘ v (p) dp, nous donne pour la détermination de v (p) l'équa- 
tion suivante: 
v'(p)(+p)—v(p)(p(a+b—2)+a—1l=0. 
D'où 
v (p) = (p + 1)°-1 p°-4. 


Pour l’annulation du second terme de (2.75), on a 
Le?! (p+ 1) p° 1 (p2+ p}le = 0, (2.76) 


où æœ et B sont les extrémités de l'intervalle d'intégration. Pour 
fixer les idées nous supposerons a >> 0, b > 0. La condition (2.76) 
est remplie pour &« = —1, 6 = 0. Donc, la première intégrale de 
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l'équation (2.74) s'écrit 
0 
ni (t)= | ep (p+-1)-t p""t dp. 
21 


En posant B = 0 et &« = —o, on obtient la seconde intégrale (du 
moins pour les t >> 0) 


0 
mt | é4(p+ 1) pa 
Dans de nombreux cas, on est obligé de choisir le chemin d'inté- 
gration dans le plan complexe. 
Considérons l'équation 
tr" (t) + Qnx' (t) + tx (t) = 0. 
Comme précédemment on trouve 
v (p) = (p° + 1)". 
La condition 
Le?! (p2 + 1)" = 0 
nous fournit @œ = —i, Pf = <+i. Donc, 
+i 
(= | et(p2+ 1) ap. 


—1 


En posant p = iu, onjobtient 
T1 (t) = à eut ({—u?)""1 du. 
1 
Séparons la partie réelle de la partie imaginaire : 
Zi(t)—à | cos ut (1 —u?)""1 du — sin ut (1—u?)""t du. 
1 =1 
La seconde intégrale est nulle, puisque l’intégrande est une fonction 


impaire. On a donc 
1 


a (t) = | cos ut (1—u?)""! du. 
1 


On obtient la deuxième intégrale de l'équation (2.74) en posant 
— —00 et Bf — +i ou —i (t > 0). En intégrant de —œ à 0 et 
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ensuite de 0 à à, il vient 
0 


1 
2 (= | et (p24+ 1) dp+i | eitt(1—u2)"-t du = 
0 


0 1 

— | e“! (u2+ 1)" 1 du + i ( cos ut ({ —u2)""1 du — 
, 
0 


1 
— \ sin ut(1—u?)""t du. 
0 


La partie imaginaire est donc égale à 5 z1 (t); la partie réelle doit 
également être solution, i.e. 
0 1 
La (t) = | e“t(u2+ 1)" du + | sin ut (1—u°)""1 du. 
0 0 0 


Dans {217] on trouvera de nombreuses questions relatives à ce 
point. 

&. Cette méthode peut être utilisée pour résoudre les équations 
différentielles aux dérivées partielles en mathématiques appliquées. 

Soit donnée nu 


Pau (a, y 2) + b(2 y 9 +e(su du f(vnt, CT) 


« d , » 
où V?u —= = ne EE = + a est l'opérateur de Laplace, (x, y, 2) 


un point d'un certain domaine et t le temps (t => 0). La condition 
aux limites s'écrit 
(4) 
œ(z, y,z)u+B(x, y, 2) = =p(x, y, 2 t), (2.78) 


. Ou 2 PR | 
où,— représente la dérivée normale. Par ailleurs, on donne des con- 
ditions initiales à l’intérieur du domaine, par exemple 


limu(xz,y,z,t)—=uo(x, y, 2), (2.79) 
t-0 
lim À u(xz,y.2,t) =u(x,y,2). (2.80) 
t-0 


Multiplions l’é équation initiale (2.77) par e-Pt et intégrons par rap- 
port à t de O0 à co. Supposons que les intégrales 


| ePlu(z,y,z,t)dt, | el — : u(z,y,z,t) dt, etc. 
Û 0 
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existent. De plus, 


[ ePiy2u dt — V? [ eP'u dt. 
0 0 


Sous les hypothèses faites sur les propriétés de la fonction inconnue 
u (x, y, z, t), on obtient à partir de (2.77), (2.79), (2.80) l'équation 
V'u*(p) + la (x, y, 2) p° + b(z, y, z) p + c(x, y, z)l u* (p) = 
— a (x, y; z) [PU (x, y, z) — WU (x, y; 2)] + 
+ b (x, y, z)u (x, y, 2) + FF (x, y, z, p). (2.81) 


La condition aux limites (2.78) se transforme en la suivante: 


œ(x, y, zu*(p)+B(x, y, 2) = (x, y,z, p). (2.89) 


Une fois qu’on a déterminé u* (p) à partir des équations (2.81) et 
(2.82), le problème consiste à tirer u (x, y, z, t) à partir de l'égalité 
u*(z,y,2, P) = | ePlu(z, y,z,t)dt. 

0 


Siu* (x, y, z, p) figure dans la table des formules, on a la solution 
aussitôt, sinon on peut l'obtenir à l’aide du théorème d’inversion 
Y+i00 

| Mu (x, y, z, À) dà. 


Y—100 


u(xz,y,2,t) — 


271i 


On rappelle (cf. $ 1) que pour calculer cette intégrale, on transforme 
l'intervalle d’intégration en un contour fermé adéquat auquel on 
applique le théorème des résidus. On remarquera que lorsqu’on forme 
l'équation auxiliaire et ses conditions aux limites et lorsqu'on 
déduit la fonction u (x, y, z, t) à partir de u* (x, y, z, p) à l’aide du 
théorème d’inversion, on impose certaines hypothèses aux propriétés 
de la fonction u (x, y, z, t). Les hypothèses telles que la possibilité 
de permuter d’un côté la transformation de Laplace, et de l’autre, 
les opérations de différentiation et de passage à la limite, ou encore 
que la solution doive être de forme déterminée, ou encore qu'elle 
soit développable en série, etc., ne sont pas restrictives du point de 
vue physique. Par ailleurs, la méthode de résolution indiquée peut 
être utilisée formellement si le résultat obtenu vérifie l'équation 
ainsi que les conditions initiales et aux limites. 

Résolvons quelques problèmes concrets. 

1°. Considérons l’équation à une dimension de la chaleur 


du d°u 
DT (2.83) 
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et supposons que sur l'intervalle 0 < x < L est posé le problème 
aux limites avec les conditions aux limites 


ule=o == Pit), Us: = Po (t) (2.84) 
et la condition initiale homogène 
u |=0 = 0. (2.85) 


Remplaçons la fonction u (x, {) par sa transformée de Laplace 
u* (x, p) = | e Plu(zx,t) dt (2.86) 
0 


qui sera la nouvelle fonction inconnue. En appliquant aux deux 
membres de (2.83) la transformation de Laplace et en supposant 
que dans (2.86) il est possible de dériver par rapport à x sous le 
signe somme, on obtient 

PE = put (2, p). (2.87) 


Appliquons la transformation de Laplace à l’équation (2.84) 
U° [oo = Pf (D), U° fret = P3 (P);, (2.88) 


où 


qi (p)= [eq (tt (k=1, 2). 
0 
Des équations (2.87) et (2.88) on déduit que 
u® (x, p}=@f (p)oi (x, p)+ à (p)o? (x, p), (2.89) 
où 
. _sSh{—z) Vp  » _ shzVp 
@Ÿ (Z, p) = ship ? (PE. 


Les fonctions of (x, p) et «2 (zx, p) sont les transformées de La- 
place des fonctions 


UE Lot 109608 [i—z 
rl Tr rl Fr 


respectivement, où 


Da(v,t)—= D, eninr-nmt 


n=—= — © 
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est la fonction de Jacobi. A partir de (2.89) déterminons uw (x, à) 
à l’aide du théorème de convolution 


u(x, t)= 75e, +=) P1 (T) dT+ 
0 


1 
1 () l— t— 
T ad, +) (0 dr. ES) 
Considérons maintenant l'équation non homogène 
ôu Gâu 
7H C2 oi 


avec les conditions initiale et aux limites homogènes 
u(z,t)l1=0=0, u(x,t)l=o=0, u(z,t)|xu = 0. (2.92) 
Supposons que 


f° (a, p}= | ef (2, 1) dt. 
0 


Appliquons la transformation de Laplace à (2.91) et (2.92) 


ce = pu*— f° (x, p), (2.93) 
u° (0, p)=u*(l, p) =0. (2.94) 


I1 est immédiat de vérifier que pour l'équation (2.93) avec les con- 
ditions aux limites (2.94), la fonction de Green est de la forme 


SGD VpshzVPp pour z<E, 

| . VrshiVp 
PES PE UD VPhEVE Dour 25E 
VrshiVp un 


et la solution de l'équation (2.93) vérifiant les conditions (2.94) 
est de la forme 


l 
u®(z, p= | ve, & p)f'(e, PE (2.95) 


0 


La fonction y* (x, £; p) est la transformée de Laplace 


v(z, 6 = [?s cn F)—*s | LE F)] 
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La formule (2.95) entraîne donc 
L 


t 
un(e, D= [Œ ve, Et-DfE Dr. (2.96) 
0 


0 


2°. Soit à trouver une fonction satisfaisant à l’équation (2.83) 
et aux conditions 


u(z, 0)—0 (z>0), u(0, t)—f(t). (2.97) 


Comme précédemment, en appliquant la transformation de Laplace 
à l'équation initiale (2.83) et en tenant compte des conditions (2.97) 
on obtient 


# 
EC = pu (x, p), u*(0, p)=f"(p). 
Mettons u* (x, p) sous la forme 


u*(xz, p)= Ce-<VP + CiexV?. 


Comme u* (x, p) est bornée lorsque z —+ ©, il vient 


: == V? 
u*(z, p}=f"(p}e-<Vr = pft(p) ——., 


d’où l’on déduit, moyennant le théorème de convolution 


= OR en 
u(x, {)= PCA | Ë e dt 
l—T)" 
2 , 
RE ff)" 


Il est aisé de voir que 


o œ 
z, 0) =0, 0, t)—=f(t = -$ dé = f (t). 
u(z, 0) =0, u (0, )— 107 | < ë* Œ = f (+) 
3°. Posons le problème suivant pour l'équation de la chaleur 
(2.83). Supposons que 0 < z << ©, 


u(z, 0) =u,, u’(0, t) — hu (0, t) (hk = const). (2.98) 
La transformation de Laplace ramène l'équation initiale (2.83) 
et les conditions (2.98) à la forme 


du* du* 


Asie » _ Achat == » 
ms = Pu* ue, | hu*. 
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Comme précédemment, la solution u* (x, p) étant bornée lorsque 
z —oo, il vient 


u* (x, p)= + +ce* Vr, 


Lee V5 A (ES + 6) = Op 
d’où 
* #0 — + 
17 (x, P) = : (1 AT d Vr)= 
= 0 4 — -xV? 0 —_— -*V?, 
pie 1+ Vr(VrEN) 
Comme 


7 {uert(— y: | 3 = (1—e-<Vr), 


HAE der ne Jp} 
en se servant de l'expression 
"(V?) ” 1 À en -E 
2 = — 0. — — e st dr 
Vr  Vr(Vr+t) Vu ): 
on obtient 
— f uo_ —h(tT-x)- _ | 
u (x, t)=wer (—= ER VARE 7= ee 


49. Soit à trouver une solution de l’équation 


2u du ee. ou 


Ge 2x — —cu =0 (0<z< OO, t > 0) (2.99) 
telle que 
u(z, 0)=0, u;(x, 0)=0, u(0, t) = (t). (2.100) 
En appliquant la transformation de Laplace, on obtient 
PE — (ap? + bp+c)u* =0, (2.101) 


_. (z, 0) = p* (p). (2.102) 


L’équation homogène (2.101) se résout de façon standard. La solu- 
tion étant bornée lorsque z —> , il est immédiat d’établir que la 


transformée de Laplace de la fonction uw (x, ft) est la fonction 


u* (x, p) = p* (phe-xVar#tbrtce, Si ac — (+) = 0, la fonction 


2 


*) L'intégration a lieu entre z et co. 
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u*(z, p) s’écrit 
b 
—x[PVa+—— 
u*(z, p)=p"* (p)e Pavel, 
d'où 
nn 
u(z,t)=e 2Va p(i—zxV a). 
Si æa—=ac— (3) #0, alors compte tenu de la relation 


b « 
e-xVaPpi+bp+e = € 2Va e-xVa == 


= co __ VE VF 
pe | Pie 20 — di, 


4 Ve 
on aura : 
0 lorsque 0OSt<z Va, 
b 
e 2Va p(t—r V a)— 
u (x, t) — = { D, J, ve V/ T2 — az? 
_:)/ © LE a DE = RE 


| lorsque t>z Va. 
5. Considérons l'équation de Volterra de deuxième espece 


Î 
PO=f (+ | (+) pr) dr (2.103) 
Û 


que l’on rencontre dans diverses branches scientifiques. Supposons 
que toutes les fonctions figurant dans cette équation admettent 
pour transformées de Laplace respectivement: 


p(p)}=ÆLlp(t)], f*(p)= LÉ (FE), k*(p)=Z£ [k()]. 
En utilisant la convolution, on a 
p* (p)=f° (p)+E* (p) p° (p), 


nf) 
PO= RG 


d’où 


et 
; o+i 
PO= | v(p)etdp. 


O—ico 
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S'agissant de l'équation (2.103), tous les noyaux itérés dépendent 
de la différence £ — +, donc la résolvante dépendra également de la 
différence ?t — +. 
L’équation de Volterra de première espèce 
Î 
1 = | AG) pr) dr 
0 


se résout de façon analogue." Cette méthode est valable aussi pour 
les systèmes d'équations intégrales de Volterra de la forme 


n { 
p=f(+S | Ant—r)qimér (=1,2,...,n). 
k=1 0 
En appliquant la transformation de Laplace aux deux membres de 
ces équations, on obtient 


gt (p)= ft (p)+ 2 Kù (p) gà (p) (i=1,2,...,n). 


La résolution de ce système d'équations du premier degré nous donne 
p* (p), et la solution du système d'équations s’écrit 
1 é' à t 
= | qi (p) «7 dp. 
ci 
Citons quelques exemples. 
1°. Equation d'Abel. En étudiant le problème de la tautochrone 


{ 
| MANEURET TES (0<a<1) (2.104) 
n ET 

Abel fut le premier, en 1826, à obtenir une équation dans laquelle 

la fonction inconnue  (t) figurait sous le signe d'intégration. 
Supposons que * (p) = £{ (t)], £{f (t)] = f* (p). Moyennant 

la transformation de Laplace, on ramène l'équation (2.104) à la 

forme 


ro) = f° (p). 
D'où 


1-a 
sn P f*{p) ___ f{) 1 pf* (p)— f (0) 
PO Fa) TU -S tra 


Il est aisé de trouver 


t 
- 1 f (0) j'(o dr 
PO= TT l'A +l ave] 
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Comme T(1—«) T'(œ) — il vient 


sin ar ? 


t 
sinar [0 jf’ AUS dt 


>] 


20. Considérons l’équation intégrale à noyau logarithmique 
{ 


| (rt) ln (t— 7) dr = j (t), (2.105) 


0 


où œ (+) est la fonction inconnue. Supposons que p* (p) = £(œ (t)}, 
fs (p) = Æ£{f (dl. Utilisons l'égalité 


C the-ch 1 
£ (| Tr (&+1) dk) = p(np+a) * 108) 

Sachant que 

£ (Int) — —+(C+in p), où C'est la constante d’Euler, on 
obtient 

L 
— p* (p)-(C+ In p)=f (p). 
Donc 
*(p)}= — pf* (p)____ p'f*(p)—f" (0) _ ___f’(0) 
PPT MpHO —— pÜnp+O)  PUnp+C) 


Compte tenu de (2.106), il vient 
t 00 £ co 
r@=—] f” (x) dx A dk —f" (0) Fr dk 
3°. Considérons l'équation 
sin { — Jo(t—T) p(Tr) dr. 
0 

Si on lui applique la transformation de Laplace, on obtient 

re VO) PO, p()= Jo (D. 


La dernière formule entraîne 
t 
sint= | Jot—T) Jo(r) dr. 


0 
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À noter que ce qui vient d’être dit se généralise aussitôt à toute 
équation intégro-différentielle de Ia forme 


ap (6) + ap? (4) +... Han () + 
7 t 
+2 [An (7) pf (1) dr f (0 
j=0 0 
et aux systèmes de telles équations. 


$ 7. Transformation de Mellin 


La transformation de Mellin 


F (s) — | f (t) tidt, s— GO + iT, (2.107) 


0 


est étroitement liée aux transformations de Fourier et de Laplace. 
La transformation de Mellin peut être appliquée avec succès 
à la résolution d’une classe de problèmes harmoniques plans dans 
un domaine aréolaire, à des problèmes de la théorie de l’élasticité, 
ainsi qu’à l’étude des fonctions spéciales, à la sommation des séries et 
au calcul des intégrales. Les théorèmes relatifs à la transformation 
de Mellin se déduisent des théorèmes correspondants des transforma- 
tions de Fourier et de Laplace moyennant un changement de varia- 
bles. Citons quelques uns d’entre eux. 
Théorème 16. Supposons que 1°-!f (t) E LI0O, +oœl, où 
la fonction f (t) est à variation bornée au voisinage du point t = t. 
A lors 
| ; o+i. 
HER O0 LL Vin ( F(s)t" ds, (2.108) 
L À 00 
oi 
où la fonction + (s) est définie par (2.107). 
Théorème 17. Supposons que F (o + iu) E L] —o, +ol 
et qu’elle est à variation bornée au voisinage du point u = t. Alors 


À 

LE (o+i(t+0)+F (o+(t—0))] = lim [ ft) +1 de, (2.109) 
— 00 1 
+ 


où 
G+ioo 
FO = [ F(s)t"ds. (2.110) 


OC—ico 


6 7] TRANSFORMATION DE MELLIN 65 


Supposons que # (s) et G(s) sont les transformées de Mellin des fonc- 
tions f (t) et g (t). I] vient aussitôt 


R-Hio0 h-bioo 


2 | FOcu-98-7 | Gas à | oe = 
k—100 k-- 100 
; co k+i00 co 
= (ro | ea-Jmase[ just. (2119 
0 h—io 0 
De façon analogue 
k-+Hi00 co CEST. 
a | soc feu | ras 
k—1i00 Û Rk— 100 


= Tews(2)#. (2.112) 
0 


Théorème 18. Supposons que 
t*-1f (HE L]0, + oo[ et G(1—k—ix)E L]—o0, +of, 
ou que 


F(k+ir)EL]—oo, +oo[ et g(t)ELI0, +of. 
La formule (2.111) est alors valable. 
L'analogue du théorème de convolution est le 


Théorème 19. Supposons que t“f (t) et t*g (t) appartiennent 
à l'ensemble L ]0, +oo et soit 


o= 10e (2) E 
0 


Alors la fonction t“h (t) € L J0, +ool et sa transformée de Mellin 
est la fonction # (s) G (s). 


Pour la transformation de Mellin appliquée à d’autres classes 
de fonctions et certaines de ses applications, voir [238], [180]. 
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CHAPITRE Ill 


TRANSFORMATION DE BESSEL 


On appelle transformation de Bessel une transformation intégrale 
de la forme 


f(Q)= | 7 K Gt) de, 
0 


où Æ (z) est une fonction de Bessel. 
Les transformations de Hankel, Meijer, Kontorovitch-Lébédev, 
etc. sont des transformations de Bessel. 


$ 1. Transformation de Hankel 


Les formules du type (1.14), (1.15), (1.16) qui développent une 
fonction arbitraire f (t) en intégrale de Fourier présentent un intérêt 
notable dans nombre de problèmes de mathématiques et de physique. 
Parmi ces développements, mentionnons le développement suivant 
des fonctions cylindriques, connu sous le nom d’intégrale de Fourier- 


Bessel : 


f(a)= | Jo(œu)udu | (0 J,(ut)tat (0<Lz<+o), (3.1) 
0 0 


où J, (x) est une fonction de Bessel, 
VS 


Théorème 1. Soit une fonction f (x) à variation bornée sur 
tout intervalle fini 10, RI et 
= 1 
| [f(xz)1z° dr 00. 


0 
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: 1 
Si V>—-5,0na 


L{f(&+0)+f(x—0)}= ( J, (zu) u du [ f(t)Js(ut)t dt. (3.2 


0 


Aux points de continuité, on a la formule (3.1). 

La validité de ce développement est démontrée sous d’autres 
hypothèses dans [238], [70], [42]. 

On appelle transformation de Hankel l'intégrale 


f(u)=68,1f(01= | FH, (ut)dt (O<u<+o). (3.3) 
0 
Le développement (3.1) entraîne la formule d’inversion 
1O = 1f(u)1= | fS(u)Ji(ut}udu (0<t<+0o). (3.4) 
0 
Observons que si la fonction f (t) est telle que f (t) = © (t£) lorsque 
t0,a+v+2>0etf(t) = O (t#) lorsque t —> 00, B + + < 0, 
alors l'intégrale (3.3) est convergente [42]. 
I1 existe un développement analogue à (3.1), soit 


| | 
f(x)= | H, (zu) (zu)° du | Y, (ut) (ut)£ f() dt, (3.5) 


0 0 


où Y , est une fonction de Bessel du deuxième ordre, H, une fonction 
de Struve. 

La formule (3.5) sert de base à l’introduction de la transforma- 
tion intégrale correspondante [222], [238], [235], [85]. 

Citons le lien existant entre la transformation de Hankel et les 
transformations multiples de Fourier. 


Soient 
{ œ © | | 
f(x, y) ue . Î [ eikx+iou fs (4, &) dÀ du), 
l © œ | | 
PO open | À etais j(, pds du. 


En passant aux coordonnées polaires moyennant les formules z + 
+ iy = rei® et À + iw — peï* on établit que les fonctions f (x, y) 


5% 
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et f* (2, w) sont liées entre elles par la relation 
oo 27 


14 k 
u(r p=s— | [ eirp cos (Ou (p, xp) p dp dy. 
7 p=0 ÿ=0 
Posant 
v'(p, +) = i-"x (p) er 
et comme 
=: 27 
Jh (2) = _ Î eiz cos 0+in8 qÿ, 
0 
on aura 
u(r, q) =ei"® | PJ n (rp) x (p) de. 
0 
Posant 
u(r, p)=© (r) ein, 
où 
D (r)= | pJn (rp) x (p) dp. 
0 
on aura 


v(p, D)= éteint | Ja (—rp)O (rrr, 
0 
x (P)= | Ja (rp) D (r)rér. 
0 


Enonçons encore quelques propriétés de la transformation de 
Hankel. 
(7 


o 1 
1°. If (a =-x (>). 
2°. Une intégration par parties donne 
d 1 4 È cs 
Ge [SE + f]= 0258, 1f (1. 


Ceci étant, on admet que 


1 
ri 20, 2420, 
di | (t— 0), : sh (t — co). 


Lu (9 --0 
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3. Egalité de Parseval. Soit 
Perf (= fu), SP,]g(t)] = g$ (u)- 
Alors 
Lee) . . œ { 
jufe du te (v>-S). 
0 0 
Les conditions de validité de la dernière formule sont énoncées dans 
[238], [42]. 
4°, Asymptotique [42]. Si f (t) = O (12) lorsque’t +0, & + v + 
+20, f(t =O(#) lorsque t +00, B + 5 < 0, alors f (u) 
(v > —1, u > 0) existe et, d’autre part, 
fi(u)=O(u%) (u—>0), &« >min(v, —6—2), 
f5(U)=O(UP) (u— 00), P&max(—+, —a—2), 


2 ? 
La formule 


_ Pu (t) u du 2 
10e] Tire ren | / (9) pu (9 dr (@<t<+o), (3.6) 


où 
pv(t) = Jy(au) Yo (ut) —Y,(au)J,(ut) (v>—+) 


est une combinaison linéaire de fonction de Bessel de première et 
de seconde espèce d’ordre v qui généralise le développement (3.1). 
Le développement (3.6) a lieu si la fonction f ({) est continue par 
morceaux et à variation bornée sur tout intervalle fini Ja, RI et 
l'intégrale 

co 1 


[IFOIÉ dt < oo. 


a 


Lorsque a —+0, la formule (3.6) se change en (3.1). 

Comme dans le cas de la transformation de Hankel, le développe- 
ment intégral (3.6) conduit à une transformation intégrale que l’on 
appelle transformation de Weber. Il y a intérêt à considérer la trans- 
formation généralisée de Weber 


f" (u) = | tC, (ut, au) f (à) dt, 
où ‘ 
Cy (2, w) = J, (2) IPY, (w) — QuY ,41 (w)] — 
— Ÿ, (2) (PJ, (w) — QwJ,: (w)l, 
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et P et Q sont des constantes arbitraires [86]. Les transformations 
de Hankel et Weber peuvent être utilisées avec succès pour résoudre 
les problèmes aux limites posés pour l'équation de Laplace et de 
Helmholtz et certains problèmes de la théorie de l'élasticité et de la 
chaleur. 


$ 2. Transformation de Meijer 


La transformation intégrale de Meijer joue un rôle important 
dans la résolution des équations différentielles de type Bessel. Elle 
est définie par la formule 


de ; 
F7 2 À K(s1) (9° 1 (D de. (3.7) 


0 


où Æ, (t) est une fonction de Mac Donald. La formule d'inversion 
s'écrit 
B+i2. 
1 ; 
lim 


Tite \ ; I y (és) (ts)? 7 (s) ds. (3.8) 


f(t)= 


Les égalités (3.7), (3.8) découlent immédiatement du développement 
correspondant d’une fonction arbitraire en intégrale de Fourier 
généralisée. 

Théorème 2. Soit f (t) une fonction de la variable réelle 1, 
0O<t< +oo, intégrable sur tout intervalle fini OT, <I<T: 
et à variation bornée au voisinage du point t — x; supposons par ailleurs 
que l'intégrale 


le-#|f{ld, B>a>0 
0 
est convergente. 
Sous ces conditions on a le développement 


B+ix 1 © 1 
LGTOHfGE-0 1 lim | I, (xs) (zs)“ ds ( K, (st) (st) f({t) dt. 
2 Ti 1.00 nn l 


(3.9) 
La fonction XÆ, (z), on le sait, est une fonction paire de v, donc 
Ja dernière formule peut encore s’écrire 
x B+ix 1 
. 2 
Zi Jim jh {I, (xs) + ZI, (zs)} (zs)” ds X 
i 


f(z+ OL f(z—0) 
2 


co 1 
x À Æ(st) (st)? f (6) dt. (3.10) 
0 
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En particulier, puisque 
1 


11(=(<) shz, 11@-(+) chz, 
. “a 


B+ià co 
-L0 1- _—0 Â à xs —s 
f(z+ À f (z ea im | € si à tf(t)dt. (3.11) 


Théorème 3. Supposons qu'une fonction f (s), analytique dans 
le demi-plan Res > a > 0, vérifie les conditions: 
10. Pour tous les t > 0 et f > « existe 


B+i2 1 


LS 
L 2 


lim | TI, (tz) (tz) 1 (z) dz, 
ep 


où —+<Re ST, 
la convergence étant uniforme en t (0 <t< b). 
2. L'intégrale 
FE 17e) 
| F1 [dl 
B— io 
est convergente. 
3. La fonction f (s) est bornée dans le demi-plan Res > B, i.e. 


| (s)1< 4, 

où À est une constante positive ne dépendant pas de s. 

49. Existe 

lim f(z+iy)=0, 

où la convergence est uniforme en toutes les valeurs réelles de y. 

Sous ces conditions, lorsque Res >, on a l'égalité (3.7), où 
la fonction f (t) est définie par (3.8). : 

Théorème 4. Soit une fonction g(s) telle que 


g(s)=f(s) + 2 cisi, 
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où f(s) remplit les conditions du théorème précédent, Re (8,) < 0 
(i=1,2,...,n). Alors 


… 1 
&=V À [ K(st)(st)?g(tdt, Resp, 
0 


où 
B+iA L 
S 
g(t)=- TE . j rt) (ts) “ g(s) ds. 


La convolution et le calcul opérationnel ont été construits pour 
la transformation de Meijer d’après le schéma de Mikusinski (cf. 
chapitre V, $ 16). 


$ 3. Transformation de Kontorovitch-Lébédev 


1. Les transformations intégrales dans lesquelles l'intégration 
s'effectue sur l’ordre des fonctions de Bessel jouent un rôle très 
important dans la résolution de certains problèmes de physique 
mathématique. Cette forme de transformations intégrales a été 
étudiée pour la première fois par Kontorovitch et Lébédev en 1938 
[120]. Elle porte aujourd'hui le nom de transformation de Konto- 
rovitch-Lébédev. Elle a été appliquée avec succès à la résolution de 
nombreux problèmes intéressants [87], [134], [179], [180]. Dans 
les transformations de Kontorovitch-Lébédev, le développement de 
type intégrale de Fourier 


zf (2) = + | Kis(æ)rshrrdr ] Ka (E)f(È) dE. (3.12) 
0 


où X,(z)est une fonction de Mac Donald. x >> 0, f (x) une fonction 
arbitraire continue avec sa dérivée et telle que z*f (x), zf (7) € 
€ L J0, +oo[ joue un rôle fondamental. Le développement (3.12) 
est valable pour des fonctions appartenant à une classe plus large 


[134], (1371. Soit 
F(x) = | f(x) Ke (x) dz. (3.13) 
Ù 


L'intégrale (3.13) définit la transformation de Kontorovitch-Lébédev. 
De (3.12) il vient immédiatement la formule d’inversion 


f(x) = A Kis(z)tshatF (x) dx (20). (3.14) 
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Sous une forme plus symétrique, (3.13) et (3.14) s’écrivent 


F(9= | 14) de (OST <+ 0), (3.15) 
0 
je À FC RTE His) HA de (O<z<+o). (346) 
0 


Quelquefois (3.13) et (3.14) se réncontrent sous la forme suivante : 


F (x) = <sh “| (x) —=<—<— Ass Re dx, (3.17) 
f(z) = | F(x) Kun(x) dr. (3.18) 
0 


Les formules analogues aux formules de Parseval de la théorie 
des séries et des intégrales de Fourier sont très utiles au calcul de 
certaines types d'’intégrales définies. Enonçons deux théorèmes 
[35], (1361. 

Théorème 5. Soit g (x) une fonction réelle arbitraire telle que 


à 
1) g(x)z * EL]0, Ho, 2) g(x)AL:10, +oof, 


ct | See 0) A dx. (3.19) 
Alors 
Î [G (x)]2 dt — | HO (3.20) 
0 0 


Théorème 6. Soient g,(x) et g;(x) des fonctions réelles 
arbitraires, vérifiant les conditions 1) et 2) du théorème précédent, 
G; (rt) et G, (rt) Les images respectivement de g, et g, par la transfor- 
mation (3.19). 

Alors * 


Ï 161 (x) Ga (x) dr = | eu (2) ga (x) dx. (3.21) 
0 0 


D'autres transformations intégrales dans lesquelles l’intégration 
s'effectue sur l’ordre des fonctions cylindriques sont considérées 
dans [1171]. 
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2. Trouver la solution de l’équation de la chaleur 
OT 1 OT 1 OT ) (3 29) 


ôT om _ 
on Ne 6 Sr+ r  ôr r?  04* 


qui vérifie les conditions 
T:=vdeist Jlorsque @— 0, 
ba | (3.23) 
T =0 lorsque =. 


Ce problème se pose lorsqu’on étudie un champ de températures 
périodique dans un domaine cunéiforme [184]. On cherchera la solu- 
tion sous la forme du produit 


T=deiviu(r, p). 
Pour la détermination de uw (r, @), on a l'équation 
Vu + ku = 0, kB= it. 
Moyennant le :hangement de variables 
z = ikr, ur, p) = V (z, ®), 


on obtient 
8®Y 1 0Y 1 02 
THE dz +7 —V =0, 


V =-1 lorsque ®—=0, 
V =0 lorsque = «. 


En appliquant la transformation de Kontorovitch-Lébédev (3.17), 
(3.18), il vient 


Va, q=e-ssme+ | A(s)sh (ps) Ki, (2) ds. 
0 
Les conditions aux limites entraînent 


\ A(s)sh(as) K;,(z) ds = —e-zsinz, (3.24) 
0 
Pour tous les B réels, on a 


| cos (Bs) Ke (2) ds = perche. (3.25) 
0 
La formule asymptotique pour de grands s 


LCI ES TA 
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entraîne aussitôt que l’intégrale de (3.25) converge pour les valeurs 


imaginaires de B — iu, pourvu que [u|+]|7| <<. En vertu 
de (3.24) et (3.25), il vient 


; 
A0 nd 2 s 
Donc 
V (z, ph=e-:sino + | ne ch (5-a) SK ;,(z) ds (S<a<T). 
0 


Bibliographie du chapitre III 
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CHAPITRE IV 


AUTRES TRANSFORMATIONS INTÉGRALES 


$ 1. Transformation de Mehler-Fock 


1. La transformation intégrale de Mehler-Fock est définie par 


© 


F(r)= | f(@)P. 


1 


1,04 720 (4.1) 


où P, (x) est une fonction sphérique de Legendre de première espèce. 
Soit f (x) une fonction réelle telle que f (x) P 1 (x) € L M, +ool. 


Alors l’intégrale (4.1) prise au sens de Lebesgue est une fonction 
réelle de t, définie pour toutes les valeurs de t > 0 (on peut égale- 
ment considérer les t << 0 avec F (—Tt) = F (x)). On calcule formel- 
lement les transformées de Mehler-Fock de fonctions quelconques 
à l’aide des représentations intégrales de fonctions de Legendre, par 
exemple 

@œ 


2 COS Ts 
de Er SPGedoe (a >0) 


(intégrale de Mehler), 


2 C COS Ts 
P = — h —————— , 
Ft. (cha)=<+c “| Fa (>0) 


et en modifiant ensuite l’ordre d'intégration. Le développement 
d’une fonction arbitraire en intégrale de Fourier 


f(z)= th mp is (z) dx f() re () dé (4.2) 


permet de trouver la formule d’inversion de la transformation (4.1). 
Formulons deux théorèmes [73]. 
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Théorème 1. Si une fonction 1 (x), définie sur l'intervalle 
[1, +oo! est telle que @ (t) =2sh _ Ÿ (ch t) possède une dérivée pre- 
mière absolument intégrable sur l'intervalle infini [0, +ool, et une 
dérivée seconde absolument intégrable sur tout intervalle fini, et si 
p (0) = 0, p (+o) = 0, alors 1 (x) peut être représentée par l'inté- 
grale 


v@= TP, (id (AKz<+o), (4.3) 
(A 
f(u)=uthur | P_1,,, 7%) dx. (4.4) 
d 2 


Théorème 2. Si une fonction f (x), absolument intégrable 
sur l'intervalle infini (0, Hoo[, possède une dérivée absolument inté- 
grable sur tout intervalle fini, et si f (0) = O0, alors f (u) peut être 
représentée par l'intégrale (4.4), où 1 (x) est de la forme (4.3). 

Ces théorèmes peuvent être démontrés moyennant d’autres 
hypothèses [134]. 

Envisageons maintenant la transformation généralisée de Mehler- 
Fock: 


va)= (Pl, (@ftu)du (A&z<+o), (45) 
0 2 


f(=+ sh (ru) (5—k+in) r (5—x—in) X 


k 
d hs 
En posant À — 0 dans les formules (4.5), (4.6), on obtient (4.3), 
(4.4). En se servant des formules 
i 


Pi. (cha)=(2xsha) ? x 
-g+iH 


(z)b(xz) dx. (4.6) 


- M meurs [+4 —kit+iu; Lex 


X csch a | + eian TE ST X 


XP 5 +84, ER; 1 — in ; — fe-scscha]}, 
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on déduit 
i 
—_ 1 
Pr. (cha)= (+ : ) (sha) ? cos (æu), 
u l (4.7) 
1 1 
?  (cha)= (+) © (ha) 7 sin (œu) | 
SRE n.) 
Si dans (4.5), (4.6) on pose À — … z = ch & et l’on tient compte 


de (4.7), on obtient : 
{ 
2 


1 : æ 
(sha) 24 (cha)= | (=) ) fs (1) cos (œu) du, (4.8) 
0 
1 
{ (u) = ) ) | (cha) (sh ni cos (uæ) da, (4.9} 
et de façon analogue pour k = — + 


er r 


Cho)? Ÿ (ch a) = (+) ( )_ sin (au) du, (4.10) 
0 


1 


oO 


Le (2) | + Cho ho zsntra)de. (4.11) 


2. La transformation intégrale de Mehler-Fock est utilisée pour 
résoudre certains problèmes de la théorie du potentiel, de la théorie 
de la chaleur, certaines équations intégrales linéaires ainsi que des 


problèmes de physique mathématique [134], [136]. Considérons à 
titre d'exemple l'équation intégrale 


FO %, (4.12) 


PR)=8() +1) = 


em 8 


« 


ou 
1<Lrz<+o, —00 << Àr < 1. 


Supposons qu’existent les transformées de Mehler-Fock des fonc- 
tions (x) et g (x) 


D(D= | pr) P_1,,,(5) de, 
1 


G(x) = EP 1 (0 dr 
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Comme 
= P 1 _. (s) 
ch nt ion 
LT) — —— ds, 
ni ) st | Z+s 


après multiplication de (4.12) par P_ te (x) et intégration sur x 


de {1 à —+o, on obtient 
D(r) =G(r) + D (1), 


d'où il vient pour —00 << Àn << 1 


D (=. 
7 chat 
La formule d’inversion donne 
p (x) = | tTth ar P 1 (2) ar. 
0 7 chat É 


On peut donner une justification rigoureuse en supposant par exem- 
ple que g (x) est une fonction continue à variation bornée sur tout 
intervalle fini J1, b[, et de plus que 

1) 8) P_ 1 @ELN, oo, 


2) G() TEL 10, oo! [134]. 


$ 2. Transformation de Hilbert 


Considérons la formule intégrale de Fourier 
1 (2)= Î {a (£) cos xt + b (t) sin xt] dt. 


: U 
Ici 


CO 


a(t) = ( Î (u) cos ut du, b( = | f (u) sin ut du. 
Posons | | 
g(x)= \ [b (t) cos zt— a (t) sin zt] dt — 
0 


=+ [ dt ( sin(u—zx)t-f(u) du. 
1 
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L'intégrale de second membre de la dernière égalité est appelée 
intégrale conjuguée de l'intégrale de Fourier et se déduit formelle- 


ment de la formule de Fourier par changement de a en b et b en —a. 
On aura toujours formellement 


l © 
ACENCE | dt { sin(u—zx)t-f(u) du = 
0 — 20 


se | En AN 2 f (u) du — 


lo TT EU—T 
_ 1(1— 
= lim — | TR jf (240) —f (x —t) dt. 
[—00 0 
D'où 
e()=+ | JEHI IE D ge, (4.13) 
0 
De façon analogue, 
f(z)— + [ sEHD EE 7, (4.14) 


0 
Les formules (4.13), (4.14) sont équivalentes aux formules 


[se œo 


_4 ADR RE f(x-+t)—f(z—t) 
g(a)= 2 V.P. 1 dt=— lim les ms dt, (415) 

_4 Feu np tou (etr+n—e(z—0) 
f@)= Ver. | dt im | de 


—œ 


(4.16) 


où V.P. signifie que l'intégrale est prise dans sa valeur principale 
au sens de Cauchy. Les deux dernières formules représentent un 
couple de transformations de Hilbert. Dans [238], [254] on trouvera 
des démonstrations rigoureuses de nombreux théorèmes se rattachant 
aux transformations de Hilbert. 


$ 3. Transformation de Laguerre 


La transformation intégrale 
T {f(t)}= | etLh(t)f(tdt=f"(n) (n=0,1,92,...), (417) 
0 


où L, (t) sont les polynômes de Laguerre d'ordre n, s’appelle transfor- 
mation de Laguerre. Celle-ci est tout indiquée pour la résolution de 
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l'équation différentielle de Laguerre 
Lz + nxz = 0, 


où 
Lit) = tx" (D) + (1 —t)zx (L). (4.18) 


L'application de la transformation de Laguerre ramène l'opération 
de différentiation £zx à une opération algébrique moyennant la for- 
mule 


T {£{xz (t)]} = —nr* (n) (n = 0, 1, 2, ...). 


Dans [159] on étudie la convolution pour les transformations 
de Laguerre 


oO T 
c(t) + | e”Ta(t) | eVircos® cos (V'itsin p) x 
(E Û 


Xb(t-+t—92 Vtt cos p) do dr 


et on construit l’appareil du calcul opérationnel pour les opérateurs 
de Laguerre £z. 


Bibliographie du chapitre IV 
[134], [238], [74], [38], [85], [159], [2541]. 
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CHAPITRE V 


CALCUL OPÉRATIONNEL 


$ 1. Notions fondamentales et propositions 


Soit’A1 l’ensemble de toutes les fonctions absolument continues 
ct, plus généralement, des fonctions à valeurs complexes, définies 
sur la section 0 < t < +o, L l’ensemble de toutes les fonctions 
définies sur la section 0 < it < + et intégrables-Lebesgue sur tout 
intervalle fini JO, Alf. 


Si une fonction F (t) € M, pour presque tous les £ existe la déri- 
[4 


vée F'(t) =f()EL et F(t) = F (0) + | f (u) du. Réciproque- 


ment, si g()€L, la fonction G (t) —  gtu) du appartient à l'en- 


0 
semble M et presque partout G’ (t) = g(t). Si f(HEL, g(t)E L, 
pour presque tous les t existe l'intégrale 


{ 
ffU—E 8 & 
0 


qui appartient à l’ensemble Z. L’intégrale 
t 


ROC HCE: 


0 


où G(t)E M, f(t)EL, appartient à l’ensemble A. Si 


t 
H(= | F(—-56E 4. 
0 
où F (EM, G() E M, alors existe H” (t) E AM. 


M est un ensemble linéaire et l’addition et la multiplication par 
un nombre y sont définies de façon naturelle. On appelle produit 
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des fonctions F () E M et G(t) € M l'expression 
t 
F(t)xG(9 = | F(t—E)G(E) dé. (5.1) 
D 
Avec une telle définition du produit point n’est besoin de distinguer 
les nombres des fonctions constantes. Voici les principales propriétés 
du produit : 

10. Si FD) EM et G(D EM, alors F() KXG(DEM. 

2°. Le produit est commutatif, i.e. F (t) x G(t) = G(t) x F (t). 

3°. Le produit est associatif, i.e. 

FOXGG) x HG) =F() % (G() x H (b). 
4°. Le produit est distributif par rapport à l'addition, i.e. 
FD xX(G(t+H()=FDxG(D+F(E) :x H (+). 

5°. Si pour tous les t € [0, +oo[ et F (t) 0 F(t) x G(t) = 0, 
alors G (t) — 0 [1651]. 

L'ensemble M forme un anneau commutatif relativement à l’ad- 
dition et la multiplication au sens (5.1). La propriété 5° traduit le 
fait que cet anneau ne possède pas de diviseurs de zéro. Tout anneau 
commutatif ne possédant pas de diviseurs de zéro peut être élargi 
à un corps de quotients. Appelons couple l’expression (F (t), G (t)), 


où F()EM,G(DEM et G(t) 0. Deux couples (F (t), G (t)) et 
(F; (t), G. (t)) sont équivalents si 


FD KG (1) = F;(t) x G (+). (5.2) 
Dans ce cas (et seulement dans ce cas) on écrira (F, G) — (F,, Gi). 
L'ensemble des couples (F, G) est divisé en classes d'équivalence. 
Comme tous les éléments équivalents à un élément d’une classe 
appartiennent à cette même classe, celle-ci est définie par un de ses 


éléments que l’on appelle représentant de la classe. Désignons par _. 


la classe qui contient le couple (F, G). En vertu de cette définition, 
. = ! si, et seulement si, F x G;, = F, x G. Définissons la somme 
{ 


et le produit des symboles . en postulant que 


Pipe Rte Ro) 

G Gi GK Gi L 

Fm Exe (5.3) 
G ” Gi G x. G: 

(G x G A0). cd 


Les seconds membres des dernières égalités ne dépendent pas du 


choix des représentants des classes F et mi. L'ensemble de tous les 


d* 
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symboles . forme un corps commutatif. Désignons-le par 9h. 


Appelons opérateurs les éléments du corps IR. Les formules (5.3) 
définissent la somme et le produit des opérateurs. Les opérateurs 


seront souvent désignés par une seule lettre, par exemple, a — _. : 
b — F., etc. Enfin, pour simplifier l'écriture, là où aucune confu- 
sion n’est à craindre, on notera ab le produit des opérateurs a et b. 
Donc, si a = _ et b = À sont des opérateurs, ab est leur produit 


H 
et de plus ab = a %x b = TR . Le symbole _ représente la divi- 


sion dans JX. Cette dernière se distingue singulièrement de la divi- 
sion ordinaire. 

Voici les principales propriétés de l’addition et de la multipli- 
cation dans le corps Y} : 

1°. ab = ba. 

20. (ab) c = a (bc). 

30. a(b + c) = ab + ac. 

On peut identifier l’ensemble des opérateurs de la forme “2 
à l’ensemble M ; en effet, à chaque opérateur 72 on peut associer 
une fonction F (t) € M. On dit alors que 9 est une extension de 


l’anneau M. Au lieu de . on écrira 1: cet opérateur s’appelle éche- 


lon unité de Heaviside. Au lieu de _. on écrira (0. 


Les constantes appartiennent au corps WW. La division des cons- 
tantes dans I se confond avec le quotient ordinaire de constantes. 


L'équation 
F(t) xz(t) = G (1) (5.4) 


n’admet pas toujours une solution dans l’anneau M. Par exemple, 
si F(t) —tet G(t) = 1, alors l'équation (5.4) prend la forme 
14 
d 
7 | GHz ® dt =1, 
0 
ou 


Î x (E) dë=1, (5.5) 


0 


et l’ensemble M ne contient pas de fonction zx (f) vérifiant l’équa- 


tion (5.5). Dans le corps M, toute équation ax = b (a 0) possède 


b 


toujours la solution zx = a-!b — és 
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Le corps M contient l’ensemble ZL. Supposons que f(t)ELZ, 

t 
alors F (t) = | Îf (u) du € M et F (0) = 0. Associons à la fonction 


0 
Îf (t) l’opérateur 2. A deux fonctions distinctes f (t) et g (t) de 


l'ensemble Z PEN les opérateurs distincts “0 et SU) ; 


puisque dans le cas contraire t x F = t KkGout x (F—6G) = 0, 
d’où F (t) = G (t). A tout opérateur 22, H (0) = 0 est associée 


une fonction À (t) = H” (t) € L. Donc, Don avons établi une corres- 
pondance biunivoque entre l’ensemble des opérateurs de la forme 
0 et l’ensemble ZL. A la somme des fonctions f (t) + g (t) est 


7Q Tea 70 


fonction f (t) par À le produit de ie 20 par À. Les pro- 
priétés énumérées permettent d’assimiler les opérateurs de la forme 
0 , F (0) = 0 aux fonctions de l’ensemble ZL et de ne faire aucune 


associée la somme des opérateurs , au produit de la 


EN entre la fonction f (t) € L et l’opérateur 70 . On écrira 
donc 
20 - r' (5 =f(. (5.6) 
Il est naturel SR fonctions les opérateurs du corps M 
qui se ramènent à la forme _ , F (0) = 0, et na gd . 
par 
Les opérateurs n’admettent pas tous une Salisations par HAN 
0 , F(0) = 0. Une 
somme de fonctions est toujours une fonction. Un exemple simple 
nous montre que le produit de fonctions n’est pas toujours une fonc- 
tion, i.e. le produit des opérateurs f (t) — Let g(t) = L(F (0) — 


= 0, G (0) — 0) ne peut pas toujours être ramené à la forme . 


opération qui consiste à le ramener à la forme , F (0) = 


l'opérateur 7 ne peut être ramené à la forme —— 


H (0) = 0. En effet, supposons que f (ft) — es , £(t) = TA 
ie. Ft) = Vi, G(t) = Vt. Alors 

{ 
+ | VS 

0 


F(DXKG(t) _ 
à t KE ’ 


FE) Xk 8} = 
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d'où 
- ELLE 
OX 0-7. 


Donc, l’ensemble des opérateurs de la forme 8 , F (0) = 0 n’est 


pas un anneau. Cependant on a le 

Théorème 1. Si l'une des fonctions f (t) ou g (t) est absolu- 
ment continue, i.e. appartient à l'ensemble M, alors le produit de ces 
fonctions est de nouveau une fonction ayant pour expression 


t 


d L 
f(x (= + | f(E—E) 8 (6) di. (5-7) 
0 
L'opérateur 72 peut être considéré comme le produit de l’opé- 


rateur Z par la fonction 20 — F (t). L'opérateur L est fonda- 


mental dans le calcul opérationnel. On le désigne spécialement par 
s | Ha e 
— La formule (5.6) s'écrit 


PF (= F'(t. (5.8) 


Ici F (HE M et F (0) = 0. Si F (t) est une fonction quelconque 
appartenant à Af, de (5.8) il vient 


PF (ti) = F°(t) + pF (0). (5.9) 


Donc, lorsque F (t) € M et F (0) = 0, le produit de la fonction 
F (t) par l’opérateur p désigne tout simplement sa dérivation. Le 
produit pa a un sens pour tout opérateur a € Y; en général pa 
est un opérateur. Pour que pa soit une fonction, il est nécessaire et 
suffisant que a — F(t)E M et F (0) = 0. L'opérateur p est appelé 
opérateur de dérivation. Si F"(t) € M, de (5.9) il s'ensuit que 


PF (&))=p'F(t) =F"(t) +pF"(0) + p°F (0). (5.10) 


D'une façon générale, lorsque F (t) possède une dérivée d'ordre n 


appartenant à l’ensemble M, l'application répétée de la formule (5.9) 
donne 


P(P'F(t) =p°F()=F®()+ pET (0) + 
+ PFeT2 (0) +... + p"F (0), (5.11) 
où p" est le produit de # opérateurs p X p X ... X p dans le corps SN. 


Si dans (5.9), on pose F (t) = e°*, il vient pe“ — ae! + p ou 
(p — a)e* = p, d'où 


P__ _ at 
— =": (5.12) 
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De façon analogue, en posant F (f) = sin wt, F (t) = cos wt dans 
(5.10), on obtient 


7 sin ot, = =: COS f. (5.13) 
L'égalité 
p (te°!) = ateñ + eñf 
entraîne 
(p— a) tet! — eñt 
ou (cf. (5.12)) 
b= PP; 
(p— a) te” — p—a ? 
d’où 
P __ #oût 
Ga 
Dans le cas général, on a 
P ”_ tn-1e-at 
(p— a) 0 T(n) ‘ (5.14) 


L'opérateur p”}, inverse de p, est de toute évidence égal à &. 
La fonction F (t) = t appartient à M. Donc de (5.1) il s'ensuit que 
quelle que soit f(t)E L 


t 
pif (t) = | 1 (dE. (5.15) 
0 


Donc, le produit de p-! par f (t) n’est autre que l’intégration. L’opé- 
rateur p”! s'appelle opérateur d'intégration. Soit p" le produit des 
n opérateurs p-! Xp? x...>xp'!; de (5.15) il vient 
PF =—. (5.16) 
La formule (5.1) donne pour f(t)E€ L 
Î 
n 1 n— œæ 
PL = pr | D" F (0 dE. (5.17) 
0 


$ 2. Opérateurs rationnels 


L'objet le plus important du calcul opérationnel est l’étude des 
opérateurs de la forme Z (p), où Z (4) est une fonction de la variable À. 


Exemple: les opérateurs _— rw , p" étudiés au paragraphe 


précédent. Dans le cas simple où Z (à) est une fonction rationnelle, 
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1.0. 
__ PO) 
Z (À) Q (À) ? 
où P (A) et Q (À) sont les polynômes 
P(à)= D an et Q(A) = D Br”, 
R—0 Rke=( 
il est naturel de définir l'opérateur Z (p) à l’aide de l'égalité 
__ P(p) 
Z(P)= 5: (5.18) 
Les opérateurs P (p) et Q (p) appartiennent au corps 9, donc leur 
quotient 0 appartiendra également au corps %K. L'opérateur 


Q (p) 
Z (p) est appelé opérateur rationnel. Il s'agit maintenant de savoir 


pour quelles fonctions rationnelles Z (à), les opérateurs Z (p) se ra- 
mènent à des fonctions. Pour cela il est nécessaire et suffisant que le 
degré du polynôme LP (À) soit inférieur ou égal à celui du polynôme 
Q (à). Donc, si n < m, il existe une fonction ® (t) telle que 


Z (p) = œ (t); (5.19) 


en outre il est aisé d’établir que œ (t) appartient à M. Si la valeur 
de l'opérateur Z (p) est connue sur l'échelon unité, la valeur 
Z (p)f(t), où f (t)E L, est donnée par la formule 


L 
Z(P1O=+ | vU—DT EE, pOEM. (5.20) 
0 


La fonction œ (ft) de l'équation (5.19) peut être déterminée par 
une décomposition de À-!Z (À) en fractions simples. On pourrait se 
servir des théorèmes de Heaviside de décomposition (cf. chapitre II, 
$ 3). Si, par exemple, les racines du dénominateur @ (À) sont sim- 
ples et Q (0) 0, alors 

…: P (0) P (à») a t 9 
Z(P)=T us FRQ Ga) (521) 
Ici Au, Àos + « + Âm SOnt les racines du polynôme Q (à). Si Q (0) = 0, 
cela veut dire que Q (p) = pQ, (p), où @, (0) 0. En se servant de 
P (p) 


(5.21), on trouve GW — % (4), et ensuite 


Z (p)= p"! X pi(t)= | P1 (u) du. 
(4) 


Bien plus complexe est le cas où les racines sont multiples. Si, 
par exemple, À = À, est une racine multiple, en décomposant 
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À"1Z (À) en fractions simples, il apparaîtra des fractions Pt | 
Donc, œ (t) contiendra des termes À,t lei! et dans le cas général, 
la solution de l'équation (5.19) sera de la forme 


p(t)= D Autrternt, (5.22) 
kr 


Réciproquement, pour toute fonction œ (t) de la forme (5.22), il 
existe un opérateur rationnel Z (p) tel que Z (p) =  (t). Les opé- 
rateurs rationnels ne parcourent pas le corps + tout entier. Ils 
forment un sous-corps de 9}, qui s'obtient par extension de l’anneau 
des fonctions (5.22). Tout élément de ce sous-corps peut être repré- 


senté par un opérateur rationnel Z (p) = 52 . Sir <m, Z (p)est 


une fonction, sinon c’est un opérateur. Le résultat précédent peut 
être généralisé à une classe plus large d'opérateurs du corps 9h. 
Cette généralisation s’opère le plus facilement à l’aide de la trans- 
formation de Laplace. 


$S 3. Opérateurs transformables-Laplace 


Soit S l’ensemble de toutes les fonctions f (t) € L telles que 
l'intégrale de Laplace 


f@= | r@etdt 
0 


soit absolument convergente, et S* l’ensemble de leurs transfor- 
mées de Laplace. Si l’opérateur a € 9 admet un représentant 
(F (t), G (t)) tel que F (D ES et G(H)ES, cet opérateur est dit 
transformable-Laplace. La fonction 


__ F*(2) 


a (z) MTEICÉ (5.23) 
où 
_ F(:) 
a — TAUX , 
F*(2)= Î F(the-ztdt, G*(2)= [ G(t)e-:t dt, 
0 ) 


est la transformée de Laplace de l'opérateur a = za . Ï est aisé 


de vérifier que la définition de & (z) ne dépend pas du choix du 
représentant (F, G) et donc a (z) est défini de façon unique par 
l'opérateur a. Cette transformation est noté symboliquement 


a a(z}. (5.24) 
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L'ensemble de tous les opérateurs transformables-Laplace du 
corps 9h est noté J%, et l'ensemble de leurs transformés %. 

Théorème 2. La transformation (5.24) établit entre les 
ensembles St et VÈ une. correspondance biunivoque qui associe à la 
somme d'opérateurs a + b la somme des fonctions a (2) + b (2), et au 
produit a X b le produit ordinaire a (z)b(z). Le zéro et l'unité de 
l'ensemble KR se transforment en le zéro et l'unité de l’ensemble R. 

Ce théorème établit l’isomorphisme des corps % et JE. Ceci 


Ld # # e e 4 ee # «+ e 
étant, l’opérateur de dérivation p — — est associé à la fonction 


& (z) = 2, i.e. p ==. Donc, pour tout opérateur rationnel, on a 


__P(p) =. P(:) e 
Z (p) — TE) — 0G = Z (2). (5.25) 
Si l’opérateur est une fonction f (£) € S, alors 
f(t) =zf* (2) = f (:). (5.26) 


Les ensembles J et Ÿt étant isomorphes l’un à l’autre, la distinction 
qui existe entre eux n’est pas essentielle dans la plupart des cas, 
ce qui permet de ne pas différencier l'opérateur p du nombre com- 
plexe z. Aussi, souvent remplacerons-nous z par p. Dans nombre de 
cas la notation de l'opérateur de dérivation et d’un nombre comple- 
xe par une même lettre simplifie l'exposé. Donc, p désignera l’opéra- 
teur de dérivation dans le corps Y et un nombre complexe dans le 
corps{). Au lieu de (5. 26) on écrira 


{ (t) = j E). (5.27) 

Dans les notations adoptées f* (p) représente la transformée de La- 

place de la fonction f (t) € S, f (p) la transformée de Laplace de l'opé- 
t 


rateur f (t) = 20, où F (1) — | f (u) du. Entre f (p) et f* (p) on 
0 


a la relation suivante f (p) — pf* (p). Dans le cas général, sia EX 
et a = AU) alors 


1 n\  F*(P) 
L'’isomorphisme (5.28) ramène l'étude de la structure du corps 
des opérateurs 9 à celle du corps ÿ. Ainsi donc, dans le corps des 


opérateurs J on peut définir un sous-corps 9 isomorphe à un corps À 
dont les éléments sont des fonctions de la variable complexe p d’un 
type défini. Tout opérateur du corps peut être considéré comme une 
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fonction de l’opérateur p, i.e. a = à (p). Cette dernière propriété 
peut être étendue au corps Ÿ! tout entier moyennant la transforma- 
tion de Laplace généralisée [46], [491]. 


$ 4. Sur la réalisation des opérateurs 
transformables-Laplace 


Le corps des opérateurs %, isomorphe au corps à, forme dans 
une classe d'opérateurs qui est vaste et importante pour les applica- 
tions. Toute fonction mesurable f ({) qui ne croît pas plus vite que 
el, i.e. vérifie pour tous les £ assez grands la condition 


|f(t)1 <Kef, (5.29) 


où À et q sont des constantes, appartient à l’ensemble S. En général, 
pour qu’une fonction f (t) € L appartienne à S, il est nécessaire et 
suffisant que pour une constante q 


lim e% | f (u) du = 0. 


{—20 0 


Donc, le corps des opérateurs 9 est une extension de l’anneau des 
fonctions vérifiant la condition (5.29), où les constantes Æ et q ne 
sont pas les mêmes pour toutes les fonctions. Si l'opérateur a est 


susceptible d’être ramené à la forme —— 0 


f (t) appartienne à l’ensemble $, alors 


a = +0 = pf° (p), (5.30) 


— f (ét) et que la fonction 


où 
f(p)= | 1@)ert at. (5.34) 
0 


La réalisation de l'opérateur est primordiale lorsqu’on résout des 
problèmes par les méthodes du calcul opérationnel. Il est indispensa- 
ble de disposer de critères qui permettraient de conclure si oui ou non 
l'opérateur considéré peut être ramené à une fonction. Les théorèmes 
et lemmes du $ 1, chap. IT, fournissent dans nombre de cas une répon- 
se à cette question. 

Ici et dans la suite, on rencontrera des fonctions multivoques 


V P; In p, arctg p, etc. Dans tous les cas, si le contraire n’est pas 
spécifié, on envisagera toujours la détermination de la fonction pour 


laquelle }/1 = 1, In1 = 0, arctg 1 — , avec arg 1 = 0, i.e. la 
détermination principale de la fonction. 
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Par exemple, l'opérateur V/p se ramène à une fonction, puisque 


4 vérifie le théorème 5 du chapitre II. Il en est de même des opé- 


rateurse-*V?, ar (0 <'u< 1) en vertu du théorème 6 du chapi- 
chyp 
e e # e 1 1 
tre II, ainsi que des opérateurs sin É In (1 + =) 7 o (=): 


1 
p (eP — 1), en vertu du théorème 3 du chapitre II, etc. 

Si une fonction & (p) est régulière au voisinage d’un point situé 
à l'infini, l'opérateur à (p) est dit régulier. 
Théorème 3. Si a (p)est un opérateur régulier et f(t)E L 
alors 


a(pf(t)= S a —f(), (5.32) 
R=0 P} 
où 
a(p)= 2x |P| > ep. 


La série (5.32) est uniformément convergente sur tout intervalle 
t 
= “= 


Par exemple, l’opérateur v = nc est régulier, donc 
1 
P_ M) 2 4 _ 1, 13 » 
VAR an (147) Tim idam ot 
1 (A2, 1.3 (214 
=1—- e +57 Z 1 — ... = Jo(Àt). 
À 


De façon' analogue, en développant e P suivant les puissances de 
n , on obtient 


Lorsqu'un opérateur se ramène à une fonction, les valeurs de 
cette dernière sont obtenues à partir de (5.31) moyennant l'intégrale 
(cf. chapitre II, $ 1) 

; Y+ico 
fO= | f'(pe" ap. (5.33) 


Ÿ—100 
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Pour obtenir une expression commode au calcul de la fonction f (ft) 
comme on l’a indiqué au chapitre II, il faut souvent changer de façon 
adéquate le chemin d’intégration dans les formules (5.33). Ceci 
étant, dans la plupart des cas, il s'avère possible d'utiliser les lemmes 
de Jordan (cf. chapitre II). On trouvera des exemples de réalisation 
des opérateurs dans [246], [50], [35], [146], [64]. 


$ 5. Transformation généralisée de Laplace 


Les propriétés standard de la transformation de Laplace entraïnent 
que S* (cf. $ 3) est un anneau par rapport à l’addition et à la 
multiplication ordinaires. L'ensemble S* ne contient pas de constan- 
tes. Cependant le produit Àf* (p) € S quel que soit À. Dans la suite, 
nous étudierons l’anneau S* sur le corps des nombres complexes. 
Soit J, l’ensemble de toutes les fonctions de S* représentables sous 
la forme e-r%g* (p), où g* (p) E S* et w > 0. Le produit e-P®g* (p) 
appartient à S* pour toute fonction g* (p) € S*. Il est évident que J, 

+ 
est un idéal *) de S*. Formons l'anneau S% — =. Ilest standard 
w 
que les éléments de S% sont des classes mixtes par rapport à J,. 
Plus exactement, l’ensemble S* se décompose en ensembles f,, tels 
que deux éléments f* (p) et f? (p) appartiennent à un même f, si, 
et seulement si, f* (p) — f? (p) appartient à J,. Dans ce cas les 
éléments f* (p) et f?(p) sont dits comparables sur l'idéal J,. On 


écrit souvent 
fi (P)= f2 (P) (Jo): 


Soient f, et g, des éléments de S$%. Il est standard que la somme 
fo + & et le produit f,g, des classes f,, et g, sont par définition les 
classes qui contiennent les éléments f* (p) + g* (p) et f* (p) g* (p), 
où f* (p) et g* (p) sont des représentants des classes jf, et g,, 1.e. 
ff EH)Ef, et g* (p) € gw- 

Remarque 1. Supposons que © > p; il est évident que 
Js € J,. Donc, si les classes f, et f, contiennent au moins un élé- 
ment commun f* (p), alors f, € f, (f, contient l’ensemble f,). 

Composons la somme directe des ensembles S%, w parcourant 
l’ensemble ordonné de tous les réels positifs, 0 << &@ << + oo. On 
sait que les éléments de la somme directe sont les systèmes {f,}, 
0 << © << + , où f, est un élément de Sé. Si {g,} est un autre 
élément de la somme directe, alors par définition 


{fo} + {go} = {fo + Lu}: 
{fo} (£w} = (fu£u}: 
À {go} = {Lo} - 


*) On appelle idéal sur un anneau S* un ensemble d'éléments jouissant 
des propriétés suivantes: 1) si a € J,,et bEJ,, alorsa+beJ,;2)sia€J,, 
b E J,,, alors le produit ab € J,. 
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Parmi] les éléments {f,} de la somme directe des ensembles S®, 
considérons ceux pour lesquels est réalisée la condition /, = f, tou- 
tes les fois que w << p. Ici f, et f, sont considérés comme des ensem- 
bles de S*. Désignons par L* l’ensemble de tous les éléments de la 
somme directe et par {f,} ces éléments eux-mêmes. On montre que 
l’ensemble L* est un anneau dont l’élément nul est {J,,}. 

Remarque 2. Supposons que dans l'anneau S* est donné 
un opérateur linéaire 7 dont les valeurs appartiennent à S*. Suppo- 
sons que l’idéal J,, est invariant par l’opérateur 7, i.e. si f* (p)E J,,, 
alors 7f* (p) € Ja. Dans ce cas l'opérateur T induit sur l’anneau 
S% un opérateur linéaire que nous désignerons également par T. 
Si fo € So, alors Tf, désigne une classe dont un représentant est 
l'élément 7f* (p), où f* (p) € fw- 

Supposons que {f*} appartient à l'anneau ZL*. Définissons 
T {f»} en supposant que T {fs} — {Tfs}. Considérons l'intégrale 


o 


fs (p) = | f(the-rtdt, f(t)EL. (5.34) 
0 
Soit f, la classe mixte de l’ensemble S$ — > qui contient la fonc- 
tion fà (p). Comme pour p > & 
e pu 
f(p)—f5(p)= | fe-rtdte-e | f(t+u)e-rtdt, (5.85) 


0 


o 
alors f5 (p) — fo (p) E Jus. Donc fi (p)E fu et fo = fp- De là il 
s'ensuit que {f,} € L*. L'élément {f,} s'appelle transformée généra- 
lisée de Laplace de la fonction f(t). La correspondance entre les en- 
sembles L et L*, engendrée par la transformation généralisée de 
Laplace sera désignée symboliquement par : 


Î (1) = Üo}- (5.36) 


Cette notation se justifie par le fait que la transformation générali- 
sée de Laplace se confond avec la transformation ordinaire si celle- 
ci existe. 

Propriétés de la transformation généralisée de Laplace 

10. Si ft) ={f} et g(t) ={g}, alors A f(t) + ng(t) = 
= À {fs} + u {go}, où À et u sont des constantes. 

20. La transformée généralisée de Laplace {f,} définit la fonction 
f (t) de façon unique à un ensemble de mesure nulle près. 

30. Théorème 4 (théorème généralisé de 
Borel). Soient f (t) = {fo} et g(t) = { go}, alors 


f(t—E) 8 (© Æ = Jo} {8u}- 


Qu) 
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40. Soient f (t) == {f.}- Pour que la fonction f (t) appartienne 
à l’ensemble S, il est nécessaire et suffisant que l’intersection de tous 
les ensembles f,, &o < © << + oc (f., est considéré comme un ensem- 
ble de S*), ne soit pas vide. 

Les autres propriétés connues de la transformation ordinaire de 
Laplace se transposent à la transformation généralisée de Laplace. 

5°. Théorème 5 (théorème de retard). Si 


f() = {fo} et 
lorsque O<i< rt, 
fe (t) = | f(t—T), lorsque O<r<t, 
alors fi (t) = {e-P'fu} où fit) = e-" {fu}. 


6°. Théorème 6 (théorème de translation). 
Sif(t) ={fu}, alors f (t) e* = {fs (p — a)}. 


fw (P— a) désigne ici la classe mixte par rapport à l'idéal J,, 
o 


dont un représentant est l'élément fé (p—a)=— [ est-r1} (t) dt où 
0 
fa (p) est un représentant de la classe f,. 


7°. Théorème 7 (théorème de similitude). Si f(t) = {fs}, alors 


. À 
f(at) = — {(+) } ; 
le (2) désigne ici la classe mixte par rapport à l'idéal J, dont 


un représentant est l'élément (+), où fo(P)E fo: 


S 6. Le corps 


Au $ 5, on a construit l’anneau L* des éléments {f,}, transformées 
de Laplace généralisées. L’anneau Z* ne possède pas de diviseurs 
de zéro. Ceci découle du théorème généralisé de Borel et du théorème 
de Titchmarch. Donc, l’anneau Z* peut être étendu à un corps de 
quotients. Procédons comme au $ 1 pour l'extension de l’anneau M 


et introduisons les symboles Del Où {£gu} = J. Par définition, on 


a Dei — fu dans le cas seulement où {f,}{ko} = {£w}{ho}. 


Ce qui veut dire que pour tous les représentants fà (p) € fu, £® (P) € 
€ Sos ko (P) E ko et hG(P)E hs, on aura 


fo(p) KG (P) = g6(P)hG(P) (Jo); >. (3.37) 


L'ensemble de tous les éléments de la forme qe constitue un 
us A] 
corps que nous noterons Y. Donc, à un isomorphisme près l’exten- 


96 CALCUL OPÉRATIONNEL [CH V 


sion de l’anneau L* à un corps de quotients coïncide avec.le corps It. 


Soit maintenant a = Fe un Opérateur quelconque du corps . 


Associons à l'opérateur a un élément du corps Mt égal à er, 


où {F,} et {G,} sont les transformées généralisées de Laplace des 


fonctions F(t) et G(t) respectivement. Appelons l'élément a = Del 


transformée généralisée de Laplace de l'opérateur a. Désignons cet- 
te correspondance par a — «a ou encore de façon plus détaillée 


F( __ {Fo} 
G (+) = {Go} ® (5.38) 


Il est immédiat de vérifier que (5.38) est une correspondance biuni- 
voque qui à une somme d'opérateurs de WW associe la somme des élé- 


ments correspondants de Y%, à un produit d’opérateurs de M le 

produit des éléments correspondants de M : à l'élément nul de M 

l'élément nul de YŸ*, à toute constante de Mt la même constante de MN, 

et enfin à l'opérateur = l'opérateur p. Si F(t) et G(t) admettent la 

transformation ordinaire de Laplace, i.e. les intégrales [F (t)e-?t dt 
0 


et j (t)e-P'dt sont absolument convergentes, alors la propriété 
0 


4° de la transformation généralisée de Laplace (cf. $ 5) entraîne 


FO _ FD > 
Guy — esp UP), 


« e e e LU = Ld 1 La e 
d’où il vient, en particulier, qu’à l’opérateur est associée la fonction 


a (p) = p. L'isomorphisme des corps 9% et JR définit la structure 

du corps des opérateurs de Mikusinski. Il s’avère que les propriétés 

des opérateurs du corps W? sont étroitement liées à celles des fonc- 
R 

tions complexes de la forme TS , où à (p) et £e (p) sont des fonc- 

8w\P 

tions entières, représentables par les intégrales 


o o 


fs(p)= | fe-rtdt et gë(p)= | gte-v'at. 


0 0 
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$S 7. Fonctions opérationnelles 


Les opérateurs dépendant d’un paramètre se rencontrent dans les 
applications du calcul opérationnel aux problèmes de physique ma- 
thématique. Dans ce paragraphe, on étudiera des opérateurs dépen- 
dant d’un paramètre réel. Si un opérateur a dépend d’un paramètre À, 
a LÀ L b, on écrira a = a (À) et on dira que a (À) est une fonction 
opérationnelle. Toute fonction opérationnelle est définie par son 
représentant (F (t), G (t)). Les fonctions F et G dépendent du parame- 
tre À, i.e. dans le cas général F —F (t; À) et G = G (t; À). De plus, 
la fonction G (t; À) n’est identiquement nulle pour aucune valeur du 
paramètre À. 

Voici des ni de fonctions opérationnelles : 


a (À) — 


re (—o<A<+ oo), 


a (A) = e-AVr= = 


2.5 
e  (—oo<À<+ oo), 


a == AVR _ 
0 pour {<À, 
= | Jo(V EX), pour 2<t CES Fe), 


0, pour {<À, 


a(A)=e"P=n(é; À) — U boue Le: (OS << + 00). 


La fonction opérationnelle a (À) est par définition continue sur l’in- 
tervalle a < À < b si existe un représentant (F, G) de la classe z 
tel que les deux fonctions F (t ; À) et G (t ; À) soient continues en les 
variables À et t dans le domaine a À << b, 0 t< + co. Par 
exemple, la fonction opérationnelle e-4? est continue sur [0, 1]. 
On remarquera que la continuité d’une fonction opérationnelle sur 
(a, b] ne signifie pas que tous ses représentants, i.e. les couples (F, G), 
sont constitués de fonctions continues. 


$ 8. Limite d’une suite d’opérateurs. 
Limite d’une fonction opérationnelle 
On dit qu’une suite d'opérateurs a, € M est convergente vers un 


opérateur a = si existent des représentants (F,, G,) tels que 


1F 


7—0394 
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2) les suites F, (t) et G, (t) convergent respectivement vers les 
limites F (t) et G (t) uniformément sur tout intervalle fini0<t<7T, 


lim Fa(=F(t), limG,(t)=G(t). 


On dit alors que l'opérateur a=À est la limite de la suite 


G 
d'opérateurs an = et l’on écrit 
lim a, = a. (5.39) 
On montre que la définition de la limite ne dépend pas du choix des 


représentants. 

Citons deux exemples *). La suite a, — cos nt est divergente 
au sens ordinaire de l’analyse classique. Au sens opérationnel, cette 
suite converge vers zéro. En effet, 
sin nt 

n lim sin nt —( 


n—+00 


et la convergence est uniforme. Donc, lim a, = 0. Il en est de même 


ñn — © 
de la suite 

… sin n{ 

A n 
An = 7 Sin ni — FE x 
F2 

On a de toute évidence 

4 

liman=———=p. 


Considérons maintenant la suite a, = ne". Il vient 
t 


— nent 
n 2 
d'où 
t 

lim net —p 
t- 

n—0 ns 
2 


Donc au sens opérationnel, la suite ne"‘ converge vers l'opérateur —p. 


*) On rappelle que dans ces exemples et dans les suivants, la multiplication 
et la division sont effectuées sur le corps M. 
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Propriétés fondamentales de la limite d'une suite. Si une suite a, 
est convergente vers une limite, chacune de ses suites partielles con- 
vergera vers la même limite. 

Si existent lim a, — a et lim b, = b, alors les suites a, + b,, 


n — © n — © 
An — On, Anbh possèdent également des limites telles que 


lim (a, +b,)=a—+b, lim a,b, = ab. 


> 00 nn 0 


En particulier, si lim a, = a et c est un opérateur quelconque, alors 


n — © 
lim ca, = ca. 
Ti 0 
Si limb, —b#0, alors 
T1 — 00 
se bh b 


En se servant de la limite d’une suite d’opérateurs, on peut introdu- 
ire la notion de limite d’une fonction opérationnelle. Plus exacte- 
ment. une fonction opérationnelle a (À) admet une limite au point 
à = À, si quelle que soit la suite À, convergente vers À, existe la 
limite lim a (4,) et cette limite ne dépend pas du choix de la suite 


Ân His Donc elle ne peut dépendre que du point À,. Dans ce cas 
on écrit 

lim a (À) = b. 

À—=ho 


Corollaire. Si la fonction opérationnelle a (À) est continue 
sur l'intervalle a < À < b, alors existera lim a (À) = a (À,) quel que 


À — À0 
soil À — Lo. 


$ 9. Dérivée continue d'une fonction opérationnelle. 
Intégrale d’une fonction opérationnelle 


Si une fonction opérationnelle a (4), a < À < b, admet un repré- 
sentant (F (t, À), G (t, À)) tel que: 
1) les fonctions F (ft, À) est G(t, À) soient dérivables par rapport 


ah, ab, et les dérivées D F(t, À) et —G,(t, à) 


appartiennent à l'ensemble W ; 

2) les fonctions F; (t, À) et G; (t, À) soient continues en t et À 
dans le domaine 0 <t< + oo, a < À < b, alors on dit que la fonc- 
tion a (à) est continüment dérivable sur l'intervalle a &< À < b. On 
appelle dérivée continue d’une fonction opérationnelle dans I la 
classe dont un représentant est le couple (F, XG—F x GirEXK QG). 


\7 
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SR Donc, 


La dérivée est désignée par a” (À) ou 


FG—FG; 


a'À — E 


(5.40) 

On montre que la définition de la dérivée ne dépend pas du choix 
du représentant de la classe. 

Propriétés de la dérivée: 

40. [a (à) +6 (2)]' = a' (à) + db" (À). 

20. [Ca(À)]j Ca’ (À), C est un opérateur constant. 

3. [a (À) b (A) = a’ (À) 8 (à) + a (A) D (À. 

40 ES a’ (À)b(2)—a (À) 200, 
° bQ) s b2 (À) 

59. Si a” (À) = 0 lorsque a À  b, alors a (4) = a (u) pour 
tous les a LA<u< b. Cette F propriété a été démontrée en 1962 
par L. Berg [121 sous certaines hypothèses. On n’en connaît pas 
encore la démonstration dans le cas général. 

Définition de l'intégrale d'une fonction opérationnelle. Si pour 
une fonction opérationnelle donnée a (à) il existe une fonction opéra- 
tionnelle À (à) telle que sur l’intervalle & < À < B À (À) possède 
une dérivée continue A” (À) et À’ (À) = a (À), alors l’opérateur À (à) 
est appelé intégrale indéfinie de la fonction opérationnelle a (à), ce 


qu’on note | a (À) dA. La propriété 5° entraîne que la fonction À (à) 


est définie à un opérateur constant près. On appelle intégrale définie 
de la fonction opérationnelle a (à), l'opérateur 

B 

| a (À) dà == A(B)— A (a). 


œ 


L'intégrale d’une fonction opérationnelle possède des propriétés 
Tor à celles de l'intégrale ordinaire. En effet, 


e. [aa 0. 
8 . 
20. | a) d? = | a (À) di. 
8 Y B 
30. | a (À) d\ = J a (À) an + | a (À) dÀ. 
œ Y 
B B 
49, 1e (à) dà =c Ï a(A)dX, où c est un opérateur ne depen- 


dant pas de À. 


B B 
5, | [a (A) b (A)] dà — [ a (À) dr + | b(à) dà. 


Le d œ œ 
6°. Si a(À) et b(A) possèdent des dérivées continues sur l’inter- 
B 


TT: 


valle a<A<b, alors Î a’ (à) b(A) dà = a(B) b(B) — a (&) b(æ)— 


B [+ 
| a (À) 6” (À) dà. 
Ë 


nfin, on peut introduire les intégrales impropres en posant 


© 


B 
au) di= lim | a (À) dà, 
[0 4 ee œ 


où la limite est comprise au sens opérationnel, i.e. pour toute suite 
de nombres $f, —> oc, existe la limite de la suite d'opérateurs 
Bn 


[a (À) dX et cette limite ne dépend pas du choix de la suite de nom- 


bres B.. 
$ 10. Fonctions en escalier 


Définition. On dit qu’une fonction f (t), 0 < t + o est 
en escalier si l’on peut partager la section ] 0, +  { en un nombre 


Fig. 4 Fig. 5 


fini ou dénombrable d'’intervalles disjoints à l’intérieur desquels a 
fonction f ({) reste constante. 
Exemple 


0, O<t<À, 
en = | _. (5.41) 


1, At. 
Le graphe de cette fonction est représenté sur la fig. 4. 
Formons l'opérateur e-P} — e-Ph, up >> À. Il est évident que 
0, pour 0<i<X, 
e-Ph— e-Ph — d pour AS{<y, 
0, pour u<t. 
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Cette fonction (fig. 5) est également en escalier. Soit 0 = À, << À << 
<< .-..<n LL... une suite monotonément croissante de 
nombres réels, lim À, — et us, WU, . . ., Uh, . . . une suite quel- 


n — 2% 
conque de nombres réels. Il est aisé d’établir que la série 


D une? = q (t) 
n- 0 


est toujours convergente. Sa somme est une fonction en escalier. En 
effet, quel que soit t fixe, on peut toujours exhiber un entier { tel 
que An tn. De (5.41) il résulte alors 


x 
P (t) F 2 Un 


lorsque Àxy <St<< À, et en particulier œ (t) = w lorsque 0 << 
<< t < À. Réciproquement, si est donnée une fonction en escalier 
œ (t) et que la fonction œ (t) soit égale à p, à l’intérieur des interval- 
les] Az, An+1 [CE —= 0, 1, 2, . . .), alors de (5.41) il s'ensuit 


p)= 2 pu (een?) (5-42) 


Ainsi l’ensemble des fonctions réelles en escalier coïncide avec 
l’ensemble de tous les opérateurs de la forme Duye-hxr, où 0 — 
R 


= À Lu... <Ln <L..., Àn — 00 et u, sont des réels. Dans 
les applications on rencontre très souvent des séries dans lesquelles 


les nombres À, A1, . . . forment une progression arithmétique À, = 
= kh (k = 0,1,2,...). Dans ce cas, 


p(t)= > ue *?. 
k=0 
En particulier, si u, — 1 (k = 0, 1, 2, . ..), alors 


p(t) — 2 e"*hP = (4 —ePh)"1, (5.43) 


Le graphe de cette fonction est représenté sur la figure 6. Si œ (£) est 
une fonction k-périodique, la représentation opérationnelle de cette 
fonction sera 
h 
= p(i—e 2") À pe? dt = pqi (p) (1 —eP*y 4. 
0 


L'ensemble de toutes les fonctions h-périodiques sur la section 
0 < t < + © se confond avec l’ensemble de tous les opérateurs de 
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la forme 


(1—e"P*)"1 pr (p), 
où 
h 


A (P) =p | p(t)eP'dt. 
0 


En introduisant la nouvelle variable £* — _ , on peut toujours 
réduire le nombre À à l’unité. Supposons que [t] désigne le plus grand 


Fig. 6 


entier positif non supérieur à {. À toute fonction continue f (f), on 
peut associer la fonction en escalier f ([tl) que l’on notera f {il 
pour abréger. Si le graphe de f (£{) est connu, on construit aisément 


01 234 5 6 7 
Fig. 7 


celui de f {t] (cf. fig. 7). Il est évident que 
ft] = > f (k) (e-"P—e-t0P) = (1 — er) > f(x)e*?. (5.44) 
k=0 k=0 


En posant f{t] el] dans (5.44), on obtient 
eÀlt] — 1=—= jt hp) = PI 
er) 3eme nr. 
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Une dérivation par rapport à À nous donne 


ar = CP) art 1) = 21 
no ap" rs (ep— e})2 ° 
En répétant cette opération, on obtient 
[1 (ft) — 1). (t)—k +1) ACER) ep 1 
k] = (eP— jh . 

ou 

[él{t—1]...[t—k+ipeAlt-h] ep 

a (5.45) 

Désignons spécialement l'opérateur 1 — e? par 
V—i—er. (5.46) 


Si f(t)—0 lorsque 10, alors 
Vilt)= fl —flt—1]. 
Sinon Vf(t)=f[t]—f{t—1] lorsque 121 et Vf{t]—f(0) lorsque 
O<t<<1. S'agissant de l’opérateur inverse + =V" on a 


œ [t] 
ai = 2 e*P{[t] = À 1). (5.47) 


Supposons que 


f(p)=V Df(kje PP et g(p)}=v à g(k)e"*?, 
p. h=0 


donc si 
f(P)=fl4 et g(p)=£g(t), 
alors 
vTHem=t-er Der à 1(e(n), 
d’où 


[t) 


LL (#) 
V'T(p)e(r)= 2 f()elt—H= D fU—k]e (x), 


k=0 
ou 


. - (t] 
f(P)g(P)=V 2 1 [k] gt —K]. (5.48) 
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Ensuite on a 
ft+nl=dt—er À ftnert (4er) À f(r)e-rr-n = 

=0 r=n 


= 7] (p)—(1 — e77) Le”? (0) + 
Len-0pf({)+...+ef(n—1)] (5.49) 
Enfin de (5.42) il s'ensuit 


ph 


LP? © . 
fltA]=(1—e +) 2 fe À, (5.50) 


Les formules (5.46), (5.48), obtenues à partir de la théorie géné- 
rale du calcul opérationnel, auraient pu servir de point de départ pour 
l'élaboration du calcul opérationnel des fonctions à argument entier, 
i.e. des fonctions définies uniquement aux points 0, 1, 2, 3,... Au 
lieu de l’opérateur de dérivation p on considère l’opérateur aux diffé- 
rences V. Son inverse, l'opérateur V”!, désigne la sommation. Au 
lieu des fonctions puissances t* il est plus commode d'introduire les 
fonctions factorielles 


LH —t(t+1)...(t+k—1), 10—1. (5.51) 
On établit sans peine que 
V{HM = on [tm S, a > 1. (5.52) 


D'où il vient 
VAI = n![t]. 
Mais de (5.43) 


eP 
= ; 

donc 

eP [tm _e-P[t+1]m 

VR nt nl : pos 
De (5.53) il résulte pour |A| << 1 

SpA n CITANT 

à (+) = À n | 

0 0 
ou = (1) ll (5.54). 


De (5.44) il vient 


VD vf tu. (5.55). 
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k 
En posant fD=+ on trouve 


e” Alt] 


VERS " 


(5.56) 

Si l’on use des tables des transformations intégrales de Laplace, 
on peut à partir de (5.56) déduire de nouvelles formules. En effet, en 
multipliant (5.56) par f (À) et en intégrant sur À entre O et + oo, 
on aura 


Lee) O0 


v e-vaf (À) = e-Autf (A) dA.. (5.57) 


Cette formule est valable si existe l’intégrale du second membre. Par 
exemple, en posant f (À) = ef, |E] < 1, il vient 


v AS 
ve = (—È) er 


Cette égalité coïncide avec (5.54). On trouvera d’autres exemples 
dans {246], [146]. 


$ 11. Equations aux différences 


Dans la théorie aux différences finies, outre l’opérateur 


Vzltl=z{tl—z{t —-1] (5.58) 

on considère souvent les opérateurs aux différences 
Az (t]=zit +1] — z{t], (5.59) 
éz(=2[t+5]-2[t-5]. (5.60) 


Pour chacun des opérateurs V, À et ô on peut composer une expres- 
sion qui contiendra la variable indépendante {t], l'opérateur aux 
différences et la fonction inconnue. Cette expression est appelée 
équation aux différences. Par exemple, 


F {[t], zft], Axft], . . ., Axft]l} = 0, 


où z {t] est la fonction inconnue, s'appelle équation aux différences 
d'ordre n à une fonction inconnue. Si dans cette équation on exprime 
les valeurs de la fonction inconnue par la formule 


k 
Az{t}= 2 (—1) 7" Cizlt+v], 


v=0 
alors l’équation aux différences d'ordre #7 peut se mettre sous la 


forme 
z{t+n]=O{{t], zft{x(t+1], ...,z{t+n—1}}. (5.61) 
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Sa solution dépend des nr valeurs initiales x (0), x (1), . ..,x (7 — 1). 
Si ces valeurs sont connues, à partir de (5.61) on définira successive- 
ment æ{(n), z (n + 1), ..., etc. La méthode opérationnelle est très 
commode pour la résolution des équations aux différences linéaires 
à coefficients constants 


z{t+n]+azfit+n—1]+...+anaz{(t+1]+aizit]=fft], (5.62) 


sous réserve que soient données les valeurs initiales de la fonction 
inconnue 


(0) =z; TA) = 2x, ::., 20-70: 


Pour passer de (5.62) à une équation opérationnelle posons 


z(p)}=(1—e?) > z(k)eP* = > z(k) er}, 


f)=({—e) 2 fer =v 2 fer, 


et, en portant dans (5.62) l'expression de x [t + n#] (5.49), on obtient 
Z(p) (e"P+ ae... anse + a,) — 
(er) {er ae. Lars) 
(em nP La et2P LE .. en +... 
+ (EP + ae) zu Herr} + f (D). 
En désignant 
Lo(er)=e"lP+aettPL... Ha, ne +a, 
La (eP) = er hr ta er tp LE, Han, 
Lr(eP)=er (k=0, 1, 2,...,n—1), 
il vient 


mx L Li (e n- 
TZ (p) Sn (1 — e?) { EE Zo + TT T1 + .….. += — L CS Zn } + 


f (p) 
+ ACOÉ (5.63) 


Donc, 


1 
z [t] = La /li+A— > Le 2. (5.64) 


En posant 


(1—eP) Ly(eP) 
= LE - = Pr (t), 
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on obtient 
n—1 


= Te + D zxqn (f] 
=0 


Le premier terme est la solution de l’équation (5.62) qui vérifie des 
conditions initiales nulles, œ, Ît] la solution de l’équation (5.62) 
qui vérifie la condition 4 (r) = 0 sir -£k et px (k) = 1. On peut 
trouver une solution zx {t] à partir de (5.63) ou (5.64) en se servant 
directement du tableau et du théorème de convolution (cf. (5.48)). 
On peut également décomposer Le en fractions simples ou se ser- 
vir du théorème de développement. Parfois il est plus facile de trou- 


ver z {t] en développant 7 en une série suivant les puissances de 


e*P, Ce développement commencera par le terme e-"?, i.e. 
a 
L (ep) a 


Si l’on connaît ce developpement, il est aisé de déduire zx [t] à partir 
de (5.64). 


Exemple. Soit à résoudre l’équation 
zt+2]— 22{t+1]+z2flt) = sino il 


avec les conditions initiales x [0] = 0, x [1] = 0. 
Il vient 


sinw[t]=(1—e?) Ÿ sinwke* — 
k=0 


L(eP) = 627 — 2e? + 1 = (eP — 1). 


Donc (cf. (5.63)), la solution s'écrit 


(eP — 1) sin w 
e®P — 2eP cos +1 ? 


sin @ 


TUE Gr rc 0 1 : 


Décomposons cette expression en fractions simples 


sin & ___ À Fa B C 
(eP—1)(e2P—2eP coSw@+1)  eP—1 eP—ei® eP—e-iw 
D'où 
A=— sin @ D — 1 — cos © + i sin © 
7 2(1—cosw) ” = &i (1 — cos &) - 


C=— Î — cos © — i sin w 
7 4i(1—coso) 
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Compte tenu de (5.44) et (5.47), on obtient 
_ [e) (t-1) 
1 _ ee eP—1 = er Ÿ et] — D ex 
h=0 k=—0 


pour & > 1. Donc, lorsque t > Ÿ on a 
[t-1 
sin @ + (1 — cos w) sin kw + sin © cos ko 


z|tl= ZU— cos a) | 2 (1 — cos w) 


k=0 
et zx (0) = 0. 

De façon analogue on peut résoudre un système d’équations aux 
différences linéaires à coefficients constants [89], [146], (2471. 


$S 12. Transformation d’Efros 


Pour la réalisation d’un opérateur on a parfois intérêt à se ser- 
vir de la transformation d’Efros (cf. également chapitre II, $ 3). 
Plus exactement, si F (p)=® (ft) et e-°1P) u (p) q (p) = © (E: t), 
alors 


u(p) FlG)N= | oO: 1) dE. 
0 


Dans {64], Efros et Danilevski ont indiqué d'’intéressantes applica- 
tions de cette transformation. On peut se servir du théorème suivant 
pour justifier la transformation d’'Efros. 

Théorème 8. Supposons que 


Co 


1) | Ip(é)e-sot dt < oo, O0 >. 


0 

2) q (p) est une fonction analytique, régulière dans le demi-plan 
Re p > ©, et vérifiant la condition Re q (p) > ©, lorsque Re p > 
> 01 > Co: 

Alors 


= TE Fi WE; t)p(®)dE8, (5.65) 


où 
W(G;t) =<+e- 8a(p). 
De (5.65), il résulte 
Fla(p}= par) | VE; 1) p(E) dE. (5.66) 
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Si l'opérateur F [q (p)l se ramène à une fonction, on peut trouver 
cette dernière en se servant de (5.66). 

L'opérateur e-*1) est souvent ramené à une fonction. Posons 
e-F4(P) —p(E; ft). Alors, au lieu de (5.66) on aura 


00 


Fla(pn=a(p) | bG: 9 pdt. (5.67) 


0 


Dans certains cas, l’opérateur q (p) peut être mis sous le signe d’in- 


tégration. 
Au chapitre IX, $ 1 figurent des exemples sur la transformation 


d’Efros. 


$ 13. Equations différentielles opérationnelles 


On appelle équation différentielle opérationnelle d'ordre n une 
expression 


F{A; z(À), z'(4), ..., 2 (2)} =0 


contenant la variable indépendante À, & << À << B, la fonction opé- 
rationnelle inconnue z (À) et ses dérivées x’ (À), x” (À), . . ., xt" (A). 
Même les types les plus simples de ces équations, par exemple les 
équations linéaires 


An (A)T (À) + an (A) TP (A)... +ao(À)z(A)=f(À), (5.68) 


où les coefficients a, (À), a1 (À), . . ., an (À) sont des opérateurs dépen- 
dant de la variable réelle À, & << À << B, sont peu étudiés. Donc, ici 
on ne considérera que des équations différentielles linéaires à coef- 
ficients constants 


ant (À) + an220 9 (À) +... + ao (à) = f (À), (5.69) 


OÙ Gp Es + + + An SOnt des opérateurs constants arbitraires. En parti- 
culier, si les opérateurs ay, &1, - .- ., an et f (À) sont des nombres, 
l’équation (5.69) se transforme en une équation différentielle ordi- 
naire d'ordre nr. La méthode de résolution de l'équation (5.69) coïn- 
cide avec les procédés standard de résolution des équations diffé- 
rentielles linéaires ordinaires à coefficients constants. D'abord on 
commence par trouver la solution générale de l’équation homogène 


ar Hanna) Li, Hat =Ù (5.70) 


sous la forme d’une fonction opérationnelle exponentielle zx (4) — 
— el, L'opérateur w est déterminé à partir de l’équation caracté- 
ristique 

anw"+ a, uwTi+.,..+ap=0. (5.71) 
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Soient w:, Wa, . . ., w, les racines de l’équation (5.71). Formons les 
fonctions exponentielles 
ent, eMM2 en. (5.72) 


Si ces fonctions sont construites, la solution générale de (5.70) sera 
l'opérateur 


z(A)= D Creer, 
k=1 


où C, sont des opérateurs arbitraires constants. Observons que a) les 
équations de la forme (5.71) ne sont pas toutes solubles sur le corps 
des opérateurs M, b) la fonction exponentielle e” n'existe pas pour 
tout opérateur w. 

L'opérateur w s'appelle logarithme si existe la fonction exponen- 
tielle ex. Par exemple, les opérateurs p et | p sont des logarithmes, 
alors que ip ne l’est pas. Si l’opérateur est logarithme et racine de 
multiplicité r de l'équation (5.71), alors chaque fonction e”, 
Aelw, ...,A"-lelw est solution de l’équation (5.70). On distinguera 
trois types d'équations différentielles selon le nombre de racines des 
logarithmes. On dira qu’une équation différentielle est 

1) logarithmique si toutes les racines de son équation caractéris- 
tique sont des logarithmes ; 

2) pure si son équation caractéristisque ne possède aucune racine 
logarithme ; 

3) mixte si existent des logarithmes parmi les racines de son 
équation caractéristique. 

En général, le nombre des opérateurs constants arbitraires figu- 
rant dans la solution générale est égal à l’ordre de l’équation. S’agis- 
sant de l'équation (5.69), on a le théorème d’unicité 

Théorème 9. Etant donnés des opérateurs V,, Vi, ..., Va 
et un point À, compris à l’intérieur de l'intervalle a < 9 << B il 
existe une fonction opérationnelle x (À) et une seule qui vérifie dans 
lu, BI l'équation (5.69) et Les conditions x (ko) = Vos &’ (ho) = Vas - - - 
anim No) Vo 

Toute solution de l'équation (5.70) peut être obtenue à partir de 
la solution générale moyennant un choix adéquat des opérateurs 
Cis Cas + + + Cn- La résolution de l’équation non homogène (5.69) 
se ramène à celle de (5.70) si l’on connaït au moins une fonction 
z, (À), solution de l’équation (5.69). Donc, pour résoudre l’équation 
non homogène, il suffit de trouver une fonction vérifiant (5.68), 
(5.69). D'une façon générale, il est difficile de trouver cette fonction. 
Dans certains cas, elle n’existe pas, dans d’autres, on peut la trou- 
ver assez facilement si le second membre de l’équation est de forme 
spéciale, en particulier, si f (À) est un polynôme ou une fonction 
exponentielle [165]. 


$ 14. Application du calcul opérationnel 
à la résolution des équations différentielles 


Soit l’équation différentielle ordinaire d’ordre rx à coefficients 
constants 


2 (t)+art(t)+...+anr()=f(t), (5.73) 
0Lt< + 00 


avec les conditions initiales x (0) = x,, x” (0) = 21, . . ., z"-b (0)— 
= Zn. Appliquant la formule (cf. (5.11)) 


20 (9 = pl [2 (02022 0. 2-0 (0)], 


on peut mettre l'équation (5.73) sous la forme 


nl) — 
L (p) [x (—%—< ma | 

a __ bi _ bn-1 

= f(t)—b : Le FT 
Tci 

n—1 
L(p)=p"+mp" + ... + an, bn = 2 Tslnth-s 
(k=0, 1,..., n—1), 

d’où 


n— 1 


OST OT tee + er — - (5.74) 


Cette formule donne l'expression de la solution de l'équation (5.73). 
On constate aussitôt que le second membre de (5.74) est une fonc- 
tion n fois dérivable, vérifiant les conditions initiales. Si les condi- 
tions initiales sont nulles x, = x, = ... — zh — 0, la solution 
prend la forme remarquablement simple 


1 
2(t)=r I). 
Les équations différentielles à argument retardé se résolvent de 


façon aussi simple lorsque leurs coefficients sont constants 


n—1 
2 ()= 2 ar (t—h)+g(t)  (0OLSt< +oo, h > 0). 
k=0 
Pour simplifier on supposera que les conditions initiales sont trivia- 
les, i.e. lorsque t < 0 on posera x (t) = x’ (t)= ... = 2" (ft) = 0. 
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Puisque z% (t — h}) = pe-"kPzx (t) la solution est 


D (= — "À 2 z (+). 


_ > agp'e "RP 


On montre que le second membre de cette expression est une fonc- 


tion nr fois dérivable, satisfaisant aux conditions initiales trivi- 
ales [50]. 


Considérons une équation aux dérivées partielles dont les coeffi- 
cients a. (x) sont des fonctions numériques de la variable z 


S 5 uv (7) ED 2j (7, ?). (5.75) 


re Oz" ot" 


Appliquant la formule (cf. (5.11)) 


OHFVu(z.t) _ _, dPu(z.t) __  v dPu(z, 0) 
EEE ôxt zx 
uns 2 Flu(r. 0) __,. OV iu (x, 0) 
P Oz 0 ot 


on ramène l'équation à la forme 


Ç Mu (z. 1 D < 9 +hu (z, 0) 
u(T, UT, 
DEATCE Du f(z t)+ Ÿ > > Don 

er un=0 v=1 k=0 


où du = au(zT; p}= 2 duvy (x) p'. Désignant le second membre de 


cette équation par O(z; p) et considérant u(z, t) comme une 
fonction opérationnelle dépendant du paramètre x, u(zx, t)=— 


=u (zx; p}=u(zx), on obtient 

Au) (x) + amau "0 (z)+ ... +au(z)=O(z; p). (5.76) 
Les coefficients a, sont ici également des fonctions opérationnelles 
dépendant de x. Donc, l’intégration de l'équation (5.75) se ramène 
à celle d’une équation différentielle opérationnelle linéaire. L'équa- 
tion (5.76) s'appelle équation transformée. En résolvant l'équa- 
tion (5.76), il importe de se servir de l’isomorphisme des corps 
et JR. Dans le corps 9, l'équation transformée (5.76) devient une 
équation différentielle linéaire ordinaire d’ordre x dont les coeffi- 
cients et le second membre dépendent d'un paramètre p qui est un 
nombre complexe. Ces équations sont bien étudiées. Soit u (zx; p) 
une solution de cette équation. Si u (x; p) € M pour des valeurs 
8—0394 
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données de x, a << x << $, cela traduit le fait que l’équation (5.75) 
possède dans le corps M la solution u (x; p), où pest l’opérateurp =. 


Appliquer le calcul opérationnel à la résolution des équations 
aux dérivées partielles, c’est : 

1) Remplacer l'équation initiale par l'équation transformée. De 
façon analogue les conditions aux limites sont remplacées par les 
conditions aux limites transformées qui seront conditions aux limi- 
tes pour la solution de u (x; p) de l’équation transformée (5.76). 

2) Trouver la solution u (x; p) de l'équation transformée qui 
vérifie les conditions aux limites transformées. 

3) Etudier la solution obtenue dans le but d'établir son apparte- 
nance au corps MW. Si uw (x; p) E M il faut voir si la solution n’est 
pas solution généralisée, ou si elle peut être ramenée à une fonction 

, CRC e « . , OMAN u(z,t) 
possédant des dérivées partielles par rapport à z et t jusqu’à ETS 
comprise. Ce dernier fait signifiera que u (x, t) vérifie l’équation 
initiale aux dérivées partielles au sens classique. 

4) Réaliser l’opérateur u (x; p), i.e. définir la fonction u (x, t) — 


u (zx; p). 

L’étude du 3) peut être singulièrement simplifié si 4) est réalisé. 
5) Démontrer que la solution w (x, t) vérifie les conditions initia- 
les et aux limites du problème. 

A titre d’exemple, considérons les équations 


p (x) ur = Po (T) Uxx + P1 (x) ux + Pe (x) u, (9.77) 
P (z) ur = Po (ZT) Uxx + Pa (T) Ux + P: (x) u (5.78) 


dans le domaine 0<z<1l,t>=>0,p (x), po (x), P:1 (x), P2 (x) étant 
des fonctions continues données sur l'intervalle 0<z< et 
p (x) >> 0. La solution u (x, t) doit posséder dans le domaine (0 < 
<z<l,t>>0) des dérivées partielles continues jusqu’au second 
ordre compris et satisfaire aux conditions initiales 


lim u(z, t)=p(x), 0<z< Il, 
t-+0 


dans le cas de l’équation (5.77) et 
limu(xz, t)=œ(x), limu,(z, t)—=v(x), 0O<t< I, 
t+0 t—+0 
dans le cas de l’équation (5.78), ainsi qu’aux conditions aux limites 


lim u(z, t)—#f(t), au,(l, t)+bu;(l, t)=cu(l, t) (5.79) 
X— +0 


pour { >> 0, où p (x), Ÿ (x) sont des fonctions continues par morceaux 
données; f (t) € S et est continue pour t > 0; a, b, c des constantes 
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données. On cherchera la solution de ces équations sous la forme 


u(x, t)=u(x; p). Les équations transformées de (5.77) et (5.78) 
eee 


Po (2) 2 +9 QE M ps(z)— pp(nlu= —p(x)po(z), (5.80) 


Po (2) + dx re Pa (x) e + [Pe (x) — p'o (x) u = 
= — pp (z) p(x)— pe (x)Ÿ (x). (5.81) 


Les conditions aux limites du problème nous donnent les conditions 
aux limites de la solution 


u(+0; p}=f(p), où f(P)=f(). 
au, (l; p)+bplu(l; p}—q(l)}]=cu(t; p) 

Théorème 10. Supposons que u (x; p) est solution de l'équa- 

tion (5.80) ou (5.81) avec la condition (5.82). Supposons que par ailleurs: 


41) Les opérateurs u (x; p), u, (zx; p) et u,. (x; p) sont réalisables 
pour 0OLr< Il. 

2) Qu’existe un nombre ©, tel que lorsque t —+ © soient remplies 
les conditions 


u (x : P) = 0 (e70'), U, (z ; P) = 0 (e°'), pe (z ; P) = 0 (e°0°) 


uniformément en x sur tout intervalle EL z< l. 
3) Qu’existe un entier k > 0 tel que 


(5.82) 


[pu (x; p}l << Q= const 


pour tous less0<r<e<l, Re p > 6; > 60. 

4) Qu'existe lim u (x; p) = g(t), t > 0 où g (t) est une fonction 
continue pour T0 et bornée lorsque t —0. 

Alors u (x; t) = u(x; p) est la solution de l'équation (5.77) ou 
(5.78) qui vérifie les conditions aux limites et les conditions initiales. 

Examinons un exemple. Trouver la solution de l'équation 


Ur = xx 


qui dans le domaine (0 £< x < 1, t > 0) vérifie les conditions initia- 
les u (x, 0) = 0, u, (x, 0) — 0 et les conditions aux limites u (0,1) — 
= 0, Eu, (l, t) = À sin ot, où £, « et À sont des constantes. 

L’ équation transformée s'écrit 
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les conditions aux limites pour uw (zx; p) sont 


a 4 _: — . __. Op 
lu (0 ; p)=0, Eu, (£; p) = p° + w e 
La solution de l’équation transformée qui vérifie les conditions aux 


limites est de la forme 


px 
sh — 
| a «a À 
u (zx; P) = pi w® pl ? b— E ? 
a 
d'où 
: pie Aa a 2b x (— 1}? sin khz-sin khat 
us eu DE TR T2 RE 
ca n=1 
où 
Ed 1 
ka=+(n——) (n=1, 2,...), © akh. 


$15. Séries asymptotiques 


D'après une définition de Poincaré, on dit que la série >, u, (2) 
0 


n= 
est le développement asymptotique d'une fonction s (z) dans un domaine 
donnée des valeurs de arg z si, quel que soit 


N 

s (z) — > un (2) 
lim — = 0 
]zi- 00 UN (z) 


En particulier, la représentation asymptotique d'une fonction f (t) 


An 


NT voudra dire que, quel que soit W, 


par la série entière Y» 
n —t{ 
N 
. 4N An \ 
lim 4 (1@— 5 Æ)-0. 
n=—0 


Le développement asymptotique est d’un usage commode lorsqu'il 
s’agit de calculer les valeurs d’une fonction pour des valeurs élevées 
. de l’argument. Aussi dans le calcul opérationnel des opérateurs se 
ramenant à des fonctions est-il important de construire le dévelop- 
pement asymptotique d’une fonction d’après sa représentation opé- 
rationnelle. Dans de nombreux cas, la représentation asymptotique de 
fonctions transformables-Laplace peut être obtenue à l’aide du 
théorème suivant [246]. 
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T hé orème Â1. Supposons que f (p) =f (t) et 

1) f (p) possède des points singuliers isolés (des pôles et des points de 
branchement), 

2) dans le demi-plan Re p< 0 la fonction —— 


vers zéro p en arg p lorsque |p| — oo, 
3) le nombre des points singuliers p — p, de plus grande partie 


Ê 2 | tend uniformément 


réelle est fini (s — 1, 2,..., let le développement de Fa Æ ou voisinage 
de p = p, est donnée par la série 


sd 2 (9) 
2 4 (p— pa) (NL ENT ELU oo). 


Ceci étant, la D asymptotique de f (t) sera 


(0) — 5 e°* 5 


(5),,° 
s={ h=—0 LP (— Sa : 


1 ms (s) 
Ta © lorsque Àx” —0, 1, 2 
Si, par exemple, tous les points singuliers de la fonction À 24 , ) 


à l'exception de p = 0, possèdent une partie réelle négative et que 
l'origine des coordonnées soit un point de branchement de premier 


ordre et enfin que —— Î Re Î) admette au voisinage de p = 0 le dévelop- 
pement 
F(p) _< D E 
P EE > GnpP 9 
DEES | 
alors 
HO — 
n=! 2 
r(—5+i): 
ou 


LR (14:35 … (261) 
f(t) à Va QU Zoh+1- 


4 


$ 16. Calcul opérationnel pour l’opérateur == tu 


La théorie du calcul opérationnel repose sur la notion de convolu- 
tion [165]. L'introduction de la convolution dans l’ensemble liné- 
aire M (cf. $ 1) confère à ce dernier une structure d’anneau commuta- 
tif sans diviseurs de zéro. Toute la théorie du calcul opérationnel de 
Mikusinski est basée sur ce fait. Dans ce paragraphe, en partant 
d'une nouvelle définition de la convolution on se propose d'élabo- 


PR ; d ,d ; 
rer un calcul opérationnel pour l’opérateur B = le On procédera 


exactement comme pour l'opérateur p — . Soit L; l’ensemble de 
toutes les fonctions f (t) définies sur la section [0, + œ{, intégrables 
sur tout intervalle fini de cette section et vérifiant la condition 
t È 
dd À = Q: 
Î = | 1/1 dt<o (5.83) 
0 0 
quel que soit 4 >0. Par Mer la fonction FA =Int appartient 


— 1)° 
—, Ê EL. Soit Af; 


l'ensemble des fonctions de la forme 
; : 

; dE 

FH=(+ 

» 


& 
0 0 


à L;, contrairement à f (6) =“ bien que 


en? 


f(u) du+C, 


où f ({) est une fonction quelconque de Z, et C une constante arbi- 
traire. Définissons dans l’ensemble Z, la convolution des fonctions 
f1 (t) E La et f: (t) € L1 à l’aide de la formule 


{ 


= | d8 | fi (n) fe LG — n) CE) dn. (5.84) 


U 


On montre que la fonction f (t) appartient à L,. L'ensemble M, est 
un ensemble linéaire. Toute fonction appartenant à l’ensemble M, 
possède presque partout une dérivée seconde F” (t). Appelons produit 
des fonctions F, (t) et F, (t) de M, l'expression 

{ 1 


FOX F0 = {14 | dE | Fi En F1(4—n) (1—H1dn. (6.85) 
Si 


t 


Fi (Et) — _ | fh\ii)dn, Fa(t) = 


de 
EF fe (n) dn 


Cm ES 
Otis 
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et 
t 1 
F0) = | 4 | fi En) fa 1(1— 1) E— Hd, 
0 0 
on aura 
{ u 
Fit) x F2(0) [ æ Î f (v) dv. (5.86) 
0 0 


De là il s'ensuit que le produit de deux fonctions de 7, appartient 
de nouveau à M,. L’addition se définit de façon naturelle sur M. 
Il est aisé de vérifier que le produit défini au sens (5.85) est commu- 
tatif, associatif et distributif relativement à l’addition. Donc, l’en- 
semble M, possède une structure d’anneau commutatif relativement 
à l'addition et la multiplication. Cet anneau ne possède pas de 
diviseurs de zéro. 

L'anneau M, peut donc être étendu à un corps de quotients 
que nous noterons Y},. Les éléments de 9, seront appelés opéra- 
teurs et désigné S 1 D 1 désig 

| gnés, comme au $ 1, par --. Donc, - désigne une 


classe de couples équivalents. L'égalité des opérateurs _. et E. 
signifie que F, XX Ge = Fe X Gi. Si dans (5.85) la fonction F, (t)=«, 
où « est une constante, alors 


a X Filt)=a {Tai nœman)- 


t 
Lx + + | ra} er, 0 


Donc, sur M, le produit d’un nombre par une fonction coïncide 
avec le produit ordinaire d’un nombre par une fonction. Si dans (5.85) 
les deux facteurs sont des nombres, alors le produit au sens (5.85) 
est confondu avec le produit ordinaire de deux nombres. De là il 


résulte que les rl de la forme a peuvent être identifiés 


aux fonctions — LF (£); en particulier, les M Je ) se confon- 
dent avec les fonctions F (t) de Mi. Dans ce cas, on rire 


FO LF(), F(HEM. 
Enfin, les opérateurs de la forme 0, F (0) = 0, s'identifient 


à des fonctions de ZL;. Plus exactement, à tout opérateur "a, 
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F (D) EM; et F (0) = 0, associons une fonction f (t) telle que 
F F0 


= f (£), (5.87) 
où | 
F (t) = Î £ f (u) du 
0 0 


Réciproquement, à toute fonction f (t) € L, correspond un opérateur 
"0 . Cette correspondance est biunivoque. Si un opérateur est sus- 


do d’être ramené à la forme ©, F(t) EM;et F (0) =0,0on 


dira que cet opérateur est une TR Désignons spécialement 
l'opérateur Z par 


LB. (5.88) 
Pour l'opérateur inverse Bi, on aurà 
+ =t. (5.89) 


Compte tenu de (5.87), on obtient pour les fonctions F (ft) E M; et 
F (0) = 0 


d dF ’ 
BP (= (t 2) (7 (0) 
ou 
BF(t)=tF"(t)+F" (+). (5.90) 
Donc, lorsque F (t) E M, et F (0) = 0 le produit BF (t) signifie 


que l’on peut appliquer à F (t) l’opérateur LT. Le produit de 


l'opérateur inverse _ par la fonction f (t) € L;, ainsi qu'il résulte 
de (5.89) et de l'égalité F (t) = t % f (t), est égal à 


t ë 
+1 = [= | f(u) du. (5 91) 
0 0 
De (5.89) et (5.91) il s’ensuit que 
1 in : 
En — TUE à (5.92) 


Alors de (5.85) et (5.92) il résulte pour f (t) € L: 


1 


4 1 : _ 
mnf® = | (— 2)" & | En) ({—m"an. (5.93) 


0 
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L'équation 
(+ )=a (5.94) 
dt dt ° 


admet pour solution les fonctions de Bessel 7, (2V At) et K, (2V ht). 
De (5.90) et (5.94) il s'ensuit que 


B{lo(2V Rt)—1]= A0 (2 V ), 
d’où 


B 
3x = lo (2 VA). (5.95} 
En posant À = io dans (5.95) et compte tenu de 
I, (2 V iwt) = ber (2 V wt) + i bei (2 V'wt), 


on trouve 
= ber (2 Var), LÉ = bei (2 V of). (5.96) 
De (5.95), on déduit 
B _— 
Bix vo (2V M). (5.97} 
(5.95) et (5.96) “és 


= + Mo (2 VM)+ Jo(2V D): | 


=7 
= (20) —Jo(2V MI. 


La plus grande partie de la théorie du calcul opérationnel, exposée 
dans les paragraphes précédents, est valable pour le corps des opé- 
rateurs Î,. En particulier, la définition de la fonction opération- 
nelle, de la limite d’une suite d'opérateurs et la notion de série opé- 
rationnelle, de la dérivation et de l'intégration d’une fonction 
opérationnelle se transposent sans changement au corps Yt,. En 
se servant de cette théorie on peut par des méthodes connues compléter 
la table des valeurs des opérateurs (5.95) à (5.98). Par exemple, en 
dérivant (5.95) par rapport au paramètre À, on trouve 


n 


(5.98) 


B=— 3: 


B { JL \" — 
NT (+) ZA CV. (5.99) 
De (5.96), il vient 
B B*° —_— 
B? w2B oB : pere 
B (= —1) = —HEw — —o bei (2 V ot). 
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d d = : — 
rl + ber (2 Vut) } = —« bei (2 Vwt). 
Donc, les fonctions ber (2V/wt) et bei (2V wt) se conduisent vis-à- 


vis de l'opérateur Jt3; exactement comme cos ot et sin wi vis-à- 
d 


vis de l’opérateur L. Pour l'opérateur B — = t-on peut construire 
le calcul opérationnel à partir d’une transformation intégrale. L’ana- 
logue de la transformation de Laplace sera ici l'intégrale 


_ {15 bei (2V wi) } = o ber (2 V'wt), 
(5.100) 


f" (B)=2 | f (4) K(2V Bt)dt. (5.101) 
0 
Si f (t) € L: et vérifie la condition : 
| f(t)1<Qe2vVt, (5.102) 
où Q et q, > 0 sont des constantes, alors dans le domaine Re VB >4% 
l'intégrale (5.101) est absolument convergente et représente une 


fonction analytique de la variable complexe B. Moyennant des hypo- 
thèses connues sur la fonction f (t), on a la transformation inverse 


= | f (8) (2 VB dB, (5.103 
L 


où le chemin d'intégration L est une parabole quelconque Re VB = 
= 4 > %o- La formule (5.103) aura lieu, par exemple, si outre (5.102), 
la fonction f (t) est assujettie à être à variation bornée au voisinage 
de tout point { = t, de la section 0 << {<< + co. Aux points de 
discontinuité, l'intégrale (5.103) est égale à + {f (£o+ 0) + f (to—0)}. 
Pour la transformation (5.101) on a le 

Théorème 12. Soient f.(t)et f.(t) des fonctions de L, véri- 
fiant la condition (5.102) (chaque fonction possède ses propres constantes 
Q et go dans (5.102)). Alors la fonction 


t 1 
f()= | 48 | HŒEm/AlA—N)(C—H)dn (6.10 
0 (Ù 


appartient à L, et remplit également une condition du genre (5.102). Si 


ft(B8)=2 | f (0) Ko (2 V Bt) dt, 
0 


1 (8)=2 | f:(0 Ko (2 V Bt) dt 
0 
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et 
f°(B)=2\ f(# Ko (2V Bt)dt, 
0 
alors 


f* (B)= f? (B) f: (B). (5.105) 


Définition. Si pour un opérateur a € 9h, il existe un re- 
présentant (F, G) tel que les fonctions F (t)et G (t) admettent la trans- 
formation (5.101), on dira que cet opérateur est transformable- 
Bessel et la fonction 


TF6 Ko(2 V Bt)dt 
a (B) = ——— (5.106) 
| G(4) Ko (2 V'Bt)dt 
0 
F 


transformée de Bessel de l'opérateur a=— 7 - Notons 1 l’ensemble de 


tous les opérateurs transformables-Bessel et 9, l’ensemble de leurs 
transformées de Bessel. Ces ensembles sont isomorphes entre eux. 


Cet isomorphisme envoie l'opérateur B — _ dans la fonction 


Ÿ ÆKo(2 V Bt) dt 
_ = B. 


Î {Ko (2 V Bt)dt 
Ô 


Dans les calculs on peut se servir de la formule 


2 | MK (2 V Bt) dt LUN (5.107) 
0 


pti 


Si l'opérateur a se ramène à une fonction, i.e. 


a= +0, F(HEM, F(0)=0, 
alors 


2 T F(D Ko (2 V En à c0 
AB) = — = 2B? | F(t)Ko(2V Bt)dt. (5.108) 
2 \tKo(2 V Bt)dt 0 
0 
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Comme 
t £ 
F (t) — ( al jf (u) du 
0 


et compte tenu de relations connues auxquelles satisfait une fonc- 
tion de Bessel 


GR (G)=—zk() ot (Ko) = —Ki(a) 


un double intégration par parties nous donne 


oO 


a(B)=—2B Î f(t) Ko(2V Bt) dt. (5.109) 


U 


On reconnaît un analogue de la transformation de Laplace-Carson. 
Moyennant cette transformation, on peut déduire de nouvelles for- 


mules pour l’opérateur B— Lt _. . L'égalité (5.109) veut dire que 


a (B) = f (t). (5.110) 


En particulier, si f (t) = {*, alors de (5.107) et (5.108) on déduit 
une généralisation de la formule (5.92) à tout ordre v 


1 +” 


= = FUL Rev>—1. (5.111) 
Si f(t)=e""*, il vient 
Be. 
Les | — du = e"tt. 
B 


Citons quelques exemples d’application de la formule (5.108). 
Des calculs directs nous donnent 


O0 


2B | K,(2 VA) Ko(2V Bt) dt = In y #. 
0 


Donc, 


B—X 
En utilisant l'intégrale connue 


B_ny/Eek (2m. (5.119 


© 


vtt gr — (0) T(v+1) 
Ï Ke (az) J (be) 24 dr 


d = 
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on généralise la formule (5.99) à un ordre quelconque v 


B _ 1 42 91/73 
D CE x } J, (LV At). (5.113) 
En particulier de (5.113) on obtient pour v = + 
B____ces2 vit A 
5.114 
VBEX a Vt ) 
En dérivant (5.111) par éd à v, on obtient 
lbB (+1) 
FA “HET [2 T(V+1) Int |. (5.115) 


L'intégrale 


| Ko(aVz?+2)z2u+i dr = 24 (2) "rA+) Ki+y (az). 
0 


entraine 
u+1 


’ sad " _—_— : À 2 Eds 
ro 2 | é— A)" Ko(2V Bt) dt=2B (+) Kim (2 VB), 


d'où 
à _. … 0 si t<X, 
2B(+)  Kim(2VB)= 4 pi ,, 6-10 
TRE ST 
ne 
En particulier, lorsque u = — 7. on déduit de (5.116) 
Er. 0 si t<X, 
Be-2VBr — 5 414 
V d = si À <t. ê c 


Multipliant la dernière égalité par une fonction arbitraire œ (à) 
(au sens ordinaire) et intégrant sur À entre 0 et co, on obtient 


VE T e-2 Ve qait | ed (5.118) 
Posant 2V à À = Ë dans l'intégrale du RS membre, il vient 
VB Ï e-iVEop (= )Eæ= 2 AE (5.119) 
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En se servant des tables de l’intégrale de Laplace, on obtien- 
drait à partir de (5.119) de nombreuses nouvelles formules. Par 
exemple, si dans (5.118) on pose 


J, (2 V ai) 
on trouve 
. = 117 EMA) D à 
B e-iVBJ,(EV a) R=+ | ViE=N 
ou encore 


B_{ EC =+f Je Va) }; 
= = L VA (—2) eù ° 
Mais la dernière intégrale qui est une convolution au sens ordinaire 


peut être calculée par la transformation de Laplace (cf. [64] page 125). 
On aura 


FE TETE Ÿ = 73 (Vai). (5.120) 


En particulier, lorsque v = 0 


= = JE(V ai). (5.121) 


En multipliant (5.121) par VE =. il vient 
B 2 = 1 
—— —J; = 
V B+a (V'at)+ 11" 


ou compte tenu de (5.114) et (5.85) 


t 1 
cos2Wat _ d d r2(V/ En 
ae (ral dE EL (V'aën) 


Cette égalité peut être mise sous la forme 


TS = (t —E) }- 


{ Î 


2Vat 
sin u du= | à ef J8(V an) an 


| VA n) GE —E) 


2 

Va 0 : 0 
De façon analogue, en élevant (5.121) au carré et tenant compte 
de (5.97), on obtient 


= { | 
10 OV a) D À à À JE (Van) IV a QG 0 EH] ên. 
0 


0 


d 
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En partant de formules identiques à (5.118), on obtiendrait de 
nouvelles relations entre les fonctions spéciales. 
Calculons à titre d'exemple la somme d’une série fonctionnelle. 


Soit k = VB + 1 — VB; on a l'identité 


œ 


HE (264264 


Compte tenu de (5.120), il vient 
Ji(Vt)—275 (Vt)+ 25 (VE) —275 (V5) +... =J0(2Vt). 


Multipliant l'égalité (5.99) par u" et sommant sur n de 0 à , on. 
trouve 


B+1° 


F ___ By B 
(B— (B— pi — — 2 —u 2 


d'où 
LS 
2 


SE (1) 7, (2VH)= (2 ATH). 


n—0 


En conclusion, établissons un lien entre la table de valeurs des 
opérateurs F (B) et la table de valeurs des opérateurs F (p), où 


p=+. Soit F(B) = pi) et F (p) = f (#). Par définition des opé- 
rateurs F(B) et F(p), on a 


F()--28 | eo K(2VFD. 


CO 


F(p)=p | f(#ert dt. 
Donc, ÿ 


2p \ P (t) Xo (2Vpi)at=p\ f (t) er! 
0 


0 
Compte tenu de l'égalité 


œ dt de _ 
(e = -z —2K9(2V pat), 


9 
où il importe de poser q = 1, il vient 


[ p(£) dt La = ET ertd 
0 0 0 
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| e"l'< dx [ p (Ex) e-6 dé — | f(the-"'dt, 
0 U 0 
d'où 
1(= | ptbe-tdæ. (5.122) 


0 


Ici o (t) appartient à L, et vérifie la condition (5.102). Ceci nous 
suggère le procédé suivant de calcul de F (B) à l’aide des tables de 
valeurs de F (p). Donc, pour trouver F (B), il faut calculer F (p) = 


— f (t), puis remplacer l'argument t de la fonction f(t) par _ et 


déterminer la valeur de l'opérateur f ( =) —  (t), ce qui donne 


F{&(B) = œ (t). Ce procédé n’est certes valable que dans le cas où la 
fonction f (t) est représentable par l’intégrale (5.122). Appliquons 


Ld # Pr B 
ce procédé au calcul, par exemple, de l’opérateur 5: On a 


p 

a Me 
—e"#,er = I, (2V at). Donc, D — — 1, (2V'at), i.e. a lieu (5.95). 
De façon analogue, on trouve 


= VAE —B) = av], (2 V at) I, (2 V'at), (5.123) 


BinËte=y/£{7,(2ya)—n(2Vat)l. (5.124 


Posant vdi dise il vient les identités 
B 


VE Fee 2(B2+1) ? 
BP 
92 me ee 


En substituant aux opérateurs leurs expressions (5.96) et (5.124), 
on trouve 


J(2Vt)A(2Vt)—-J:(2Vt)1:(2Vt)+ 

+ JS (2VT)15(2V7)—J:(2V7) 1 (2VT) +... = + bei (25), 
J(2Vt)L(2Vt)—-27,(2Vt)2(2Vt)+ 

+927, (2Vt)Z(2Vt)—2J,(2Vt) 1 (2Vt)+...—=ber(2Vt). 
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Avant d'achever ce bref exposé sur la théorie du calcul opération- 

9 La d d 9 9 : 
nel pour l’opérateur B — ton observera qu’en plus de l’appli- 
cation de cette théorie au calcul des intégrales définies et à la somme 
de séries fonctionnelles, on aurait pu citer des exemples d'intégration 


d'équations différentielles. L’équation différentielle la plus simple 
à résoudre est 


L(+t- )s()=f (6), (5.195) 


où L'(4) = À" + at +... + a, est un polynôme à coefficients 
constants a;. Si les conditions initiales sont nulles, i.e. 


“A (0) = 0, Bz (t) |t=0 = [tx” (t) + Z (£)]:=0 = 0, 
B?z (t) |1=0 = (2x7) (4) + 4tz" (t) + 2x” (t)]=0 = 


Vs Ho « A fée =0 — 0, 
en substituant à l'opérateur — — tel opérateur B l'é SE (5.125) 


s'écrit L (B) x (t) = f (t) et sa sr sera x ({) — f(t). On trou- 


AE) 
ve l'expression de TE f(t) en décomposant TT en fractions 
simples de la même façon que dans le calcul opérationnel de Heavi- 
side. Si les conditions initiales x, ne sont pas nulles, on fait la subs- 
titution suivante (cf. 5.90)) 


d d 
pri LA re À (1) = B [x (£) — Lo]. 
L'application successive de cette formule donne 


(+ + e—) z(t)= B“x (t)— Brya— Bas — 0e... — B'xo. 
Exemple: soit à résoudre l'équation 
(tr) z 31 z(t)+27(t)=0 5.126) 
at) ?( a a © () = (5. 


avec la condition x (0) = 0, Bz (t)|:.0 = &. Mettons cette équation 
sous la forme opérationnelle 


B°z(t) — Ba — 3Bz(t) + 2r (t) = 0 
(B? — 3B + 2)zx(t) = «B. 


ou 
D'où (cf. 5.95)) 
D) = ET = DT pre (0 (2V 24) —2(2Vt)]. 


On résout de façon analogue certains types d’équations aux 
dérivées partielles. 


9—0894 


$ 17. Sur une généralisation du calcul 
opérationnel 


De ce qui précède il suit que le calcul opérationnel pour l’opéra- 
teur p = _ peut être construit par une extension du corps des fonc- 


tions sur lequel le produit est défini par 


F (t)»G (+) =+ | F (t— Et) G(E) dt. (5.127) 
0 


Le calcul opérationnel pour l'opérateur B te se construit 


de façon analogue si le produit est défini par 
t 1 


F(eG (= tr | dE | FIG—H(—nIG(EMan. (5.128 


0 0 
Supposons qu'une fonction & (t) est définie pour tous les t & 0, 


5 0 
qu’elle possède des dérivées de tous les ordres quel que soit t > 0, 
o (t) > 0 lorsque t >> 0 et qu’existent les intégrales 


o* (p)= o (t)e-P' dt, (5.129) 
w* (B) = 2 | o (t) Ko(2 V Bt) t. (5.130) 
0 


Soit L,, l’ensemble de toutes les fonctions f ({) définies sur la 
section [0, c[ et vérifiant la condition 
t 


| [f( 10 ® Œ<o, V>0. (5.131) 
0 


De façon analogue soit B,, l’ensemble de toutes les fonctions f (t) 
définies sur la section {[0, œ[ et vérifiant la condition 


L 


3 
[& [ifmlo(man<o, v>0. (5.132) 
0 


0 


Définissons dans les ensembles L,, et B, l'opérateur U en posant 


= [ete @ dé, fELe (5.133) 
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et 
t 1 


Uf=—T | dE fŒn)w[(t—E) (1—n)le (En) dn, FEB. (5134) 


0 


L'opérateur U possède un inverse U-!. En effet, de la condition 
Uf = 0 qui est valable pour Vt>> 0, il résulte qu’au moins l’une 
des fonctions f (t) w (t) ou w (t) est presque partout nulle. Or © (t) > 
> 0, donc f (t) — 0 presque partout. On remarquera que pour que 
l'opérateur inverse Ut existe, il suffit que © (t) 0 presque par- 
tout sur la section [0, cf. 

Soient UL,, l’ensemble de toutes les fonctions F (t), t > 0, 
représentables sous la forme 


F(t) = Uf, f€eL,; (5.135) 
UB,, l’ensemble de toutes les fonctions F (t), t => 0, représentables 
sous la forme 

F(t) = Uf, fEeB,. (5.136) 


On se servira des DA PIODS 


h (t) — OI [re-DrBoe-0b& (5.137) 


et 
t 


1 
hO= TT | & | FLE —E (1m x 


X g (En)@[(t—E) (1 —n)]® (En) an. (5.138) 
De SL il s'ensuit que Sc tous lest > 0,ona 


[neLE&< | |f @ le © & | | 8) le (@®) dE. 


Donc, siféL,etgeL,,, : fonction À (t) définie par (5.137) appar- 
tient également à L,.. 


S'agissant de l’égalité (5.138), pour tous les t > 0, on a une rela- 
tion analogue (cf. [52]): 


4 E 
Î[Æ [12 6m) 1e (m dn< 


0 0 


nt 


t 
<{ À Lente an Ÿ & À 1eme en dn 
0 ( 


Donc, sif£EB,etge Bo alors la Le h (t) définie par (5.138) 
appartient également à B,. Par conséquent, les ensembles L,, et B, 


9e 
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sur lesquels l’addition et la multiplication, définies respectivement 
par (5.137) et (5.138), sont comprises au sens ordinaire, forment des 
anneaux commutatifs. De là il résulte que les ensembles UL, et 
UB, sont également des anneaux si l’addition est prise au sens ordi- 
paire et la multiplication M à l’aide des égalités 


F (+ G()=U TE | Fe-D6bee-D0@æ (5.139) 


où F(t)EUL,, G(t)EUL,, ou bien 


t 


1 
FE) + G(9 UT | & | F [(t—È) (1 — n)1 


X G(Ën)w[(t—Ë) (1 —n)]o (En) dn, (5.140) 
où F(H)EUB,, Gt) E UB,. 
En effet, en mettant (5.137) et (5.138) sous la forme 


kh = Jg, 
(5.139) et (5.140) deviennent 
F (t)=G (t) = UT'FG. (5.141) 
Par hypothèse, F = Uf et G = Ug, donc 
F(t)=G(t) = Ur! (Uf-Up). (5.142) 
Prouvons maintenant que dans les anneaux L, et B, on a 
U (fg) = Uf-g = f-Ug. (5.143) 


Démontrons tout d’abord (5.143) pour l'anneau Z,. Soit 


RO re D g (Bu (t—Do(E) dE fe. 
HR : 
h (E)w(t— ©) o (®) dt — 


IDE LITE CN ECEE 


g (v) © (v) du ot-5 560 o (E— v) dé = 
0 
t t— 


= À g () © () | & (£—v— u) f (u) w (u) du. 


0 0 


ue ©‘ ve 
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Mais 
Uf=-T | fOo(t-DoŒ)Æ=F(E, (5.144) 
0 
donc 
{ 


frDo(t—-D0® &= (ete) F(t—ve(t—v) dr. 


0 0 
d'où 
Uh= TT h (E)o (t—E) o () dé = 
= | Fe-Drte6-netæ-ur g. 


Le meilleur moyen de Ar (5.143) pour l’anneau B, est de 
se servir de ue (cf. [5] 


© (zy) k (zy) — Ja [rem eter-8 (1 —1n)] © {(zy —E) (1 — n)] X 


0 

V 

X w (Er) dn= | FE gl(z— 9 (mio (Em ol(z—9 (y—n)]ddn. 

0 
D'où il vient 


X ]|}/ X Y 
| [An o(moltz—E5 (4m dan= Î Lot(r—E5) (y—miéan x 
0 0 
+ 0 0 
x | À (u0) o (uv) 8 [(E—u) (n—v)] à 1(E—u) (n—v)] du dv = 
0 oO 


X y 
= | À 7 En) G1(z— 5) (y—n)Jo (En) o [(z—8) (y — n)] dE dn. 
0 
où 
x 1! 
C= ser | g (En) © I(z — 8) (y — n)] © (En) dE dn = 
XV 


= [a | g (En) © [(zy— 5) (1 — n)] © En) én. 


0 
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Posant | 
H (xy) mn ( ÉTTTTEE (y— n)] (En) dé dn — 
# : 
= | dE h (En) w [(zy —#) (1 — n)] & (En) dn, 
on trouve . 
HN = | [ Em EL — D (4m x 


X &[(z—E) (y—n)] (En) dé dn — 


Î dé | f En) G L(zy —®) (1 —n)] © [(zy—8) (1— n)] © En) dn 


w = 


En remplaçant dans les dernières égalités zy par la variable f, on 
obtient 


{ 1 
H()=—T | 28 | Fr) G LED (1 — I 


0 


x 1(E— (1 — lé (En) dn = f6, 


H(t)=— 


ot 


a | k (En) & [(t—#) (1 — n)] © (En) dn = Uh, 


pale 


Ou tre) 


CO=T | & Î & En) w ((t— 8) (1—n)1 (En) dn = Us. 


0 


uw n 


D'où il résulte que 
Uh = U (8) = ÿG = fUg. 
De (5.142) on obtient pour les ensembles UL,, et UB,, 
F (t)°G (t) = UT'FG = Ut (Uf-Ug) = fUg = Ufs. 


Donc, si F et G appartiennent à UL,, ou UB,,, le produit F = G 
appartiendra respectivement à UL,, ou ‘UB.. Ainsi, les ensembles 
UL,, et UB,, sur lesquels le produit est défini par (5.139) et (5. 140) 
sont des anneaux. Ces anneaux ne possèdent pas de diviseurs de zéro, 
ils peuvent donc être étendus à des corps de quotients que nous note- 
rons À (UL4) et R (UB,). Comme la fonction constante f = 1 
appartient à L,, et B,, la fonction 

l 


0= 5 | o(t—u)o(u) du EUL,,; (5.145) 
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et la fonction 
t 


1 
c()=T | & [olt—E (nie (EmdneuBs. (5.146) 
0 


0 


De toute évidence on peut abréger l'écriture des égalités (5.145), 
(5.146), i.e. © (t) = U1. D'où 


C(t)sF(t)—=U"1(U1-F(t))=F(t1=UF(t. 
Donc, dans les anneaux UL,, et UB, le produit par © (t) coïncide 
avec l’action de l’opérateur U, i.e. 
o (t) = U. (5.147) 
L'opérateur inverse appartiendra au corps R (UL,) ou R (UB:) 
et U-1 — 1/0 (t). Désignons spécialement cet opérateur par 


1 
= (5.148) 


On appellera fonction tout opérateur de la forme 


F( 
su = CF (6). 


où F (t) E UL, ou F (t)E UB,. Comme on a toujours F = Uf — 
= (1/Q)f, l'égalité! 


F(?) 
so == (5.149) 
a lieu quelle que soit la fonction F (t) de l’ensemble UL,, ou UB.. 
Dans les corps R (UL.,) et R (UB,) on peut introduire des trans- 
formations intégrales étroitement liées à celles de Laplace et de 
Bessel-Meijer. 
Considérons tout d’abord le corps R (UL,). Un opérateur F/G € 
€ R (UL.) est par définition transformable si existent les intégrales 


F*(p) = | F(t)œ(t)e?t a, 
si L (5.150) 
c°(p}= { Got et dt. | 
0 


À un tel opérateur correspond une fonction de la variable complexe p, 
plus exactement : 
F __ F*(p) 


On montre que l’ensemble de tous les opérateurs transformables 
est un sous-corps de R (UL,) isomorphe au corps de toutes les fonc- 


136 CALCUL OPPRATIONNEL [CH. V 


tions de la forme F* (p)/G* (p) sur lequel l'addition, la multiplication 
et la division sont comprises au sens ordinaire. Ceci étant, à la trans- 
formée de la fonction F (f) € UL,, est associée la fonction de la 
variable complexe 


Ÿ F(é) © (4) e7pt dt 


F(t) = * / 


=" F°)=F(p) (5.152) 


j w(t)e-ptat 
0 


On voit aisément que cette formule est valable pour la fonction 
f (t) E Ls, si seulement existe l'intégrale 


| f(t)w (t)e-?t dt. 
0 
Il est immédiat d'établir que 
Qt 0) 1 
0() * &w(p)ow*(?)  &°(p) 
Donc, 
1 1 : 
Qt 50 — 0 (D). (5.153) 


On obtient des formules analogues pour le corps R (UB.). Sup- 
posons que 


F*(B)=2 F(bo(t)K, (CV à, | 
| (5.154) 


G*(B)=2 


O8 SOS? 9 


G(t)w(t)Ko(2V Bt) dt. | 


On dira que l'opérateur 

+ ER (UBu) 
est transformable si existent les intégrales (5.154). L'ensemble de 
tous les opérateurs transformables forme un corps qui est isomorphe 


à celui de toutes les fonctions de la variable complexe B de la for- 
me F* (B)/G* (B). Cet isomorphisme est défini par la correspondance 


_F(#) __ F*(B) 
T0 = FE. (5.155) 


L'opérateur U — 1/Q est associé à la fonction 


—=2 | o (t) Ko (2 V Bt) dt =w* (B). (5.156) 
0 
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Par ailleurs, si F (t) E UB,, alors 


1 & 
F()=-=7 F*(B). (5.157) 
La formule (5.157) est valable pour toute fonction f (f) de B, si 
existe l'intégrale 
[ f(t)o(t) Ko (2 V Bt) dt. 


0 


Le corps des opérateurs R (UL.) est isomorphe au corps Ÿ des 
opérateurs de Mikusinski [165]. En effet, à tout opérateur 


F 
7 ER(UL) 


on peut associer l'opérateur 
f (t) w (t) 
g (e) © (t) ? 


où les produits f (t)w (t) et g (t) w (t) sont pris au sens habituel. 
Montrons que cette correspondance 

ES À fo 

GT g 8 
définit un isomorphisme. Pour cela, il suffit de montrer que l’anneau 
des fonctions L., est isomorphe à l’anneau Z des fonctions sommables 
sur tout intervalle 0 <<{t << T dans lequel le produit est égal à la 


convolution ordinaire 
t 


n (= (1-5 E 8. 
0 
Si à une fonction j (t)LE L,, on associe la fonction f (t)o (t) € Z, i.e. si 
l'on pose f (t) <> f(t)o (t), il est évident qu’une somme de fonctions 
se transforme en une somme, et le produit 


h=je= TT | 16 -D 8 Eo(t—E) 0 (D 


en le produit des opérateurs jf et g. Ce qui démontre l’isomorphisme 
des corps R (UL,) et . 

On a la même proposition pour le corps R (UB.). Ce corps est 
isomorphe au corps M, introduit dans [52]. L’isomorphisme des 
corps R (UL.) et M, d’une part, et R (UB.) et M3, d'autre part, 
n’amoindrit pas leur rôle dans la théorie du calcul opérationnel. 
Une étude détaillée des opérateurs dans les corps R (UL,.) et R (UB,.) 
peut conduire à de nouveaux résultats en analyse, et notamment en 
théorie des fonctions spéciales. 
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Traitons quelques cas particuliers. Supposons que dans l’anneau 
UL, la fonction & (t) = t"-1T (n), nr > 0. L'ensemble UL,, sera 
composé de toutes les fonctions > oo sous la forme 


F()=U0i= per | fŒ(—EA EI dE. (5.158) 


L'ensemble L, est constitué de toutes les fonctions vérifiant la 
condition 
{ 


Î IF IE" 1 dE < co. 


0 
Supposons que f(t)—t*, où œ>—n. Alors 


UE — To | ES (t—E)"-1 EN GE. (5.159) 
D'où 
œ _,atn Fia+n) ; 
Ut" =t Tatin) : (5.160) 


Donc, toute fonction t8, B + n > 0, appartient à l’ensemble UL.. 
Lorsque & = 0, on déduit de (5.160) 


nr 
HUE Fu 
Donc 
4 y FT (a) 
Pi À T (2n) 
et 
r (2n) 
Sè T'(n)47 


De (5.158) il vient 
1  dn 
ni jm À IUf = f (t). 


De là il résulte que sur UL,, l'opérateur Q coïncide avec l'opérateur 
différentiel 
| pm. 
: ant 7 
Donc, sur l’anneau LU, le produit est égal à 


F()G (D) = FT Jre-pc@e-pea (5.161) 
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Comme les fonctions #& et £6 (&« > 0, B>>0) appartiennent à UL,, 
il vient 
4 


æœ ,4B _ 1 dn 1 a+n-1 pP+n-1 
L'si mi Tam T'(n) [ (t —E) ë dE. 
0 
D'où 
1 
1 - dn 9n— n— = 
RE = es RSS 1 | (A—u)2t"- 1 uBn-1 qu, 
0 
ou 
œ _,a+8 l(a+n)T (B+n) 
et = l'(a+B+n)T (nr) : 
Donc 
EE 5.169 
_—— _ . 
F(a+n) T(B+n) C(a+B+n)T(n) * ( ) 
Cette formule nous permet de calculer facilement 
1 —_ D (2) n 
D'où il résulte 
- __ (5.164) 


G —a = l(n) à F(n+kn) © 


Si dans l'anneau UB, on pose 


ni 
« (£) TER) ‘ n>0, 


l’ensemble UB, sera composé de toutes les fonctions représentables 
sous la forme 


t 1 
FO SU = pr | ("Et dE | 5 (En) (1 — m1 nt an. 
: ° (5.165) 


L'ensemble B,, sera composé de toutes les fonctions vérifiant la condi- 
ti 
ion | . : 
Î À | 1 f6m1 nt dn < co. (5.166) 
0 0 


Si f(t) =t", alors pour &œ+n>>0 la fonction {€ B, et 
i 
| TN n— = 
= (t— EL ESt ‘& | (—n)"tnétr-1 dn. 


0 


1440 CALCUL OPÉRATIONNEL (CH. V 


Donc 


Ut it" Far (5.167) 


Pour &æ—=0 on obtient 


D'où 


et 


___ T2(2n) 


De (5.165) on déduit que sur l’ensemble UB, l'opérateur Q 
coïncide avec l’opérateur différentiel 


(a) e (ge 


Donc, sur l’anneau UB, le produit est égal à 
1 jdn n/din 1 
F(t)}»-G(t)=— ini (5) t (+) T3 (n) X 
t 1 
x À dB Ê FI(—E) (1—n)1 6 (En) @— EE (4 mt net an. 
0 


(1 


(5.169) 
En particulier, lorsque F° = #, G =- æ, a+n>0, B+n>0, il vient 
CE (5.170) 
Pi(a+n) L[i(B+n)  T(a+B+n) l(n) ° ° 
L'égalité 
Le T3 (7) 2” 
Q  T2(2n) 
et (5.170) entraînent 
14 _  T2(n) ï 
Donc 
Q Q rs (n) tknah 
Q—a l3(n +kn) LE 


Par des méthodes standard on établirait de nombreuses formules 
analogues à (5.164) et (5.172), i.e. on calculerait des opérateurs du 


$ 17] SUR UNE GENERALISATION DU CALCUL OPÉRATIONNEL 141 


type f (QG). Souvent on simplifie singulièrement ce problème en 
utilisant les transformations intégrales 


f(p) = | 10 terrt at 
0 


et 


ss 


f(B = | f ()2-1K, (2 V Bi) dt. 


Les formules opérationnelles obtenues permettent en particulier 
de résoudre des équations différentielles de la forme 


Oz (t)+aTiz(t)+...+a,z(t)=f(t), (5.173) 
où (2 est diriger 


1 = 1 pd \n nf d \"ns 
nt rt OÙ ni (+) ; (+) FE 


Certains cas particuliers de l'opérateur © sont intéressants. Par 
exemple, dans R (UB.) l'opérateur (2 se transforme pour nr = 1 en 
l'opérateur de Bessel 


et pour z — 2 en l'opérateur 
sr (gr) ei ie 
(nd 
me rs Fe +8 FF TT OR Sr E 


Observons que mul Pal par l'opérateur (5.174) 


(5.174) 


1 d 
admet pour solutions 
B ei (at) — SEA TS 
B er (at) = ber £ V'at) | 
x ei (at) = kei e V'at) 
ker (2 V'at) 


x er (at) — … 
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Ces fonctions jouent un rôle important dans le calcul opérationnel 
du cas particulier examiné de l'opérateur . Par exemple, de (5.172) 
il vient que 


5 =B ei (Van). 


En procédant comme pour l’opérateur S — 1/0 (ft), on pourrait 
étudier d’autres opérateurs, par exemple, l'opérateur 1/4%, &« >> 0. 
À un tel opérateur est associé un certain opérateur différentiel du 
type ©. Ainsi dans le corps R (UL.), où w (t) = t"-1/T (n), l’opéra- 
teur 1/1”, m étant un entier positif, est associé à l'opérateur diffé- 
rentiel 
1 d'\m ns [ (n) 
tn-1 (5) l'(ntm)° 
Dans le corps R (UB,), où w(t) = t"-1/T?(n), l’opérateur 1/£”, 
m étant un entier positif, est associé à l'opérateur différentiel 

1 d \m m/f d\m nu  T?(n) 
7 (PC 

En conclusion, on remarquera que la théorie du calcul opération- 
nel de l'opérateur S peut être développée bien plus amplement que 
nous venons de le faire, en particulier dans les problèmes relatifs 
aux fonctions opérationnelles, aux séries opérationnelles, à la réso- 
lution des équations différentielles et à la théorie des fonctions 
spéciales. 
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TABLES DES FORMULES 


CHAPITRE VI 


LISTE DES NOTATIONS DES FONCTIONS SPÉCIALES 
ET DE CERTAINES CONSTANTES 


6.1 arccos zx — _ Intz+y/z—1 


6.2 arcsin = In (iz+ Y 1—727° 


4 —ir 


6.3 arctg z=+ In 


1+iz 
ui z—i 
6.3a arccotg 2=7 In pau 
6.4 Archzæz— In (xz-t V z° —1) 
6.5 Arshz=in(z+Wz+1) 
1, 1+ 
6.6 Arthz=- In 1 
6.6a Arcth 2 In 2 T+ 
6.7 (@)n=a(a+1) ... (a+n—1); n=1,2, ... 
6.8 (a)o = 1 
_ F(a+w) 
6.9 Gh=e 
a\ _ F(a+1) 
10 | G=ronre-sen 
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me 
6.11 | B(k)— | cos Pap Jo 
û VT=F sing 
{ 
6.12 | B(z, n= | Ex-1 ({—EJU-1 Œ — 
0 


T (2) T (y) 
T(z+y) 


6.13 bei, (z)=Im{J, (it V i x)], où Im (z) désigne la partie imaginaire du 
nombre complexe z 


6.14 ber, (z)=Re[J,(i V à z)], où Re(z) désigne ls partie réelle du nom- 
bre complexe z 


_ ce —\n an 
6.15 bei (z)=beiz=Im{/o(Vi = Ÿ EÛTE (5) : 
n=0 


_ = A 
6.16 bero (x) =berz—Re[/o(W z)]= }, co(s)" 
n=0 


6.17 = —[" (1)= —14(1)—0,577215665 … 
T (nr + 2v) 1  1—7Z7 
: V 1 — _— D pe res 
6.18 Ch (z) = D'(r+1)T (2) 2F1 (r+2v, n,V+T ? 2 |] 
6.19 C(= = | = 608 u du 
| V 2a L u 
x | x 
| e7iu { (2) 
6.20 C(z)—iS (x) = Va du — 3 [ H=+ (u) du 
0 0 “ 
oO 
6.21 Ceon (Z; 9) = >. AG) cos 2kz est la fonction de Matieu 
k=0 


6.22 Clon+! (z, q) = > AG D cos (2k+1)3 
ke=0 
6.23 Cesn (Z- 9) = Cean (iz, 9) 
6.24 Ceon+s1 (2: 9) = Cean+s (iz, 9) 
C chu—1 
6.25 chi @=in yz+ | à 
0 


du 


e* + ex 


6.26 chr— 5 


CH. VI) 


6.27 cos ie 
2 
œ < 
6.28 Gi()=—| St du= in yz— | 
x 0 
6.29 Ci (x) = —ci (x) 
cs 
; " sin? @ cos? p dp 
6.30 cw=| SRE FT 
0 ({—X2sin2 p)° 
7 
6.31 | D(k)— (see 
| V 1—ksint 
Le 
6.32 Dn(z)=e * Hen(z); n—0,1,2, 
Div A ” 
6.33 | D,(z)—2% ?z 2W, . (+) 
A 
6.34 e—2,718281828 ..… 
Le 
z" 
6.35 ex=expz= Ÿ nl 
n=0 
_ 
6.36 E (= | V  1—k2 sin" q dy 
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